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Preface 

This  report  is  a  product  of  the  U.S.  Army  Corps  of  Engineers’  Risk  Analysis  for  Water  Resources 
Investments  Research  Program.  The  program  is  managed  by  the  Institute  for  Water  Resources,  which  is 
a  unit  of  the  Water  Resources  Support  Center.  The  report  was  prepared  to  fulfill  part  of  several  work  units 
in  the  research  program.  These  work  units  focused  on  developing  and  applying  the  concepts  of  risk 
preference  and  risk  communication  to  water  resources  issues.  The  report  conforms  to  the  basic  planning 
model  and  to  the  risk  and  uncertainty  analysis  recommendations  presented  in  Economic  and  Environmental 
Principles  and  Guidelines  for  Water  Related  Land  Resources  Implementation  Studies  (P&G). 

The  risk  analysis  framework  encompasses  the  four  basic  steps  in  dealing  with  any  risk: 
characterization,  qualification,  evaluation,  and  management.  The  purpose  of  conducting  these  analyses  is 
to  provide  additional  information  to  both  Federal  and  non-Federal  partners  on  the  engineering  and  economic 
performance  of  alternative  investments  that  address  water  resources  problems.  The  goal  is  to  produce  better 
informed  decisions  and  to  foster  the  development  of  the  idea  of  rational  joint  consent  by  all  parties  to  an 
investment  decision. 

This  report,  entitled  Risk-Based  Evaluation  of  Flood  Warning  and  Preparedness  Systems,  represents 
a  synthesis  and  elaboration  of  three  earlier  technical  reports'  to  the  Institute  for  Water  Resources  prepared 
by  Environmental  Systems  Modeling,  Inc.  The  results  presented  here  have  as  a  unifying  theme  that  design 
and  evaluation  of  structural  and  nonstructural  measures  for  flood  mitigation,  including  flood  warning  and 
preparedness  systems,  is  an  integrative,  holistic  process  that  requires  an  understanding  of  the  contribution 
each  type  of  measure  makes  to  the  performance  of  the  overall  system.  The  models  rely  on  concepts  of 
multiobjective  decisionmaking,  tradeoff  analysis,  and  the  risk  of  extreme  events.  This  report  is  divided  into 
OVERVIEW  and  TECHNICAL  sections.  Each  of  the  four  OVERVIEW  sections  summarizes  in  a 
nontechnical  style  a  methodology  developed  for  the  integration  of  flood  warning  and  preparedness  systems 
into  the  design  and  evaluation  process.  The  four  methodologies  are  (1)  integration  of  structural  measures 
and  flood  warning/preparedness  systems,  (2)  multiobjective  decision-tree  analysis,  (3)  performance 
characteristics  of  a  flood  warning  system,  and  (4)  selection  of  optimal  flood  warning  threshold.  Each 
OVERVIEW  section  describes  the  main  features  of  the  model,  case  study,  or  example.  The  four 
TECHNICAL  sections  correspond  to  the  sections  of  the  OVERVIEW  and  contain  the  mathematical  details 
that  would  be  needed  in  an  application  of  the  methodologies.  In  addition  to  being  a  consultant  for  this 
report.  Prof.  Roman  Krzysztofowicz  is  the  sole  author  of  the  TECHNICAL  section  of  Part  3-Performance 
Characteristics  of  a  Flood  Warning  System;  the  OVERVIEW  section  of  Part  3  is  excerpted  and  edited  from 
the  same  TECHNICAL  section.  The  contribution  and  description  of  case-study  data  in  Section  4— Selection 
of  Optimal  Flood  Warning  Threshold— also  is  adopted  from  work  of  Krzysztofowicz. 


'Performance  Characteristics  of  a  Flood  Warning  System  and  Selection  of  Optimal  Warning 
Threshold  (September  1990);  Case  Studies  in  Selecting  Optimal  Flood  Warning  Threshold  (March 
1992);  and  Integration  of  Structural  Measures  and  Flood  Warning  Systems  for  Flood  Damage 
Reduction  (March  1992) 
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Mr.  Stuart  A.  Davis,  of  the  Institute  for  Water  Resources’  Technical  Analysis  and  research 
Division,  provided  technical  review  for  this  report.  Dr.  Eugene  Z.  Stakhiv,  Chief  of  the  Policy  and  Special 
Studies  Division,  and  Dr.  David  A.  Moser,  of  the  Technical  Analysis  and  Research  Division,  reviewed 
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Introduction 


Each  of  the  four  methodologies  developed  in  this  report  contributes  an  important  dimension  to  risk-based 
evaluation  for  flood  damage  reduction  -  which  is  incomplete  without  accounting  for  both 

structural  and  nLstructural  measures.  The  unifying  theme  of  these  results  is  that  the  design  and  evaluation 
of  structural  and  nonstructural  measures  for  flood  mitigation,  including  flood  warning  and  Preparedness 
systems  is  an  integrative,  holistic  process  that  eventually  must  build  on  an  understamding  of  the  contribution 
of  each  type  of  measure  to  the  performance  of  the  overall  system.  Furthermore,  the  design  of  flood  mitigation 
is  tied  to^multiple  objectives  of  minimizing  cost  and  risk  and  maximizing  performance.  Consideration  of  the 
risk  of  extreme  events  is  an  essential  element  in  the  evaluation  of  design  tradeoffs. 

The  four  methodologies  developed  here  for  the  modeling  and  evaluation  of  flood  warning  and  preparedness 
systems  are; 

(1)  Integration  of  structural  measures  and  flood  waming/preparedness  systems, 

(2)  Multiobjective  decision-tree  analysis, 

(3)  Performance  characteristics  of  a  flood  warning  system,  and 

(4)  Selection  of  optimal  flood  warning  threshold. 

The  assumptions,  main  functions,  and  limitations  of  the  four  methodologies  are  summarized  in  Table  1 . 

Multiple  Objectives 

The  single-objective  models  that  had  been  advanced  in  the  fifties,  sixties,  and  seventies  are  today  considered 
by  many  to  be  unrealistic,  too  restrictive,  and  often  inadequate  for  most  real-world  complex  problems.  The 
proliferation  of  books,  articles,  and  conferences  and  courses  during  the  last  decade  or  two  on  what  has  come 
to  be  known  as  multiple-criteria  decisionmaking  (MCDM)  is  a  vivid  indication  of  this  somber  realization  and 
of  the  maturation  of  the  field  of  decisionmaking  [see  Chankong  and  Haimes  1983].  In  particular  an  optinium 
derived  from  a  single-objective  mathematical  model,  including  that  which  is  derived  from  a  decision  tree 
often  may  be  far  from  representing  reality  -  thereby  misleading  the  analyst(s)  as  well  as  decisionmaker(s). 
Fundamentally,  most  complex  problems  involve,  among  other  things,  Ae  minimization  of  costj  the 
maximization  of  benefits  (not  necessarily  in  monetary  values),  and  the  minimization  of  risks  of  various  km  s. 
For  example,  decision  trees,  which  are  a  powerful  mechanism  for  the  analysis  of  complex  problems,  can 
better  serve  both  the  analysts  and  the  decisionmakers  when  they  are  extended  to  deal  with  the  above  multiple 

objectives. 

Impact  Analysis 

On  a  long-term  basis,  managers  and  other  decisionmakers  are  often  rewarded  not  because  they  have  made 
many  optimal  decisions  in  their  tenure;  rather,  they  are  honored  and  thanked  for  avoiding  adverse  and 
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Table  1.  Assumptions,  Main  Functions,  and  Limitations  of  the  Four  Methodologies 


Assumptions 


Main  Functions 


Limitations 


Integration  of  structural 
measures  and flood 
warning/ 

preparedness  systems 

Knowledge  of  flood 
frequency,  discharge, 
stage,  and  damage 
relationships  for  various 
combinations  of  structural 
and  flood 

waming/preparedness 

systems. 

Determine  the  optimal 
design  options  from  among 
alternative  combinations  of 
structural  and  flood 
waming/preparedness 
measures  in  a 
multiobjective  framework, 
including  cost,  the 
expected  flood  loss,  and 
risk  of  extreme  floods  . 

No  operational  issues 
associated  with 
waming/preparedness 
systems  and  stmctural 
measures  are  considered. 

Multiobjective  decision- 
tree  analysis 

Knowledge  of  the 
probabilities  for  the 
underlying  distributions  of 
water  level.  Knowledge  of 
severity  of  loss  with 
alternative  decisions  at 
various  time  stages. 

Determine  the  optimal 
sequential  decisions  in  an 
individual  flood  event 
based  on  the  observation  of 
water  stage. 

No  flood  forecast  is  taken 
into  consideration. 

Performance 
characteristics  of  a  flood 
warning  system 

Knowledge  of  the  joint 
probability  description  of 
flood  forecast  and  actual 
flood  crest. 

Provide  an  evaluation 
model  of  the  performance 
of  a  flood  forecast  system. 

In  particular,  the  ROC 
curve  characterizes  the 
tradeoff  between  the 
probabilities  of  detection 
and  false  warning. 

Interactions  between 
successive  flood  events 
through  the  dynamics  of 
the  community  response 
fraction  are  not  taken  into 
account. 

Selection  of  optimal  flood 
warning  threshold 

Knowledge  of  the  joint 
probability  description  of 
flood  forecast  and  actual 
flood  crest.  Knowledge  of 
the  loss  to  the  community 
associated  with  flood  stage. 
Knowledge  of  the 
dynamics  of  the 
community  response 
fraction. 

Find  the  optimal  threshold 
level  at  which  to  issue  a 
flood  warning  in  order  to 
balance  the  desire  for  high 
present-flood-loss 
reduction  with  the 
possibility  of  high  future 
flood  loss  being  inevitable. 

The  derived  optimal 
threshold  may  not  be 
stationary;  i.e.,  the  optimal 
threshold  may  vary  in 
different  flood  events  even 
if  the  community  response 
fraction  is  the  same. 
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catastrophic  consequences.  If  one  accepts  this  premise,  then  the  role  of  impact  analysis  -  studying  and 
investigating  the  consequences  of  present  decisions  on  future  policy  options  --  might  be  as  important,  if  not 
actually  more  so,  than  generating  an  optimum  for  a  single-objective  model  or  identifying  a  Paret^ptimum 
set  (a  Pareto-optimum,  or  non-inferior,  alternative  cannot  be  improved  in  any  one  objective  without  seeing 
a  corresponding  loss  with  respect  to  one  or  more  other  objectives)  for  a  multiobjective  model.  Certainly,  when 
the  ability  to  generate  both  is  present,  having  an  appropriate  Pareto-optimum  set  and  knowing  the  impact  of 
each  Pareto-optimum  on  future  policy  options  should  enhance  the  overall  decisionmaking  process. 

The  Risk  of  Extreme  and  Catastrophic  Events 

Risk  which  is  a  measure  of  the  probability  and  severity  of  adverse  effects,  has  until  recently  been  commonly 
quantified  via  the  expected-value  formula.  This  formula  essentially  precommensurates  events  of  low 
frequency  and  high  damage  with  events  of  high  frequency  and  low  damage.  Although  learned  students  of  risk 
analysis  recognize  the  disparity  between  the  above  fallacious  representation  of  extreme  and  catastrophic 
events  and  the  perception  of  these  events  by  individuals  or  the  public  at  large,  many  continue  to  use  this 
approach.  The  trend,  however,  is  moving  toward  the  conditional-expected-value  approach,  where  extreme 
and  catastrophic  events  are  partitioned,  isolated,  quantified  in  terms  of  conditional  expectation  (e.g.,  using 
concepts  from  the  statistics  of  extremes),  and  then  evaluated  along  with  the  common  expected  value  of  risk 
or  damage  [Asbeck  and  Haimes  1984;  Haimes  1988;  Karlsson  and  Haimes  1988]. 

The  partitioned  multiobjective  risk  method  (PMRM)  developed  by  Asbeck  and  Haimes  [1984]  separates 
extreme  events  from  other  noncatastrophic  events,  and  thus  provides  the  decisionmaker(s)  with  additional 
valuable  and  useful  information.  In  addition  to  using  the  traditional  expected  value,  the  PMRM  generates  a 
number  of  conditional  expected-value  functions,  termed  here  risk  functions,  which  represent  the  risk,  given 
that  the  damage  falls  within  specific  probability  ranges  (or  damage  ranges). 

Combining  either  a  conditional  expected  risk  function  or  the  unconditional  expected  risk  function  with  the 
cost  objective  function  creates  a  set  of  multiobjective  optimization  problems  in  which  the  tradeoffs  between 
cost  and  the  risk  arising  from  the  various  ranges  of  damage  are  analyzed.  This  formulation  offers  more 
information  about  the  probabilistic  behavior  of  the  problem  than  the  single  multiobjective  formulation  that 
minimizes  only  the  cost  and  the  expected  damage.  The  tradeoffs  between  the  cost  function  and  any  risk 
function  allow  decisionmakers  to  consider  the  marginal  cost  of  a  small  reduction  in  the  risk  objective,  given 
a  particular  level  of  risk  assurance  for  each  of  the  partitioned  risk  regions,  and  given  the  unconditional  risk 

function. 

Flood  Forecasting  and  Warning/Preparedness  Systems 

Flood  damage  reduction  can  be  provided  by  either  structural  or  nonstructoal  measures  or  a 
combination  of  both.  Structural  flood  damage  reduction  measures,  such  as  an  increase  in  dam  height,  affect 
the  flood-frequency  relationship.  Nonstructural  measures,  such  as  a  flood  waming/preparedness  system,  do 
not  have  an  impact  on  the  flood-frequency  relationship;  however,  they  modify  the  flood-damage  relationship. 

The  benefits  of  flood  forecasts  have  been  studied  and  systems  approaches  to  flood  forecasting  have 
been  pursued  by  many  research  scholars  for  more  than  twenty  years  [NACOA  1972;  Bhavnagri  and 
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Bugiiarello  1965;  Bock  and  Hendrick  1966;  Day  and  Lee  1976;  Lee  et  al.  1975;  Sniedovich  et  al.  1974; 
Sniedovich  and  Davis  1977].  Curtis  and  Schaake  [1988]  evaluated  flood  warning  benefits  both  on  a  national 
(or  regional)  scale  and  on  a  specific  site  problem.  Prediction  models  for  loss  of  life  from  floods  were  studied 
by  Lee  et  al.  [1986]  and  Shabman  [1987].  Barrett  et  al.  [1988]  developed  categories  for  flood  warning  systems 
based  on  types  of  flood  forecasting  systems  and  flood  response  systems. 

Predicting  the  future  behavior  of  a  time-dependent  random  variable  is  a  major  research  task  in  the 
theory  and  applications  of  stochastic  processes.  Critical  events  occur  when  the  level  of  the  random  variable 
crosses  a  given  high  level  (e.g.,  flooding  level).  An  alarm  is  set  off  when  the  random  variable  exceeds  a 
specified  threshold  level.  An  alarm  system  is  considered  optimal  if  it  detects  catastrophes  with  an  acceptable 
level  of  probability  and  at  the  same  time  yields  a  minimum  expected  number  of  false  alarms  [Lindgren  1 979, 
1980,  1985;  de  Mar^  1980].  The  paper  by  de  Mar6  [1980]  indicates  that  when  Judging  the  performance  of 
an  alarm  system,  it  is  not  very  interesting  to  know,  in  the  mean,  how  close  the  prediction  is  to  the  actual 
process;  however,  it  is  important  for  a  system  to  be  able  to  detect  catastrophes  without  causing  too  many  false 
alarms. 


In  a  series  of  papers,  Krzysztofowicz  and  his  colleagues  [Alexandridis  and  Krzysztofowicz  1985; 
Ferrell  and  Krzysztofowicz  1983;  Krzysztofowicz  1983a,  b;  1985;  Krzysztofowicz  and  Davis  1983a,  b,  c, 
d;  1984]  conceptualized  the  flood  forecast-response  process  in  the  form  of  a  total  system.  This  system  is 
defined  as  a  cascade  coupling  of  two  components:  (1)  the  forecasting  system,  which  includes  data  collection, 
flood  forecasting,  and  forecast  dissemination;  and  (2)  the  response  system,  which  encompasses 
decisionmaking  and  action  implementation.  Based  on  the  above  mathematical  description  of  the  physical 
flood  forecast-response  process,  Krzysztofowicz  and  his  colleagues  establish  performance  measures  of  flood 
warning  systems. 

Pate-Comell  [1986]  presents  a  method  for  assessing  the  performance  of  the  forecasting  system  and 
human  response,  given  the  memory  that  people  have  kept  on  the  quality  of  previous  alerts.  The  tradeoff 
between  the  rate  of  false  alerts  and  the  length  of  the  lead  time  is  studied  to  account  for  the  long-term  effects 
of  "crying  wolf"  An  explicit  formulation  of  benefits  from  warning  systems  is  derived  under  the  above 
considerations. 

Toward  Implementation  of  the  Methodologies 

An  immediate  and  most  worthwhile  challenge  is  the  refinement  of  the  four  methodologies  of  this 
report  for  the  operational  setting.  For  instance,  a  decision  support  system  for  the  risk-based  evaluation  of 
flood  warning  systems  might  be  developed  to  integrate  these  methodologies  in  a  framework  consistent  with 
Corps  of  Engineers  planning  procedures  [HEC  1988]. 

Organization  of  the  Report 

The  body  of  this  report  has  two  major  types  of  subdivisions:  the  OVERVIEW  and  the  TECHNICAL 
sections.  Each  of  the  four  sections  subtitled  OVERVIEW  summarizes  in  a  nontechnical  style  a  methodology 
developed  for  the  integration  of  flood  warning  systems  into  the  design  and  evaluation  process.  Each 
OVERVIEW  section  describes  the  main  features  of  the  model,  case  study,  or  example.  The  four 
TECHNICAL  sections  correspond  to  the  sections  of  the  OVERVIEW  and  contain  the  mathematical  details 
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that  would  be  needed  in  an  application  of  the  methodologies.  The  OVERVIEW’S  present  an  excerpted  group 
of  the  figures  and  tables  used  in  the  TECHNICAL  sections. 


Risk-Based  Evaluation  of 
Flood  Warnine  and  Preparedness  Systems 


Parti 

Integration  of  Structural  Measures 
And  Flood  Warning  Systems: 

Technical 

Introduction 

In  most  cases,  the  maximum  flood  loss  reduction  can  be  only  achieved  through  an  optimal 
combination  of  both  structural  and  nonstructural  flood  damage  reduction  measures,  since  the  adoption  of 
integrated  measures  will  certainly  enlarge  the  feasible  region  of  flood  damage  reduction  measures  when 
compared  with  situations  where  only  structural  or  nonstructural  measures  alone  are  considered.  Structoal 
measures  include  the  construction  of  reservoirs,  leVees,  and  flood  walls.  Nonstructural  measures  include 
floodplain  land  use  planning,  flood  insurance,  flood  warning  systems,  floodproofing,  and  permanent 
relocation  Various  flood  damage  reduction  measures  prevent  inundation  of  the  floodplain  m  different  ways 
and  have  different  impacts  on  the  flood  damage-frequency  relationship.  A  structural  measure,  such  as  an 
increase  in  reservoir  height,  affects  the  frequency-discharge  relationship;  levees  and  flood  walls  confine  the 
discharge  within  certain  channels,  thus  changing  the  relationship  between  discharge  and  elevation;  most 
nonstructural  measures,  such  as  a  flood  warning  system,  modify  the  stage-damage  relationship. 

The  idea  of  combining  both  structural  and  nonstructural  measures  in  flood  damage  reduction  is  not 
new.  Various  research  results  have  been  reported  that  combine  structural  measures  with  nonstoctural 
measures,  such  as  zoning,  floodproofing,  and  flood  insurance.  Readers  can  refer  to  Thampapillai  and 
Musgrave  [1985],  which  provides  a  comprehensive  survey  in  reviewing  integrated  structural  and 
nonstructural  measures  in  flood  damage  mitigation.  To  date,  however,  no  other  research  work  on  combining 
structural  measures  with  flood  warning  systems  has  appeared  in  the  literature. 


Issues  of  both  design  and  operation  are  involved  in  structural  measures  as  well  as  nonstructural  ones. 
Building  a  reservoir  is  a  structural  measure  in  flood  damage  reduction.  Determination  of  the  height  of  the 
reservoir  is  a  design  issue,  while  determination  of  the  amount  of  the  release  on  a  monthly  or  daily  basis  is  an 
operational  issue.  Installing  a  flood  warning  system  is  a  nonstructural  measure  in  flood  damage  reduction. 
Determination  of  an  acceptable  reliability  of  a  warning  system  is  a  design  issue  with  a  consideration  of  the 
system  cost,  while  determination  of  the  flood  warning  threshold  for  various  flood  events  is  an  operational 
issue.  It  is  important  to  note  that  operational  issues  can  only  be  addressed  in  a  framework  of  dynamic 
optimization.  For  example,  different  levels  of  flood  warning  thresholds  will  cause  different  probabilities  of 
missed  forecast  and  false  alarm,  thus  affecting  the  fraction  of  the  community's  future  response.  In  this  part, 
we  consider  only  the  design  options  for  both  structural  measures  and  flood  warning  systems;  thus,  building 
on  and  extending  the  existing  methodology  of  computing  flood  loss  for  a  given  structural  measure  developed 
by  the  Army  Corps  of  Engineers  to  incorporate  flood  warning  systems. 

For  the  computation  of  flood  loss  for  a  given  flood-control  structural  measure,  a  widely-used 
procedure  developed  by  the  Army  Corps  of  Engineers  investigates  the  relationships  between  discharge  vs. 
frequency,  discharge  vs.  elevation,  and  damage  vs.  elevation,  such  that  the  damage-frequency  curve  can  be 
generated  for  an  average  annual  flood  loss.  An  integration  approach  has  been  developed  in  this  report  to 
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combine  the  calculation  of  flood  loss  reduction  through  flood  warning  systems  with  the  calculation  of  flood 
loss  for  a  given  flood-control  structure,  thus  facilitating  the  evaluation  of  combined  structural  measures  and 
flood  warning  systems  in  reducing  flood  loss.  This  new  concept  is  demonstrated  in  an  example  problem. 

Description  of  the  Integrated  Approach 

In  the  procedure  for  computing  flood  damage  that  has  been  developed  by  the  Army  Corps  of 
Engineers,  four  functional  relationships  (or  curves  ~  Fig.  1-1)  are  needed  to  completely  quantify  each 
alternative  of  structural  measures: 

1)  curve  of  frequency  vs.  discharge, 

2)  curve  of  elevation  vs.  discharge, 

3)  curve  of  elevation  vs.  damage,  and 

4)  curve  of  frequency  vs.  damage. 

Note  here  that  the  fourth  curve  can  be  derived  if  the  other  three  are  known.  In  general,  the  relationships  of 
frequency  vs.  discharge,  elevation  vs.  discharge,  and  damage  vs.  elevation  are  constructed  from  real  data  such 
that  the  curve  of  frequency  vs.  damage  can  be  derived  in  order  to  compare  the  expected  flood  damage  for 
structural  measures. 

The  approach  of  discrete  enumeration  of  all  possible  combinations  of  both  structural  and 
nonstructural  measures  is  adopted  in  our  development.  Assume  that  there  are  N  feasible  alternatives  of 
structural  measures  and  M  feasible  designs  of  flood  warning  systems.  Therefore,  there  are  (N  +  1)(M  +  1) 
combinations  of  flood  damage  reduction  alternatives,  and  this  includes  one  do-nothing  option,  N  options 
involving  only  structural  measures,  M  options  involving  only  flood  warning  systems,  and  NM  options 
involving  a  combination  of  both  a  structural  measure  and  a  flood  warning  system. 

In  this  part  we  subscribe  to  a  premise  that  the  introduction  of  a  flood  warning  system  will  not  affect 
the  relationships  between  the  frequency  and  discharge  and  between  the  elevation  and  discharge.  It  will, 
however,  alter  the  curve  of  elevation  vs.  damage,  thus  changing  the  relationship  between  frequency  and 
damage. 


To  evaluate  the  flood  loss  reduction  by  installing  a  flood  warning  system,  the  concept  of  categoiy- 
unit  loss  function  detailed  by  Krzysztofowicz  and  Davis  [1983]  is  adopted.  The  main  modification  is  that  the 
notation  0,  which  was  originally  used  in  Krzysztofowicz  and  Davis  [1983]  as  the  response  degree  of  an 
individual,  is  used  in  this  study  to  represent  the  fraction  of  people  in  a  community  who  respond  to  flood 
warnings. 

The  cost  function  of  evacuation,  C^,  is  assumed  to  be  a  linear  function  of  the  response  fraction 
Ce=MC0  (11) 

where  MC  is  the  maximum  evacuation  cost  for  the  community  when  a  full  response  is  present. 
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Figure  1-1.  Curves  for  Relating  Flood  Frequency,  Discharge,  Elevation,  and  Damage 


Assume  that  the  elevation  of  the  floodplain  zone  under  consideration  is  y  and  the  flood  stage  is  h. 
The  flood  loss  function  without  a  warning  system  is  essentially  given  by  the  curve  of  elevation  vs.  damage 
(Figure  1-lc)  for  each  given  structural  measure.  Alternatively,  the  flood  loss  function  without  a  warning 
system  L„oCan  be  expressed  by 

L,„=MD6(h-y)  (1-2) 

where  MD  is  the  maximum  possible  damage  of  the  community  due  to  a  flood  of  the  highest  magnitude  ^d 
6(h  -  y)  is  the  unit  damage  function  specifying  the  fraction  of  MD  which  occurs  when  the  depth  of  flooding 
is  (h  -  y). 

The  flood  loss  with  a  warning  system  L„is  assumed  to  be  of  the  following  form: 

L„  =  MC  0  +  MD[1  -  0  MR(h  -  y)]6(h  -y)  (1-3) 

where  MR(h  -  y)  is  the  unit  reduction  function  specifying  the  reduction  of  the  maximum  flood  loss  MD 
when  the  depth  of  flooding  is  (h  -  y)  and  full  response  of  the  community  is  made,  i.e.,  0  =  1. 
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In  summary,  the  flood  loss  reduction,  Lrd,  can  be  expressed  as  the  difference  between  L„oand 

Lri,=  0  MD  MR(h  -  y)6(h  -  y)  -  MC  0  (1.4) 

The  value  of  the  maximum  flood  loss  for  a  community,  MD,  and  the  functional  form  of  the  unit 
damage  function  can  be  obtained  from  the  curve  of  elevation  vs.  damage  when  structural  measures  are 
evaluated.  The  only  additional  information  required  to  calculate  the  relationship  of  elevation  vs.  flood  loss 
reduction  through  a  flood  warning  system  is  the  value  of  maximum  evacuation  cost,  MC,  the  value  of 
response  fraction  in  the  community,  0,  and  the  unit  reduction  function,  MR(h  -  y).  The  resulted  curve  of 
elevation  vs.  flood  loss  reduction  can  be  viewed  as  a  function  parametrized  by  the  response  fraction  0.  We 
should  note,  however,  that  the  flood  loss  reduction  Lrb  is  a  linear  function  of  the  response  fraction  0. 

Reducing  the  value  of  damage  in  the  curve  of  damage  vs.  elevation  for  each  structural  measure  by 
LRjj(h;0)  for  each  given  value  of  elevation  h  yields  a  new  relationship  between  elevation  and  damage  when 
a  flood  warning  system  is  introduced.  Setting  0  equal  to  one  yields  a  maximum  achievement  of  flood  loss 
reduction.  Combining  this  new  curve  of  elevation  vs.  damage  with  the  other  two  curves  of  frequency  vs. 
discharge  and  elevation  vs.  discharge  provides  us  with  a  new  relationship  between  frequency  and  damage  for 
a  combined  structural  measure  and  a  flood  warning  system. 

A  recent  report  by  Jack  Faucett  Associates  [1990]  provides  procedures  for  calculating  the  cost  and 
benefits  of  flood  warning  systems,  which  is  useful  in  determining  the  unit  reduction  function  MR(h  -  y)  and 
in  evaluating  the  tradeoff  between  the  cost  and  the  flood  loss. 

Although  the  relationship  of  damage  vs.  frequency  provides  the  most  complete  evaluation  for  each 
flood  damage  reduction  alternative,  it  is  necessaiy  to  compress  information  to  generate  a  risk  measure  when 
various  flood  damage  reduction  alternatives  are  compared.  The  most  commonly  used  risk  measure  is  the 
expected  value  of  the  flood  loss.  Although  the  expected-value  approach  indicates  the  central  tendency  of 
flood  damage  of  each  flood  damage  reduction  alternative,  it  fails  to  separate  the  extreme  catastrophic  flood 
events  from  the  rest.  The  partitioned  multiobjective  risk  method  (PMRM)  [Asbeck  and  Haimes  1984]  adopts 
the  concept  of  conditional  expectation,  which  enables  us  to  isolate,  quantify,  and  evaluate  the  impact  of  each 
flood  damage  reduction  alternative  on  extreme  catastrophic  flood  events. 

Multiobjective  analysis  will  be  performed  in  this  study  to  evaluate  the  various  flood  damage 
reduction  alternatives.  There  are  three  objective  functions.  In  consistency  with  the  notations  used  in  PMRM 
[Asbeck  and  Haimes  1 984],  for  each  flood  damage  reduction  alternative  we  use  f,  to  denote  the  cost,  fj  the 
expected  damage,  and  fjfa)  the  conditional  expectation  of  extreme  floods  whose  return  periods  are  greater 
than  — — .  Both  fj  and  fjCa)  can  be  derived  from  the  curve  of  damage  vs.  frequency  for  each  flood  damage 
reduction“altemative: 

f5=meanof{L}  (1.5) 

and 

f/a)  =  mean  of  {L  |  return  period  of  L  (1 .6) 

A  flood  damage  reduction  option  may  have  different  impacts  on  the  expected  flood  loss  and  the  expected 
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flood  loss  with  floods  whose  return  period  exceed  certain  threshold  level.  This  framework  will  provide  more 
decision  aids  to  determine  the  optimal  flood  damage  reduction  strategy.  The  added  tradeoff  information 
between  the  cost  and  the  expected  extreme  flood  loss  will  explicitly  address  public  concerns  about 
catastrophic  flood  loss. 


Example 

This  section  develops  an  example  to  illustrate  the  integrated  approach  developed  in  Section  11.  The 
following  four  studies  performed  for  or  by  the  Corps  of  Engineers  provided  the  basic  data  for  this  example 
problem: 

1)  Allegheny  River  and  Eldred  Brook,  McKean  County,  Pennsylvania,  1977 

2)  Youghiogheny  River  at  Connellsville,  Pennsylvania— three  reports,  1979, 1980,  and  1985 

3)  South  Branch  Potomac  River  at  Petersburg,  West  Virginia,  1990 

4)  South  Fork  and  South  Branch  Potomac  Rivers  at  Moorefield,  West  Virginia,  1990 

After  reviewing  these  four  sets  of  documents,  the  study  undertaken  for  the  South  Fork  and  South 
Branch  Potomac  Rivers  at  Moorefield,  West  Virginia,  in  1990  was  selected  as  the  basis  for  the  development 
of  the  example. 

T.ocal  Flood  Protection  at  South  Fork  and  South  Branch  Potomac  Rivers  at  Moorefield.  West  Virginja 

The  documents  provided  for  this  study  were  a  reconnaissance  report  dated  September  1 987  and  an 
Integrated  Feasibility  Report  and  Environmental  Impact  Statement  dated  March  1990.  The  latter  consists  of 
a  main  report  and  13  appendices  [U.S.  Army  Corps  of  Engineers  1990]. 

The  following  extracts  from  the  main  report  [U.S.  Army  Corps  of  Engineers  1990]  provide  some 
background  for  this  example: 

The  town  of  Moorefield  in  Hardy  County,  West  Virginia,  is  subject  to  flooding  from  the 
South  Fork  and  South  Branch  Potomac  River.  Serious  floods  have  occurred  in  March  1936, 

June  1949,  and  November  1985  . 

In  response  to  the  flooding  problem,  the  Corps  of  Engineers  and  the  Interstate  Commission 
on  the  Potomac  River  Basin  initiated  a  cost-shared  feasibility  study  in  February  1988  to 
identify  and  evaluate  possible  solutions . 

A  range  of  possible  structural  and  nonstructural  measures  was  examined.  These  measures 
included  levees,  floodwalls,  channel  improvements,  bridge  modification,  and  nonstructural 
alternatives.  The  most  effective  measures  were  combined  into  plans  for  comparison  to  the 
without  project  condition . 
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In  the  example  developed  in  the  following  sections,  two  given  structures  of  flood  damage  reduction 
will  be  investigated.  Plan  1  is  the  zero-cost  plan,  that  is,  the  without-project-condition  alternative.  Plan  4  is 
a  structural  plan  and  includes  levees  and  floodwalls  to  protect  residential  areas,  industrial  plants,  businesses, 
schools,  and  commercial  areas  in  both  North  and  South  Moorefield.  A  detailed  description  of  this  plan  is 
provided  on  page  67  of  the  main  report  [U.S.  Army  Corps  of  Engineers  1990]. 

Development  of  the  four  functional  relationships  for  Plan  1 

For  plan  1,  three  of  the  four  functional  relationships  were  available  from  data  provided  either  in  the 
main  report  or  in  the  appendices.  The  process  for  development  of  the  functional  relationship  was  to  take  the 
original  data  and  perform  a  regression  analysis  in  order  to  obtain  the  functional  relationship.  This  process  was 
carried  out  for  all  the  three  curves  for  which  the  data  were  available.  The  fourth  relationship  was  obtained 
by  use  of  the  other  three  relationships. 

Discharge  versus  Elevation 

The  data  for  this  curve  were  extracted  from  Plate  A17  in  Appendix  B  [U.S.  Army  1990].  The  data 
thus  obtained  are  shown  in  Table  1-1  and  the  corresponding  plot  is  shown  in  Figure  1-2. 


Table  1-1.  Discharge  vs.  Elevation  for  Plan  1 


Discharge 

(1000  cfs) 

Elevation 

(feet) 

10.7 

812.6 

15.0 

814.1 

23.0 

816.3 

31.7 

817.4 

44.3 

818.2 

63.3 

819.1 

96.0 

820.5 

111.0 

821.1 
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Figure  1-2  Discharge  vs.  Elevation  for  Plan  1 


Based  on  this  data,  a  regression  analysis  was  performed  and  the  following  functional  relationship  was 
obtained: 

E  =  1.28767  (—^ - 10.7)°'''“’ -  0.15423  +  812.6  (1-7) 

1000 

where  E  is  the  flood  elevation  in  feet,  and  D  is  the  discharge  in  cfs.  Figure  1-3  shows  the  comparative  plots 
of  the  actual  vs.  fitted  curves,  indicating  that  the  derived  functional  relationship  is  an  adequate  model  of  the 
actual  system. 

Frequency  versus  Discharge 

The  data  for  this  curve  were  obtained  from  Appendix  B  of  the  main  report  [Table  1,  page  9,  U.S 
Army  1990].  The  data  thus  obtained  are  shown  in  Table  1-2  and  the  corresponding  plot  is  shown  in  Figure 
1-4. 
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Figure  1-3.  Actual  vs.  Fitted  Curves  for  Discharge  vs.  Elevation— Plan  1 
Table  1-2.  Frequency  vs.  Discharge  for  Plan  1 


Discharge 

(1000  cfs) 

Flood 

Frequency 

10.5 

0.200 

14.6 

0.100 

22.6 

0.040 

31.2 

0.020 

43.5 

0.010 

62.0 

0.005 

94.0 

0.003 

109.0 

0.002 
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Figure  1-4.  Frequency  vs.  Discharge  for  Plan  1 


Based  on  this  data,  a  regression  analysis  was  was  performed  and  the  following  functional  relationship  was 
obtained: 


D  =  10® 


0.00211 


>  2  0.62411 


+  0.00582 


(1.8) 


where  D  is  the  discharge  in  cfs,  and  F  is  the  flood  frequency.  Figure  1-5  shows  the  comparative  plots  of 
the  actual  vs.  fitted  curves,  indicating  that  the  derived  functional  relationship  is  an  adequate  model  of  the 
actual  system. 


Frequency  versus  Damage 

The  data  for  this  curve  were  obtained  from  Appendix  I  of  the  main  report  [Table  1-23,  page  1-35, 
U.S.  Army  Corps  of  Engineers  1990].  The  data  thus  obtained  are  shown  in  Table  1-3  and  the  corresponding 
plot  is  shown  in  Figure  1-6. 
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Figure  1-5.  Actual  vs.  Fitted  Curves  for  Frequency  vs.  Discharge— Plan  1 
Table  1-3.  Frequency  vs.  Damage  for  Plan  1 


Damage 

Flood 

($  million) 

Frequency 

0.000 

0.100 

0.000 

0.050 

0.712 

0.040 

3.354 

0.020 

4.864 

0.013 

5.430 

0.010 

6.601 

0.005 

7.141 

0.003 

7.682 

0.002 
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Figure  1-6.  Frequency  vs.  Damage  for  Plan  1 

Based  on  this  data,  a  regression  analysis  was  performed  and  the  following  functional  relationship  was 
obtained. 

L  =  -56.65856  +  9.02988  (1-9) 

where  L  is  the  flood  damage  in  $  million,  and  F  is  the  flood  frequency.  Figure  1-7  shows  the  comparative 
plots  of  the  actual  vs  fitted  curves,  indicating  that  the  derived  functional  relationship  is  an  adequate  model 
of  the  actual  system. 


Elevation  versus  Damage 


Since  there  were  no  available  data  relating  the  flood  elevation  to  the  flood  damage,  the  three 
functional  relationships  obtained  earlier.  Equations  (1.7)-{1.9),  were  used  to  derive  the  required  functional 
relationship.  The  resulting  equation  is; 


L  = 


-56.65856 


0.00211 

1  [('  £+0.15423-812.6) 
1000  [\  1.28767  J 

1/0.41567 

+  10.7 

1  -  0.00582 

0.59494/0.62411 


+  9.02988  (1  10) 

where  L  is  the  flood  damage  in  $  million,  and  E  is  the  flood  elevation  in  feet.  The  resulting  plot  is  shown 
in  Figure  1-8. 
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Development  of  the  Functional  Relationships  for  Plan  4 

For  plan  4,  an  alternative  elevation-vs. -discharge  curve  was  provided  for  the  elevation  in  the 
channel  [Plate  A17,  Appendix  B,  U.S.  Army  Corps  of  Engineers  1990].  However,  the  damage  depends 
upon  the  relationship  of  discharge  vs.  elevation  in  the  floodplain.  Therefore,  we  assume  data  for  this 
relationship.  The  curve  obtained  this  way  is  then  used  to  derive  the  other  three  curves.  The  main 
relationship  that  we  are  interested  in  is  the  frequency  vs.  the  damage  since  this  will  enable  us  to  compute 
the  mean  damage  and  the  conditional  mean  damages.  Therefore,  for  plan  4,  only  the  development  of  the 
elevation-vs-  discharge  and  the  frequency-vs. -damage  curves  need  to  be  shown. 

Discharge  versus  Elevation 

The  data  for  the  discharge-vs. -floodplain  elevation  were  assumed  from  the  data  provided  for  the 
discharge-vs.-chaimel  elevation.  The  data  thus  obtained  are  shown  in  Table  1-4  and  the  corresponding  plot 
is  shown  in  Figure  1-9.  Figure  1-10  shows  the  comparative  plots  of  floodplain  elevation  vs.  discharge  for 
plan  1  and  plan  4. 


Table  1-4.  Discharge  vs.  Elevation  for  Plan  4 


Discharge 

(1000  cfs) 

Floodplain 

Elevation 

(feet) 

10.7 

812.6 

15.0 

814.1 

23.0 

816.3 

31.7 

817.0 

44.3 

817.5 

63.3 

817.8 

96.0 

818.5 
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Figure  1-9.  Discharge  vs.  Floodplain  Elevation  for  Plan  4 
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Based  on  this  data,  a  regression  analysis  was  performed  and  the  following  functional  relationship  between 
discharge  and  elevation  was  obtained  for  plan  4: 

■  0.31598 

10.7  -  0.11404  +  812.6  (1.11) 


E  =  1.55092 


D 


1000 


where  E  is  the  flood  elevation  in  feet,  and  D  is  the  discharge  in  cfs.  Figure  1-11  shows  the  comparative 
plots  of  the  actual  vs  fitted  curves,  indicating  that  the  derived  functional  relationship  is  an  adequate  model 
of  the  actual  system. 


Figure  1-11.  Acmal  vs.  Fitted  Curves  for  Discharge  vs.  Elevation— Plan  4 


Frequency  versus  Damage 

Since  the  main  curve  of  interest  is  the  fi:equency-vs. -damage  curve,  this  curve  was  developed  using 
Equations  (1.9),  (1.10),  and  (1.11).  (Note  that  the  assumption  here  is  that  the  elevation-vs.-damage  and  the 
frequency-vs. -discharge  relationships  do  not  change  due  to  the  construction  of  the  levee,  floodwalls,  and 
other  structural  measures  that  constitute  plan  4.)  If  the  design  that  provides  for  a  50-year  level  of  protection 
is  selected  for  plan  4,  then  the  frequency-vs. -damage  curve  is  truncated  at  that  level.  The  resulting 
frequency-vs.-damage  curve  for  plan  4  is  shown  in  Figure  1-12.  Figure  1-13  shows  the  comparative  plots 
of  frequency  vs.  damage  for  plan  1  and  plan  4. 
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Figure  1-13.  Comparative  Frequency  vs.  Damage  Curves  for  Plan  1  and  Plan  4 


Risk-Based  Evaluation  of 
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Since  the  expression  obtained  for  the  frequency- vs. -damage  curve  is  very  complex,  the  results 
obtained  from  that  expression  are  used  as  input  for  a  regression  analysis  in  order  to  obtain  a  simplified 
expression.  A  two-part  equation  was  used: 

(1)  a  straight  line  was  used  for  the  the  return  period  between  50  and  75  (equivalent  to  frequencies 
between  0.02  and  0.0133)  as  given  by  Eq.  (1.12a), 

L  =  521.70  (0.02 -F)  (l-12a) 

where  L  is  the  flood  damage  in  $  million,  and  F  is  the  flood  frequency;  and 

(2)  regression  analysis  was  used  for  a  return  period  greater  than  75  years,  as  given  in  Eq.  (1 . 12b), 

L  =  -50.45478  (F)°-^ -I-  8.49407  (1.12b) 

Figure  1-14  shows  the  comparative  plots  of  the  acmal  vs.  fitted  curves,  indicating  that  the  derived  functional 
relationship  is  an  adequate  model  of  the  actual  system. 


Figure  1-14.  Actual  vs.  Fitted  Curves  for  Frequency  vs.  Damage— Plan  4 
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Tncorporation  of  a  Flood  Warning  System 

Assume  that  one  design  option  for  a  flood  warning  system  is  available  to  be  added  to  the  flood 
damage  reduction  measures.  This  will  increase  two  more  integrated  flood  damage  reduction  measures,  plan 
1  +  flood  warning  system  and  plan  4  +  flood  warning  system.  The  alternative  of  adopting  plan  1  +  flood 
warning  system  is  essentially  an  alternative  of  adopting  only  a  nonstructural  measure,  since  plan  1  is  the  zero- 
cost  plan  or  the  without-project-condition  alternative.  The  alternative  of  plan  4  +  warning  system  is  an  actual 
Integrated  measure. 

Introduction  of  a  flood  warning  system  will  change  the  relationship  between  elevation  and  damage. 
Specifically,  the  flood  loss  reduction  is  given  by  Equation  (1.4)  as  a  function  of  the  flood  elevation. 

The  value  of  the  maximum  flood  loss  for  a  community,  MD,  and  the  functional  form  of  unit  damage 
function  can  be  obtained  from  Equation  (1 .4).  In  this  specific  example,  we  notice  that  the  base  elevation  of 
the  floodplain  zone  is  812.5  feet,  the  maximum  flood  loss  for  the  community,  MD,  is  equal  to  $9.02988 
millions,  and  the  functional  form  of  the  unit  damage  function  is 


b(h-y)  =  1  -  6.27456 


[  0.00211 

1  fj'  /j-812.44577] 
1000  [(  1.28767  J 

2.40576  1 

+  10.7j 

1  -  0.00582 

0.95326 


(1.13) 


The  maximum  evacuation  cost  for  the  community,  MC,  is  assumed  to  be  equal  to  $50,000  and  the 
unit  reduction  function  is  assumed  to  be 

MR(h  -  812.44577)  =  0.25  +  0.04(h  -  812.44577)  -  0.00333(h  -  812.44577)'  (1.14) 

Since  the  flood  loss  reduction  LRnin  Equation  (1.4)  is  a  linear  function  of  the  response  fraction  0,  we 
only  study  the  case  of  maximum  flood  loss  reduction  when  a  full  response  is  present,  i.e.,  0  =  1.  The  cases 
where  a  frill  response  is  not  present  can  be  easily  found  by  interpolating  the  curves  of  elevation  vs.  damage 
with  and  without  a  flood  warning  system. 

In  summary,  the  loss  reduction  when  a  flood  warning  system  is  introduced  can  be  expressed  as  a 
function  of  flood  elevation 


9.02988  -  56.65856 

i 

0.00211 

0.95326 

1  I 

^  /i-812.44577'| 

'  '  '  '  1  ^  Iv*/ 

^  -  0.00582 

1000  I 

1.28767  j 

*  {0.25  +  0.04(h  -  812.44577)  -  0.00333(h  -  812.44577'}  -  0.05  (1.15) 

Subtracting  LRufrom  the  damage  coordinate  in  Figure  1-8  (Equation  1.10)  for  each  elevation  yields 
the  curve  of  damage  vs.  elevation  when  a  flood  warning  system  is  introduced.  Figure  1-15  shows  this  curve. 
Figure  1-16  shows  the  comparative  plots  of  damage  vs.  elevation  with  and  without  a  warning  system. 


24 


Bisk-Based  Evaluation  of 
Flood  Warning  and  Preparedness  Systems 


Figure  1-16.  Comparative  Plots  of  Elevation  vs.  Damage  for  Plan  1 
With  and  Without  Flood  Warning  System 
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Combining  Figure  1-3  (Equation  3.1),  Figure  1-5  (Equation  1.8)  and  Figure  1-15  (Equation  1.15 
subtracted  from  Equation  1.10)  gives  us  the  curve  of  damage  vs.  frequency  in  Figure  1-17  for  the  option 
of  plan  1  +  flood  warning  system. 

Since  the  expression  obtained  for  the  frequency-vs. -damage  curve  is  very  complex,  the  results 
obtained  from  that  expression  are  used  as  input  for  a  regression  analysis  in  order  to  obtain  a  simplified 
expression.  This  regression  analysis  gives  us  the  following  functional  relationship: 

L  =  -32.13928  (F)'”^'’* -1-  6.08851  (1.16) 

Figure  1-18  shows  the  comparative  plots  of  the  actual  vs  fitted  curves,  indicating  that  the  derived  functional 
relationship  is  an  adequate  model  of  the  actual  system. 

Similarly,  combining  Fig.  1-9  (Eq.  1.11),  Fig.  1-5  (Eq.  1.8)  and  Fig.  1-15  (subtracting  Eq.  1.15 
from  Eq.  1.10)  gives  us  the  curve  of  damage  vs.  frequency  in  Fig.  1-19  for  the  option  of  plan  4  -t-  flood 
warning  system. 

Since  the  expression  obtained  for  the  frequency-vs. -damage  curve  is  very  complex,  the  results 
obtained  from  that  expression  are  used  as  input  for  a  regression  analysis  in  order  to  obtain  a  simplified 
expression.  A  two-part  equation  was  used: 

(1)  a  straight  line  was  used  for  a  return  period  between  50  and  75  (equivalent  to  frequencies 
between  0.02  and  0.0133)  as  given  by  Eq.  (1.17a), 

L  =  340.58  (0.02  -  F)  (1.17a) 

where  L  is  the  flood  damage  in  $  million,  and  F  is  the  flood  frequency;  and  (2)  regression  analysis 
was  used  for  a  return  period  greater  than  75  years,  as  given  in  Eq.  (1.17b), 

L  =  -29.39980  (F)'’^^^ -f  5.46498  (1.17b) 

Figure  1-20  shows  the  comparative  plots  of  the  actual  vs  fitted  curves,  indicating  that  the  derived  functional 
relationship  is  an  adequate  model  of  the  actual  system. 

Computation  of  Conditional  and  Unconditional  Expected  Damages 

Plan  1 


Recall  from  Equations  (1.5)  and  (1.6)  that  the  conditional  expected  damage  f4(L|a)  is  an  average 
of  damage  over  the  risk  of  extreme  events  given  that  the  ath  percentile  of  damage  is  exceeded,  and  that  the 
unconditional  expected  damage  f5(L)  is  the  overall  average  damage.  In  order  to  compute  the  conditional  and 
unconditional  expected  damages,  we  must  obtain  the  probability  density  function  (pdf)  of  damages.  We  can 
obtain  this  from  the  frequency-vs. -damage  curve.  This  relationship  is  given  by  Eq.  (1.9),  which  is 
simplified  and  shown  in  Eq.  (1.18): 
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Figure  1-19.  Frequency  vs.  Damage  for  Plan  4  with  Flood  Warning  System 


Figure  1-20.  Acmal  vs.  Fitted  Curves  for  Frequency  vs.  Damage— Plan  4  with  Flood  Warning  System 
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where  L  is  the  flood  damage  in  $  million,  F  is  the  flood  frequency,  a  equals  -56.65856,  b  equals  0.59494, 
and  c  equals  9.02988. 


Let  /(L)  be  the  pdf  of  L;  then  from  Equation  (1 . 1 8)  we  can  obtain 


./(L)  =  F  = 


(1.19) 


The  probability  of  flood,  p(L),  is  given  by  the  area  under  this  curve.  This  area  is  computed  by  integrating 
it  from  0  to  c  (note  that  c  is  an  upper  bound  for  this  equation  since  M  is  equal  to  c  when  F  approaches  0); 


The  mean  value  of  the  damage,  {fsiL)},  can  be  computed  as  follows: 


/5(L)  =  J  Lf(L)  dL 

0 

=  $0.3772  million 


(1.21) 


The  conditional  expected  value,  {f  4(Ll  a},  is  given  by 


I  LfiL)  dL 
ML\a)  =  - 

jmdL 

K 


where  a  is  the  partition  point  on  the  probability  axis,  and 

^■1(1  *b)lb\  bl{Ub) 


l-'O  ^ 


[a+p(L)-l][l+b] 
a  b 


+  c 


(1.22) 


(1.23) 
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is  the  corresponding  partitioning  point  on  the  damage  axis.  For  a  =  0.99,  we  have  L.  =  5.7716  and 
/4(L|a=0.99)  =  $6.6568  million 

For  a  =  0.9,  we  have  L„  =  1.3385,  and 

/4(L|a=0.9)  =  $3.4281  million 


Plan  4 


We  adopt  the  same  approach  used  in  the  previous  section  for  plan  1  in  order  to  compute  the 
required  measures  for  plan  4.  The  first  step  is  to  obtain  the  probability  density  function  (pdf)  of  damages. 
We  can  obtain  this  from  the  frequency-vs. -damage  curve.  This  relationship  is  a  two-part  equation  and  is 
given  by  Equations  (1.12a)  and  (1.12b).  Equation  (1.12b)  is  simplified  and  shown  in  Equation  1.24: 

L  =  a(F)'’-l-c  (1.24) 

where  L  is  the  flood  damage  in  $  million,  F  is  the  flood  frequency,  a  equals  -50.45478,  b  equals  0.53290, 
and  c  equals  8.49407. 

Using  Equations  1.12a  and  1.12b  we  obtain  the  pdf  of  damage, /(L),  as 


/(L) 


0.02  - 


521.70 


if  L  <  3.4780 


if  L  >  3.4780 


(1.25) 


The  probability  of  flood,  p(L),  is  given  by  the  area  under  this  curve.  This  area  is  computed  by  integrating 
it  from  0  to  c  (note  that  c  is  an  upper  bound  for  this  equation  since  M  is  equal  to  c  when  F  approaches  0): 


P(L) 


3.4780 

/  {“•“ 

0  '• 


0.0580  +  0.0229 


0.08089 


(1.26) 
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The  mean  value  of  the  damage  {f  5(L)}  can  be  computed  using  Equation  (1.21)  as 
f5(L)  =  $0.2035  million 

The  conditional  expected  value  {/4(Lla)}  is  given  by  Equation  (1.22);  for  a  =  0.99,  we  have 


|[a+/?(L)-l-0.0580][l+&] 
I  a  b 


d-c 


(1  +&)/hl  &/(!+&) 


+  c 


4.7345 


(1.27) 


and 

f4(L|a=0.99)  =  $5.7046  million 
and  for  a  =  0.9,  we  have  L.  =  0,  and 


f4(L|a=0.9)  = 


3.4780 

I  |0.02- 


521.70 


dL  + 


3.4780  ^ 

3.4780  ^  c 

f  jo.o2--^u. 


1/6 


dL 


Vb 


dL 


$2.5153  million 


(1.28) 


Plan  1  +  Flood  Warning  System 

We  adopt  the  same  approach  used  in  the  previous  section  for  plan  1  in  order  to  compute  the 
required  measures  for  plan  1  +  flood  warning  system.  The  first  step  is  to  obtain  the  probability  density 
function  (pdf)  of  damages.  We  can  obtain  this  from  the  frequency-vs. -damage  curve.  This  relationship  is 
given  by  Equation  (1.16),  which  is  simplified  and  shown  in  Equation  (1.29): 


L  =  a(F)‘  +  c  0-29) 

where  L  is  the  flood  damage  in  $  million,  F  is  the  flood  frequency,  a  equals  -32.13928,  b  equals  0.54301, 
and  c  =  6.08851.  Using  Equation  (1.20)  we  obtain 

p(L)  =  0.10008 
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The  mean  value  of  the  damage  {(/‘jCL)}  can  be  computed  using  Equation  (1.21)  as 
/5(L)  =  $0.1586  million 

The  conditional  expected  value  {(ffL  |  a)}  is  given  by  Equation  (1 .22)  and  the  partition  point  {L  J  is  given 
by  Equation  (1.23).  For  a  =  0.99,  we  have  L„  =  3.3816  and 

/4(L|a=0.99)  =  $4.0863  million 

For  a  =  0.9,  we  have  L„  =  0.0018,  and 

/4(L|a=0.9)  =  $1.5862  million 


Plan  4  +  Flood  Warning  System 

We  adopt  the  same  approach  used  in  the  previous  section  for  plan  4  in  order  to  compute  the 
required  measures  for  plan  4  +  flood  warning  system.  The  first  step  is  to  obtain  the  probability  density 
function  (pdf)  of  damages.  We  can  obtain  this  from  the  frequency-vs. -damage  curve.  This  relationship  is 
a  two-part  equation  and  is  given  by  Eqs.  (1.17a)  and  (1.17b).  Equation  (1.17b)  is  simplified  and  shown  in 
Eq.  (1.30). 

L  -  a(F)'’  +  c  (1.30) 

where  L  is  the  flood  damage  in  $  million,  F  is  the  flood  frequency,  a  equals  -29.39980,  b  equals  0.51246, 
and  c  equals  5.46498. 


Using  Equations  (1.17a)  and  (1.17b)  we  obtain  the  pdf  of  damage, /(L),  as 


/(L)  = 


r 


< 


V 


jo. 02 


340.58 


if  L  <  2.2705 


if  L  >  2.2705 


(1.31) 


The  probability  of  flood,  p(L),  is  given  by  the  area  under  this  curve.  This  area  is  computed  by  integrating 
it  fi-om  0  to  c  (note  that  c  is  an  upper  bound  for  this  equation  since  M  is  equal  to  c  when  F  approaches  0). 
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P(L) 


0.03784  +  0.01424 


0.05208 


(1.32) 


The  mean  value  of  the  damage  {fs(L)}  can  be  computed  using  Equation  (1 .21)  as 
fsCL)  =  $0.0839  million 

The  conditional  expected  value  {ffh\a)}  is  given  by  Equation  (1.22).  For  a  =  0.99,  we  have 


La  = 


|[«+p(L)-l-0.03784][l+&] 
a  b 


d-c 


(1  +&)/&]  bl{Ub) 


+  c 


2.6307 


(1.33) 


and 


/4(L|a=0.99)  =  $3.3478  million 
and  for  a  =  0.9,  we  have  =  0,  and 


f4(L|a=0.9) 


2.2705 

I  “ 


2.2705 


$1.6116  million 


(1.34) 


Tradeoff  Analysis 

Once  the  conditional  and  unconditional  expected  values  for  the  different  plans  are  computed,  we 
can  perform  a  tradeoff  analysis  in  terms  of  costs  and  damages.  The  compiled  results  are  shown  in  Table 
1-5. 
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Table  1-5.  Summary  of  Results:  Tradeoffs  Among  Cost,  Expected  Damage  (fj),  and 
Risk  of  Extreme  Events  (f4)  for  the  Four  Alternatives 


Average 
Annual  Cost 
($  million) 

/5(L) 

/4(L|a=0.9) 

/4(L|(X=0.99) 

Plan  1 

0.000 

0.377 

3.428 

6.657 

Plan  l-fW 

0.050 

0.159 

1.586 

4.086 

Plan  4 

0.865 

0.204 

2.515 

5.705 

Plan  4+W 

0.915 

0.084 

1.412 

3.348 

Since  plan  1  is  the  option  of  doing  nothing,  it  does  not  have  an  associated  cost.  The  average  annual  cost 
for  plan  4  is  given  as  $0,865  million  for  a  50-year  level  of  protection  [Table  1-18,  page  1-29,  U.S.  Army 
Corps  of  Engineers  1990].  The  average  annual  cost  of  the  flood  warning  system  is  assumed  to  be  $50,000. 
Figure  1-21  shows  the  resulting  tradeoffs  when  using  the  expected  value  alone.  The  solid  line  shows  the 
Pareto  optimal  frontier.  Figures  1-22  and  1-23  show  the  tradeoffs  for  a  =  0.9  and  0.99,  respectively. 
Figure  1-24  shows  these  tradeoffs  together.  Note  that  the  option  of  plan  4  without  the  warning  system  is 
noninferior,  or  a  viable  option,  in  considering  only  structural  measures;  the  same  option  becomes  inferior 
when  considering  the  warning  system  options.  Plan  1,  plan  1 +W,  and  plan  4-l-W  constitute  the  noninferior 
set  of  options  in  the  combined  structural/nonstructural  analysis.  The  combined  analysis  of  structural  and 
nonstructural  measures,  incorporating  the  risk  of  extreme  events,  clearly  demonstrates  the  relative 
inefficiency  of  plan  4  without  the  warning  system,  an  important  result  that  would  not  have  come  from  a 
traditional  approach. 
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Figure  1-21.  Cost  vs.  Damage  Tradeoff  for  the  Expected  Value  (fj) 
(Note  that  there  are  three  Pareto  optimal  (efficient)  alternatives.) 
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Figure  1-23.  Cost  vs.  Conditional  Expected  Damage  (f,)  Tradeoff  for  a  =  0.99 
(Note  the  three  Pareto  optimal  (efficient)  alternatives.) 


/5(L) 

/4(LI  a=0.9) 


Damages  ($  million) 


Plan  4+W 
Plan  4 


Plan  1+W 
Plan  1 
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Figure  1-24.  Tradeoffs  Among  Cost,  Expected  Damage  (fj),  and  Conditional  Expected 

Damage  (Q  for  the  Four  Alternatives 

(Risk  of  Extreme  Events  (f,)  evaluated  at  two  partitioning  levels  (ct  =  0.9,  a  —  0.99).) 
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Conclusions 

Integrated  flood  damage  reduction  alternatives  using  both  structural  measures  and  flood  warning 
systems  will  reduce  the  vulnerability  of  a  community  to  flood  damage  and  add  resiliency  because  of  added 
redundancy,  thus  providing  more  options  in  decreasing  expected  flood  loss  with  tradeoff  consideration  of  the 
associated  costs.  The  approach  in  this  report  builds  on  the  previous  Corps  of  Engineers  work  and  is  very  easy 
to  adopt  and  implement.  Given  the  analysis  for  each  structural  and  flood  warning  system,  the  combined 
analysis  is  simpler.  The  additional  data  requirement  is  minimal.  The  incorporation  of  the  measure  of  the  risk 
of  extreme  events  along  with  the  expected  flood  loss  in  a  multiple-objective  framework  offers  deeper  insight 
in  determining  the  best  flood  damage  reduction  strategy. 
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Parti 

Multiobjective  Decision  Tree  Analysis: 

Technical 

Introduction 

Decision-tree  analysis  has  emerged  over  the  years  as  an  effective  and  useful  tool  in  decisionmaking. 
More  than  two  decades  ago,  Howard  Raiffa  [1968]  published  the  first  comprehensive  and  authoritative  book 
on  decision-tree  analysis.  Ever  since,  its  applications  to  a  variety  of  problems  fi’om  numerous  disciplines 
have  grown  by  leaps  and  bounds  [see  Sage  1977  and  Hamburg  1970].  Advances  in  science  and  in  scientific 
approaches  to  problem  solving  are  often  made  on  the  basis  of  earlier  works  of  others.  In  this  case,  the 
foundation  for  Raiffa*  s  contributions  to  decision  tree  analysis  can  be  traced  to  the  works  of  Bernoulli  on 
utility  theory  [see  von  Neumann  and  Morgenstem  1944;  Edwards  1967;  Savage  1954;  Adams  1960;  Arrow 
1963;  Shubik  1964;  Luce  and  Suppes  1965;  and  others].  This  chapter,  in  an  attempt  to  build  on  the  above 
seminal  works,  extends  and  broadens  the  concept  of  decision-tree  analysis  to  incorporate:  (a)  multiple, 
noncommensurate  and  conflicting  objectives,  (b)  impact  analysis,  and  (c)  the  risk  of  extreme  and 
catastrophic  events.  Indeed,  the  current  practice  often  involves  one-sided  use  of  decision  trees  —  optimizing 
a  single-objective  function  and  commensurating  infrequent  catastrophic  events  with  more  frequent 
noncatastrophic  events  using  the  common  unconditional  mathematical  expectation. 


Multiple  Objectives 

The  single-objective  models  that  had  been  advanced  in  the  fifties,  sixties,  and  seventies  are  today 
considered  by  many  to  be  unrealistic,  too  restrictive,  and  often  inadequate  for  most  real-world  complex 
problems.  The  proliferation  of  books,  articles,  and  conferences  and  courses  during  the  last  decade  or  two 
on  what  has  come  to  be  known  as  multiple-criteria  decisionmaking  (MCDM)  is  a  vivid  indication  of  this 
somber  realization  and  of  the  maturation  of  the  field  of  decisionmaking  [see  Chankong  and  Haimes  1983]. 
In  particular,  an  optimum  derived  from  a  single-objective  mathematical  model,  including  that  which  is 
derived  from  a  decision  tree,  often  may  be  far  from  representing  reality  -  thereby  misleading  the  ^alyst(s) 
as  well  as  the  decisionmaker(s).  Fundamentally,  most  complex  problems  involve,  among  other  things,  the 
minimization  of  costs,  the  maximization  of  benefits  (not  necessarily  in  monetary  values),  and  the 
minimization  of  risks  of  various  kinds.  Decision  trees,  which  are  a  powerful  mechanism  for  the  analysis 
of  complex  problems,  can  better  serve  both  the  analysts  and  the  decisionmakers  when  they  are  extended 
to  deal  with  the  above  multiple  objectives. 


Impact  Analysis 

On  a  long-term  basis,  managers  and  other  decisionmakers  are  often  rewarded  not  because  they  have 
made  many  optimal  decisions  in  their  tenure;  rather,  they  are  honored  and  thanked  for  avoiding  adverse 
and  catastrophic  consequences.  If  one  accepts  this  premise,  then  the  role  of  impact  analysis  --  studying  and 
investigating  the  consequences  of  present  decisions  on  future  policy  options  —  might  be  as  important,  if 
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not  actually  more  so,  than  generating  an  optimum  for  a  single-objective  model  or  identifying  a  Pareto 
optimum  set  for  a  multiobjective  model.  Certainly,  when  the  ability  to  generate  both  is  present,  having  an 
appropriate  Pareto  optimum  set  and  knowing  the  impact  of  each  Pareto  optimum  on  future  policy  options 
should  enhance  the  overall  decisionmaking  process  within  the  decision-tree  framework. 


The  Risk  of  Extreme  and  Catastrophic  Events 


Risk,  which  is  a  measure  of  the  probability  and  severity  of  adverse  effects,  has  until  recently  been 
commonly  quantified  via  the  expected-value  formula.  This  formula  essentially  precommensurates  events 
of  low  frequency  and  high  damage  with  events  of  high  frequency  and  low  damage.  Although  learned 
students  of  risk  analysis  recognize  the  disparity  between  the  above  fallacious  representation  of  extreme  and 
catastrophic  events  and  the  perception  of  these  events  by  individuals  or  the  public  at  large,  many  continue 
to  use  this  approach.  The  trend,  however,  is  moving  toward  the  conditional-expected-value  approach,  where 
extreme  and  catastrophic  events  are  partitioned,  isolated,  quantified  in  terms  of  the  conditional  expectation 
(e.g.,  using  concepts  from  the  statistics  of  extremes),  and  then  evaluated  along  with  the  common  expected 
value  of  risk  or  damage  [Asbeck  and  Haimes  1984;  Haimes  1985;  Karlsson  and  Haimes  1988a,  1988b]. 

The  partitioned  multiobjective  risk  method  (PMRM)  developed  by  Asbeck  and  Haimes  [1984] 
separates  extreme  events  from  other  noncatastrophic  events,  and  thus  provides  the  decisionmaker(s)  with 
additional  valuable  and  useful  information.  In  addition  to  using  the  traditional  expected  value,  the  PMRM 
generates  a  number  of  conditional  expected- value  functions,  termed  here  risk  functions,  which  represent 
the  risk,  given  that  the  damage  falls  within  specific  probability  ranges  (or  damage  ranges).  Assume  Aat  the 
risk  can  be  represented  by  a  continuous  random  variable  X  with  a  known  probability  density  function 
p^(x;Sj),  where  ^  (j  =  l,...,q)  is  a  control  policy.  The  PMRM  partitions  the  probability  axis  into  three 
ranges.  Denote  the  partitioned  points  on  the  probability  axis  by  Cj  (i  =  1,2).  For  each  a,  and  each  policy 
Sj,  it  is  assumed  that  there  exists  a  unique  damage  by  such  that 


Px(Pij;Sj)  =  a.  (2.1) 

where  P,  is  the  cumulative  distribution  function  of  X.  These  Py  (with  Poj  and  Pjj  representing,  respectively, 
the  lower  bound  and  upper  bound  of  the  damage)  define  the  conditional  expectation  as  follows: 


fi(Sj)  =  E{x  I  p,(x;Sj),  X  e  [P,.2j,Pi.,j]} 


or 


fi(Sj)  = 


Pi-.j 

/  xp^(x;s.)dx 

P.-..i 

/  p^(x;Sj)dx 
Pi-2J 


(i  =  2,3,4;  j  =  l,...,q) 


(i  =  2,3,4;  j  =  l,...,q) 


(2.2) 


(2.3) 
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where  fj,  fj,  and  f4  represent  the  risk  with  high  probability  of  exceedance  and  low  damage,  the  risk  with 
medium  probability  of  exceedance  and  medium  damage,  and  the  risk  with  low  probability  of  exceedance 
and  high  damage,  respectively.  The  unconditional  (conventional)  expected  value  of  X  is  denoted  by  fjCsj). 
The  relationship  between  the  conditional  expected  values  (fj,  fj,  fj  and  the  unconditional  expected  value 
(fj)  is  given  by 

fsCsj)  =  02f2(Sj)  +  ejfjCSj)  +  04f4Sj)  (2.4) 

where  0i(i  =  2,3,4)  is  the  denominator  of  Eq.  (3).  From  the  definition  of  Py,  it  can  been  seen  that  0i  ^  0 
is  a  constant,  and  02  +  03  +  04  =  1- 

Combining  either  the  generated  conditional  expected  risk  function  or  the  unconditional  expected 
risk  function  with  the  cost  objective  function,  f,,  creates  a  set  of  multiobjective  optimization  problems: 

min  [f„fj‘  (i  =  2,3,4,5)  (2-5) 

where  the  superscript  t  denotes  the  transpose  operator.  This  formulation  offers  more  information  about  the 
probabilistic  behavior  of  the  problem  than  the  single  multiobjective  formulation  min  [fi.fj]'.  The  tradeoffs 
between  the  cost  function  f,  and  any  risk  function  f;  (i  e  {2, 3,4,5})  allow  decisionmakers  to  consider  the 
marginal  cost  of  a  Small  reduction  in  the  risk  objective,  given  a  particular  level  of  risk  assurance  for  each 
of  the  partitioned  risk  regions  and  given  the  unconditional  risk  function,  fj.  The  relationship  of  the  tradeoffs 
between  the  cost  function  and  the  various  risk  functions  is  given  by 

1Ai5  =  02/^12  +  03/^13  +  04^^14  (2-6) 

where 

=  -afM(i  =  2, 3,4,5)  (2-7) 

and  02,  03,  and  Q  are  as  defined  earlier.  A  knowledge  of  this  relationship  among  the  marginal  costs 
provides  the  decisionmakers  with  insights  that  are  useful  for  determining  an  acceptable  level  of  risk. 


Methodological  Approach 

Extension  to  Multiple  Objectives 

Similar  to  the  decision-tree  in  conventional  single-objective  analysis,  a  multiobjective  decision  tree 
(Fig.  2-1)  is  composed  of  decision  nodes  and  chance  nodes.  Each  pair  of  an  alternative  and  a  state  of  nature, 
however,  is  now  characterized  by  a  vector-valued  performance  measure. 

At  a  decision  node,  usually  designated  by  a  square,  the  decisionmaker  selects  one  course  of  action 
from  the  feasible  set  of  alternatives.  We  assume  that  there  are  only  a  finite  number  of  alternatives  at  each 
decision  node.  These  alternatives  are  shown  as  branches  emerging  to  the  right  side  of  the  decision  node. 
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The  performance  vector  associated  with  each  alternative  is  written  along  the  corresponding  branch.  Each 
alternative  branch  may  lead  to  another  decision  node,  a  chance  node,  or  a  terminal  point. 

A  chance  node,  designated  by  a  circle  on  the  tree,  indicates  that  a  chance  event  is  expected  at  this 
point;  that  is,  one  of  the  states  of  nature  may  occur.  We  consider  two  cases  in  this  study:  a)  a  discrete  case, 
where  the  number  of  states  of  nature  is  assumed  finite;  and  b)  a  continuous  case,  where  the  possible  states 
of  nature  are  assumed  continuous.  The  states  of  nature  are  shown  on  the  tree  as  branches  to  the  right  of  the 
chance  nodes,  and  their  known  probabilities  are  written  above  the  branches.  The  states  of  nature  may  be 
followed  by  another  chance  node,  a  decision  node,  or  a  terminal  point. 

Allowing  for  the  evaluation  of  the  multiple  objectives  at  each  decision  node  constitutes  an  important 
feature  of  the  approach  presented  here.  This  is  a  significant  extension  of  the  average-out-and-folding-back 
strategy  used  in  conventional  single-objective  decision-tree  methods. 

To  allow  for  this  extension,  we  first  define  a  k-dimensional  vector- valued  performance  measure 
associated  with  an  action  a„  and  a  state  of  nature  0„  as  follows, 

r(a„,0„)  =  [r,(a„,0„),r2(a„,0„),...,rk(a„,0jr  (2-8) 

A  point  r  =  [ri,r2,...,rj‘  in  the  objective  function  space  is  said  to  be  noninferior  if  there  does  not 
exist  another  feasible  point  r'  =  [r'i,r'2,...,r'J‘ such  that 

r\  <  r,  (i  =  1,2,.. .,k)  (2.9) 

with  at  least  one  strict  inequality  holding  for  i  =  l,2,...,k. 

The  sequential  structure  of  multiobjective  decision  trees  necessitates  the  introduction  of  a  vector  of 
operators  that  combines  the  vectors  of  performance  measures  of  successive  decision  nodes.  Let  o  denote 
a  k-dimensional  vector  of  binary  operators  which  are  to  be  applied  to  elements  corresponding  to  the  same 
components  of  any  two  vectors  of  a  performance  measure.  For  example,  if 

r,  =  [2,3]',  r^  =  [-3,2]',  o  =  (+,•) 


then 


Ti  0  r^  =  [2-3,  3*2]t  =  [-1,6], 

The  solution  procedure  for  multiobjective  decision  trees  is  stated  in  three  steps: 

Step  1.  Chart  the  decision  tree  for  the  problem  under  study. 

Step  2.  Assign  an  a  priori  probability  or  calculate  the  posterior  probability  for  each  chance  branch. 
Assign  the  vector-valued  performance  measure  for  each  pair  of  an  alternative  and  a  state  of  nature.  (Or  map 
the  vector-valued  performance  measure  to  each  of  the  terminal  points  of  the  tree.) 
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Step  3.  Start  from  each  terminal  point  of  the  tree  and  fold  backward  on  the  tree. 

At  each  decision  node,  n,  and  at  each  branch  emerging  to  the  right  side  of  the  decision  node,  find 
the  corresponding  set  of  vector-valued  performance  measures,  r(ai"),  for  each  alternative,  aj,  and  identify 
the  set  of  noninferior  solutions  by  solving 


r"  =  min  U  r(ai")  (2. 10) 

i 

Remark:  In  multiobjective  decision-tree  analysis,  instead  of  having  a  single  optimal  value  associated  with 
a  single-objective  decision  tree,  we  have  r",  a  set  of  vector-valued  objective  values  of  noninferior  decision 
alternatives  at  decision  node  n. 

At  each  chance  node  m  and  at  branches  emerging  to  the  right  side  of  the  chance  node,  find  the 
corresponding  set  of  vector-valued  performance  measures,  r™,  for  each  state  of  nature  0  j“,  and  then 
calculate  the  vector-valued  expected-performance  measure,  or  other  specified  vector-valued  "risk" 
performance  measure,  which  is  denoted  by  r“ 

r“  =  min  E*{rj“}  (2.11) 

j 

Remarks: 

a)  In  single-objective  decision-tree  analysis,  there  is  no  choice  process  at  the  chance  nodes,  since 
only  an  averaging-out  process  takes  place  there.  In  multiobjective  decision-tree  analysis,  a  set  of  Pareto 
optimum  alternatives,  rj™,  is  associated  with  each  branch  emerging  from  chance  node  m.  If  each  set  of 
Pareto  optimal  solutions  rj”  has  elements,  then  there  exist  Pjldj™}  combinations  of  decision  rules  needing 
to  be  averaged-out,  and  a  vector  minimization  must  be  performed  to  discard  from  further  consideration  the 
resulting  inferior  combinations. 

b)  The  superscript  s  in  E’  denotes  the  sth  averaging-out  strategy;  in  particular,  E’  (for  s  =  5) 
denotes  the  conventional  expected-value  operator,  and  E'*  (for  s  =  4)  denotes  the  operator  of  conditional 
expected  value  in  the  region  of  extreme  events  (which  will  be  discussed  in  detail  in  a  later  section). 

Step  3  is  repeated  until  the  set  of  noninferior  solutions  at  the  starting  point  of  the  tree  is  obtained. 
Impact  of  Experimentation 

The  impact  of  an  added  piece  of  information  (obtained,  e.g.,  through  experimentation)  on  different 
objectives  is  now  addressed,  and  the  value  of  the  information  is  quantified  by  a  vector-valued  measure.  In 
conventional  decision-tree  analysis,  whether  or  not  an  experiment  should  be  performed  depends  on  an 
assessment  of  the  expected  value  of  experimentation  (EVE),  which  is  the  difference  between  the  expected 
loss  without  experimentation  and  the  expected  loss  with  experimentation.  If  the  EVE  is  negative, 
experimentation  is  deemed  unwarranted;  otherwise,  the  experiment  that  yields  the  lowest  loss  is  selected. 
In  multiobjective  decision-tree  analysis,  the  monetary  index  does  not  constitute  the  sole  consideration; 
rather,  the  value  of  experimentation  is  judged  in  a  multiobjective  way  where,  in  many  cases,  the  noninferior 
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frontiers  generated  with  and  without  experimentation  do  not  dominate  each  other.  The  added 
experimentation  in  these  cases  reshapes  the  feasible  region  (and  thus  the  noninferior  frontier)  and  generates 
new  and  better  options  for  the  decisionmakers  (Fig.  2-2). 

Example  for  the  Discrete  Case 

Problem  Definition 

The  example  problem  discussed  here  concerns  a  simplified  flood  warning  and  evacuation  system. 
Three  possible  actions,  (a)  evacuation,  (b)  issuing  a  flood  watch,  and  (c)  doing  nothing,  are  under 
consideration.  There  are  cost  factors  associated  with  the  first  two  options.  The  decision  tree  covers  two  time 
periods,  and  the  cost  associated  with  each  option  is  a  function  of  the  period  in  which  the  action  is  taken. 
The  complete  decision  tree  for  the  problem  is  shown  in  Fig.  2-3.  The  following  assumptions  are  made: 

a)  There  are  three  possible  actions  with  associated  costs  for  the  first  period: 

1)  issuing  an  evacuation  order  at  a  cost  of  $5  million  [EVl], 

2)  issuing  a  flood  watch  at  a  cost  of  $1  million  [WAl],  and 

3)  doing  nothing  at  no  cost  [DNl]. 

b)  For  the  second  period  the  actions  and  the  corresponding  costs  are: 

1)  issuing  an  evacuation  order  at  a  cost  of  $3  million  [EV2], 

2)  issuing  a  flood  watch  at  a  cost  of  $0.5  million  [WA2],  and 

3)  doing  nothing  at  no  cost  [DN2]. 

c)  The  flood  stage  is  reached  at  water  flow  (W)  =  50,000  cfs. 

d)  There  are  three  underlying  probability  density  functions  (pdfs)  for  the  water  flow: 

1)  W  ~  lognormal  (10.4,1),  represented  as  LN„ 

2)  W  ~  lognormal  (9.1,1),  represented  as  LNj,  and 

3)  W  ~  lognormal  (7.8,1),  represented  as  LN,. 

The  prior  probabilities  that  any  of  these  pdfs  is  the  actual  pdf  are  equal. 

e)  There  are  four  possible  events  at  the  end  of  the  first  period: 

1)  A  flood  (W  >  50,000  cfs)  occurs. 
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Figure  2-3.  Decision  Tree  for  the  Discrete  Case 
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2)  The  water  flow  is  greater  than  that  of  the  previous  period  (15,000  cfs  s  W  ^  50,000  cfs) 
represented  as  Wl.  ’ 


3)  The  water  flow  is  in  the  same  range  as  that  of  the  previous  period  (5000  eft  <:  W  ^  15,000  cfs) 
represented  as  W2.  ’ 


4)  The  water  flow  is  lower  than  that  of  the  previous  period  (W  ^  5000  eft),  represented  as  W3. 

f)  L  =  7  and  C  =  $7,000,000  are  respectively  the  maximum  possible  loss  of  lives  and  properties,  given 
no  flood  warning.  All  other  costs  are  shown  in  Fig.  2-3. 

Calculation  of  Probabilities  for  the  First  Period 

Chance  node  Cl 


To  calculate  the  probabilities  of  a  flood  or  no  flood  event  at  the  end  of  the  second  period  (see  Fig. 
2-4),  we  use  the  facts  that  the  possible  pdf  of  the  water  flow  (W)  is  LNj  with  probability  1/3,  i  =  1,2,3, 
and  that  the  flood  stage  is  at  W  =  50,000  cfs.  The  probability  of  a  flood  event  can  be  calculated  as  follows: 

3 

prob.(flood)  =  52  prob.(flood  |  Lni)prob.(LNi) 

1=1 

3 

=  E  (1/3)  prob.(x  ^  50,000cfs  I  LNj)  (2.12a) 

where 


Equation  (12b)  is  converted  into  a  standard  normal  distribution  by  using 

exp[-  {ln(x)  -  nl}^/2afldx 


prob.(x  ^  50,000cfs  |  Lnj)  = 


50,000 


(2.12b) 


Figure  2-4.  Averaging  out  at  Chance  Node  Cl  (Discrete  Case) 
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z  =  [ln(x)  -  AiJ/o, 


(2.13) 


yielding 


prob.(x  ^  50,000cfs  |  Lnj)  — 


{ln50000-ii.i/a^) 


\/2n 


(2.14) 


Equation  (14)  is  evaluated  using  standard  normal  distribution  tables.  This  yields 
prob. (flood)  =  prob.(x  ^  50,000)  =  0.1271 
Chance  nodes  C2  and  C3 

Nodes  C2  and  C3  each  present  four  possible  events  at  the  beginning  of  the  second  period:  namely, 
a  flood  event,  the  water  flow  is  higher,  the  water  flow  is  the  same,  and  the  water  flow  is  lower 
orevious  period  (see  Fig.  2-5).  The  distribution  of  water  flow  at  the  end  of  the  first  ^ 

Lsumption  d.  The  probability  of  each  event  is  calculated  using  Eqs.  (2-12),  (2-13)  and  (2-14)  with  modi  e 

integral  intervals: 

prob.(flood)  =  prob. (50,000  :£  x  ^  ~)  =  0.1271 
prob.(higher)  =  prob.(15,000  ^  x  ^  50,000)  =  0.2466 
prob.(same)  =  prob.(5000  ^  x  s  15,000)  =  0.2686 


prob. (lower)  =  prob.(0  <  x  ^  5000)  —  0.3577 


Calculation  of  Probabilities  for  the  Second  Period 

Regardless  of  whether  a  watch  action  (WAI)  or  do  nothing  (DNl)  actfon  was  taken  at  the  M 
period,  toe  possible  actions  must  he  considered  at  the  second  period  -  evacuate,  issue  Motteflotd 
watch  do  nothtog  Depending  on  the  acUons  taken  in  the  first  and  the  second  periods  and  on  the  wara  flow 
rt*e^o“  period,  Vint  values  of  the  expected  losses  tor  e^h  of  *e  KnnM  chance  n^s  are 
calculated  Three  equally  probable  underlying  pdfs  for  the  water  flow  prevail  m  the  first  period.  At  Ae  end 
of  Ae  first  period,  after  measurmg  Ae  water  flow  Ae  posterior  probabilities  for  each  of  Aese  pdfs  are 

calculated  using  Bayes'  formula: 


prob.  (Wj  1  LN)  prob.  (LN^ 
prob.  (Lnj  |  Wj)  =  -j- 

^  prob.iWj  1  LN)  prob.  (LN^ 


(2.15) 
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(1.4;  5,600.000 


FLOOD 


(0.529;  3,611.663) 
(0.3058;  1,264.579) 


D2  * 

WA2  ^ 

V  HIGHER 

K  ' 

(0.1529;611,663) 

(0.3050:764,579) 


DN2  (0.4588:1,376.241) 


vSAME 


tO.7;  2,800.00 

^  — •  (0;0) 

p  (1.4;  3.500,000 

^  (0;0) 

6,300.000 

-•(0;0) 

p  (0.7;  2,000, 00( 


Ey2^ 

WA2 


^(0;0) 

(1.4;  3.500.000 


LOWERS 


p  (2.1:6,300,000 

^  (0:0) 
p  (0.7;  2,800.000 


D4  _ Ten 


(0;0) 

(1.4:3.500,000 

(0;0) 

(2.1;  6,300,000 


Figure  2-5.  Second  Stage  Tree  Corresponding  to  Chance  Node  C2  (Discrete  Case) 


Risk-Based  Evaluation  of 


where  prob.(LNi)  =  1/3  and  Wj  is  given  in  assumption  e,  and  prob.(V(  ]  LI^  )  is  calculated  using  Eqs. 
(2.12),  (2.13),  and  (2.14).  Then,  the  probability  of  a  flood  event  at  any  chance  node  is  calculated  as 

3 

prob.  (flood  1  Wj)  =£  prob.  (Flood  |  LN^  prob.  (LN,  |  Wj)  2.16 


For  example, 

prob.(flood  I  higher)  =  prob.(flood  |  LN.)  *  prob.(LN,  |  higher)  +  prob.(flood  |  LN^)  * 
prob.(LN2  I  higher)  +  prob. (flood  |  LN,)  *  prob.(LN3  |  higher) 

The  values  of  prob.(flood  j  LN^)  (i=  1,  2,  3)  are  calculated  using  Eqs.  (2.12),  (2.13)  and  (214)  ^d  the 
values  of  Pr(LN,  |  higher)  (i  =  1,  2,  3)  are  calculated  using  Eq.  (15).  Therefore,  from  Eq.  (2.16), 

prob. (flood  I  higher) 

=  0.3372  *  0.603129  +  0.0427  *  0.351699  +  0.0013  *  0.045172  =  0.218451 


Similarly, 

prob.(flood  1  same)  =  0.100545 
prob. (flood  1  lower)  =  0.021444 

The  required  value  of  the  loss  vector-valued  functions  is  then  computed  by  multiplying  the  flood  probability 
by  the  damage  vector.  Consider,  for  example,  arc  EV2  corresponding  to  decision  node  D2  in  Fig.  2-5, 


Lbv3  1 02  =  0.21845 1  *0.7  =  0.1529 


Cbv2|d2  =  0.218451  *  2,800,000  -I-  3,000,000  =  3,611,663 


Table  2-1  presents  the  values  of  the  loss  vectors  for  the  second-period  decision  arcs.  Folding  back 
at  each  decision  node,  the  vector-valued  functions  are  compared,  and  all  dominated  (inferior)  solutms  are 
eliminated.  Consider,  for  example,  decision  node  D2.  The  vector  corresponding  to  the  decision  DN2  is 
inferior  to  the  vector  corresponding  to  the  decision  WA2. 


0.3058 

1,264,579 


W.42 


< 


0.4588 

1,376,241 


DN2 


Table  2-2  presents  the  noninferior  decisions  for  the  second-period  decision  arcs.  Averaging-out  at 
the  chance  nodes  for  the  first  period,  each  noninferior  decision  corresponding  to  each  arc  is  multiplied  by 
the  probability  for  that  arc,  yielding  a  single  decision  rule  for  the  first-  period  decision  node.  For  example, 
we  have  18  different  combinations  at  WAl,  one  of  which  is  (EV2  ]  higher,  EV2  |  same,  EV2  j  lower). 
The  value  of  the  loss  vector  for  this  combination  is: 
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Table  2-1.  Expected  Value  of  Loss  Vectors  for  the  Second-Period  Decision  Arcs 

(Discrete  Case) 


Node 

Arc 

L 

c 

D2 

EV2 

0.1529 

3,611,663 

VA2 

0.3058 

1,264,579 

DN2 

0.4588 

1,376,241 

D3 

EV2 

0.0704 

3,281,526 

WA2 

0.1408 

851,908 

DN2 

0.2112 

633,434 

D4 

EV2 

0.0150 

3,060,043 

WA2 

0.0300 

575,054 

DN2 

0.0450 

135,097 

D5 

EV2 

0.3058 

3,917,494 

WA2 

!  0.4588 

1,570,410 

DN2 

1.5292 

1,529,157 

D6 

EV2 

0.1408 

3,422,289 

WA2 

0.2112 

992,671 

DN2 

0.7038 

703,815 

D7 

EV2 

0.0300 

3,090,065 

WA2 

0.0450 

605,076 

DN2 

0.1501 

150,108 

C2 

F 

0.1779 

711,760 

C3 

F  1 

0.8897 

889,700 

L  —  loss  of  lives 
C  -  cost  ($) 

Table  2-2.  Noninferior  Decisions  for  the  Second-Period  Decision  Nodes  (Discrete  Case) 


Node 

Noninferior  decisions 

D2 

EV2,  WA2 

D3 

EV2,  WA2,  DN2 

D4 

EV2,  WA2,  DN2 

D5 

EV2,  WA2,  DN2 

D6 

EV2,  WA2,  DN2 

D7 

EV2,  WA2,  DN2 

52 
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0.1779  +  0.1529  *  0.2466  +  0.0704  *  0.2686  +  0.0150  ♦  0.3577 
(1,000  +  711.76  +  3,611.633  *  0.2466  +  3,281.526  *  0.2686  +  3,060.043  *  0.3577)1,000 

0.2399  ■ 

"  4,578,391 

where  the  first  and  second  elements  represent  a  loss  of  lives  of  0.2399  and  a  cost  of  $4,578,391, 
respectively.  Table  2-3  presents  the  values  of  the  vector  of  objectives  for  the  first-period  decision  node. 
Note  from  Table  2-3  that  a  total  of  nine  noninferior  decisions  are  generated  for  action  WAl.  Similarly, 
there  are  five  noninferior  solutions  for  action  DNl  (by  self-comparison  of  all  vectors  for  that  action),  ^d 
fourteen  noninferior  solutions  after  comparing  all  decisions  for  all  actions  (see  Fig.  2-6).  Fig.  2-7  depicts 
the  graph  of  all  noninferior  solutions. 

Extension  to  Multiple-Risk  Measures 

Multiobjective  decision-tree  analysis  calls  for  the  adoption  of  multiple-risk  measures.  Often,  the 
expected  value,  by  itself,  provides  insufficient  information  for  risk  management.  The  expected  value  of 
adverse  effects,  which  has  been  most  commonly  used  in  conventional  decision-tree  analysis,  is  in  many 
cases  inadequate,  since  this  scalar  representation  of  risk  commensurates  events  that  correspond  to  all  levels 
of  losses  and  to  their  associated  probabilities.  The  common  expected-value  approach  is  particularly  deficient 
for  addressing  extreme  events,  since  these  events  are  concealed  during  the  amalgamation  of  events  of  low 
probability  and  high  consequence  with  events  of  high  probability  and  low  consequence.  The  synthesis  of 
several  approaches  -  single-objective  decision-tree  analysis,  multiobjective  optimization,  the  partitioned 
multiobjective  risk  method  (PMRM),  and  the  statistics  of  extremes  -  has  led  to  the  development  of  (he 
multiobjective  decision-tree  method.  This  new  form  of  decision-tree  analysis  can  handle  different  risk 
functions,  including  the  common  expected  value,  the  conditional  expected  value  for  extreme  events,  and 
the  event  with  maximum  probability,  thus  providing  decisionmaker(s)  with  more  comprehensive  knowledge 
and  a  robust  decision  policy. 

Determining  the  folding-back  strategy  associated  with  conditional  expected  values  is  substantially 
different  from  such  an  operation  using  the  conventional  expected  value.  Unlike  the  latter,  which  is  a  linear 
operator,  the  conditional  expected-value  operator  is  nonlinear.  This  nonlinearity  represents  an  obstacle  in 
decomposing  the  overall  value  of  the  conditional  expected  value  and  in  calculating  it  at  ^fferent  decision 
nodes.  Thus,  in  calculating  conditional  risk  functions  fi,  all  performance  measures  at  the  different  branches 
are  mapped  to  the  terminal  points  where  the  partitioning  is  performed. 

In  order  to  develop  a  folding-back  strategy  for  the  conditional  expected  value  f4  (the  schemes  for 
fj  and  fj  are  similar  and  thus  are  omitted  here),  some  properties  in  a  sequential  calculation  of  f4  will  be  first 
discussed. 

Consider  a  two-stage  decision-tree  problem  with  a  damage  function  f(a,,6i,a2,02),  where  Uj  is  the 
action  at  stage  j  and  6j  is  the  state  of  nature  at  stage  j  (j  =  1  and  2).  The  optimal  value  of  f4  is  given  by 
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Figure  2-6.  Decision  Tree  for  the  First  Stage  (Discrete  Case) 
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Figure  2-7.  Pareto  Optimal  Frontier  (Discrete  Case) 
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Table  2-3.  Decisions  for  the  First-Period  Decision  Node  (Discrete  Case) 


First- 

period 

decision 

1  Second-period  decision 

1  Loss  vector 

Higher 

Same 

Lower 

L 

C 

*  EVl 

- 

- 

- 

0.0000 

5,177,940 

•*  WAl 

EV2 

EV2 

EV2 

0.2399 

4,578,391 

*  WAi 

EV2 

EV2 

WA2 

0.2452 

3,689,511 

*  WAl 

EV2 

EV2 

DN2 

0.2506 

3,532,138 

WAl 

EV2 

WA2 

EV2 

CO 

oo 

m 

CM 

o 

3,925,796 

*  WAl 

EV2 

WA2 

WA2 

0.2641 

3,036,916 

*  WAl 

EV2 

WA2 

DN2 

0.2695 

2,879,543 

WAl 

EV2 

DN2 

EV2 

0.2777 

3,867,113 

WAl 

EV2 

DN2 

WA2 

0.2830 

2,978,233 

*  WAl 

EV2 

DN2 

DN2 

0.2884 

2,820,860 

WAl 

WA2 

EV2 

EV2 

0.2776 

3,999,600 

WAl 

WA2 

EV2 

WA2 

0.2829 

3,110,720 

WAl 

WA2 

EV2 

DN2 

0.2883 

2,953,347 

WAl 

WA2 

WA2 

EV2 

0,2965 

3,347,005 

*  WAl 

WA2 

WA2 

WA2 

0.3018 

2,458,125 

*  WAl 

WA2 

WA2 

DN2 

0.3072 

2,300,752 

WAl 

WA2 

DN2 

EV2 

0.3154 

3,288.323 

WAl 

WA2 

DN2 

WA2 

0.3207 

2,399,442 

*  WAl 

WA2 

DN2 

DN2 

0.3261 

2,242,070 

DNl 

EV2 

EV2 

EV2 

1.0136 

3,880,297 

DNl 

EV2 

EV2 

WA2 

1.0190 

2,991,417 

DNl 

EV2 

EV2 

DN2 

1.0566 

2,828,675 

DNl 

EV2 

'  WA2 

1  EV2 

1.0325 

3,^7,701 

*  DNl 

EV2 

WA2 

WA2 

1.0379 

2,338,821 

DNl 

EV2 

WA2 

DN2 

1.0756 

2,176,079 

DNl 

EV2 

DN2 

EV2 

1.1648 

3,150,115 

DNl 

EV2 

DN2 

WA2 

1.1702 

2,261,235 

DNl 

EV2 

DN2 

DN2 

1.2078 

2,098,493 

DNl 

WA2 

EV2 

EV2 

1.0513 

3,301.506 

DNl 

WA2 

EV2 

WA2 

1.0567 

2,412,626 

DNl 

WA2 

EV2 

DN2 

1.0943 

2,249,884 

DNl 

WA2 

WA2 

EV2 

1.0702 

2,648,910 

*  DNl 

W/V2 

WA2 

WA2 

1.0756 

1,760,030 

*  DNl 

WA2 

WA2 

DN2 

1.1132 

1,597,288 

DNl 

WA2 

DN2 

EV2 

1.2025 

2,571,324 

DNl 

WA2 

DN2 

WA2 

1.2079 

1,682,444 

*  DNl 

WA2 

DN2 

DN2 

1.2455 

1.519,702 

DNl 

DN2 

EV2 

EV2 

1.3153 

3,291,333 

DNl 

DN2 

EV2 

WA2 

1.3207 

2,402,453 

DNl 

DN2 

EV2 

DN2 

1.3583 

2,239,711 

DNl 

DN2 

WA2 

EV2 

1.3342 

2.638.737 

DNl 

DN2 

WA2 

WA2 

1.3396 

1,749,857 

DNl 

DN2 

WA2 

DN2 

1.3772 

1,587,115 

DNl 

DN2 

DN2 

EV2 

1.4665 

2,561,151 

DNl 

DN2 

DN2 

WA2 

1.4719 

1,672,271 

*  DNl 

DN2 

DN2 

DN2 

1.5095 

1,509,529 

*  -  noninferior  decisions 
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i  J  f(ai,0i,a2,02)P(01.02  I  ai,a2)deide2 

.*  _  min  f(ai,0i,a2,02)^P~^Ca) _ _ _ _ _ 

4  “ai,a2  J  T  p(0  1 ,02  I  ai  ,a2)d0  id02 

f(ai,0i,a2,02)^P-ka) 

(2.17) 


where  a  is  the  partitioning  point  on  the  probability  axis.  Rewrite 

p(0„03  I  =  p(02  1  0.,a„a3)p(0,  |  a^)  (2-18) 

The  fact  that  an  action  at  a  subsequent  stage  does  not  affect  the  state  of  nature  at  a  previous  stage  is  used 
in  Eq.  (2.18).  Consequently,  the  optimization  problem  in  Eq.  (2.17)  can  be  evaluated  in  a  two-stage  form. 


1  [  J  f(ai,0i,a2,02)p(92  I  ®l»ai,a2)d02  ]  p(0i  I  ai)d0i 

.*  _  min  f(ai  ,0i  ,a^02)^P'na) _ _ _ 

4"  ai,a2  7  i  J  p(02  I  0  I,ai,a2)d02  ]  p(0l  •  ai)d0i 

f(ai  ,0 1  ,a2,02)^P'  *  (a) 


(2.19) 

The  optimization  problem  in  Eq.  (2.19)  is  nonseparable.  To  separate  the  objective  function  with  respect  to 
stages,  it  is  thus  necessary  to  record  two  numbers  at  each  stage  --  the  values  of  the  numerator  and  the 
denominator  for  each  optimal  conditional  expected  value.  A  more  serious  problem  related  to  the 
decomposition  of  Eq.  (2.19)  is  its  nonmonotonicity.  This  can  be  easily  observed  by  the  fact  that 
minimization  of  a(*)/b(«)  does  not  necessarily  lead  to  the  solution  of  minimization  of  [c+a(«)]/[d-l-b(«)] 
where  c  and  d  are  two  constants  and  b  and  d  are  positive.  The  only  exception  to  the  above  claim  holding 
is  the  case  where  b  remains  a  constant.  The  following  simplification  will  be  introduced  to  make  stagewise 
calculation  of  the  value  of  the  conditional  expectation  f4  possible.  From  the  definition,  we  have  P[f(0i,02) 
^  P'(a)]  =  a.  When  the  value  of  0,  is  fixed,  P[f(0i,02)  ^  P’(a)  |  0,]  is  not  necessarily  equal  to  a.  In  order 
to  have  a  common  denominator,  we  introduce  a  set  of  PfXa)  to  keep  P[f(0i,02)  s  Pr‘(a)  1  0i]  =  a,  where 
Pi  is  the  conditional  cumulative  distribution  function  of  02,  given  value  of  0,.  When  we  fold  back,  this 
simplification  yields 

/  P[f(0.,02)  Pf*  (a)  I  0,]  p  (0,)d0i  =  a 

e  (z.zu; 


In  summary,  we  should  adhere  to  the  following  rules  when  calculating  the  conditional  expected  value  in 
the  folding-back  procedure  of  decision  trees: 

1)  Partition  and  calculate  f4  at  terminal  points  according  to  the  conditional  probability  density 
function. 
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2)  Fold  back  and  perform  at  each  chance  node  the  operation  of  the  expected  value. 

Note  that  although  reducing  the  variance  (the  uncertainty)  of  the  risk  may  not  contribute  much  to 
reducing  the  expected  value  fj,  it  often  markedly  reduces  the  conditional  expected  value  f4  associated  with 
extreme  events  (see  Fig.  2-8).  Two  benefits  that  result  from  additional  experimentation  include  reducing 
the  expected  loss  and  reducing  the  uncertainly  associated  with  decisionmaking  under  risk.  However,  in  most 
cases,  these  two  dual  aspects  of  experimentation  conflict  with  each  other.  The  general  framework  of 
multiobjective  decision-tree  analysis  proposed  here  provides  a  medium  with  which  these  dual  aspects  can 
be  captured  by  investigating  the  multiple  impacts  of  experimentation. 


Example  Problem  for  the  Continuous  Case 
Problem  Definition 


The  problem  developed  in  the  previous  example  for  the  discrete  case  is  modified  here  to  handle 
continuous  loss  functions  and  extreme  random  events.  The  main  diflference  between  the  discrete  and  the 
continuous  cases  lies  in  calculating  the  damage  vector  for  the  terminal  nodes,  which  can  be  determined 
using  the  expected  value  f5(*)  and/or  the  conditional  expected  value  ft(*).  The  subsequent  computations 
are  similar  to  those  carried  out  for  the  discrete  case.  Consequently,  assumption  f  for  the  discrete  case  is 
modified  as: 

I 

f)  L  and  C  are,  respectively,  the  possible  loss  of  lives  and  cost,  given  that  no  flood  warning  is 
issued;  they  are  linear  functions  of  the  water  flow  W.  All  other  costs  (as  shown  in  Fig.  2-4)  are 
given  in  terms  of  the  loss  functions  L  and  C,  where  L  =  W  *  Lp,  and  Lp  =  0.0001,  C  =  W  *  Cp, 
and  Cp  =  100 


The  complete  decision  tree  for  this  case  is  shown  in  Fig.  2-9.  The  loss  functions  L  and  C  are  calculated 
using  the  unconditional  expected-value  function  f5(*)  and/or  the  conditional  expected-value  function  f  4(«). 
The  unconditional  expected  loss  f5(*)  is  given  by 


fsC) 


50000 


1 

-exp 


1 

'  ln(T^0  -  fY 

2 

o 

dW 


1  -  3) 


10.82  -  IX. 


exp(^  +  0^/2) 


1  -  3» 


10.82 


(2.21) 


where  Pf  is  equal  to  L  f  or  C  f  when  Eq.  (2.21)  is  used  to  calculate  §  for  loss  of  lives  or  monetary  costs, 
respectively.  The  conditional  expected  loss  f4(*)  is  given  by 


/4(*) 


P,  [l  -  4>  [«■'(«)  -  o) 


] 


exp(/i  -1-  o^/2)/(l  -  a) 


(2.22) 
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Pdf 


Figure  2-8.  Variance  and  Region  of  Extreme  Events 
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where  Pf  is  equal  to  L  f  or  C  f  when  Eq.  (2.22)  is  used  to  calculate  for  loss  of  lives  or  monetary  costs, 
respectively,  and  a  is  the  partitioning  point  on  the  probability  axis,  which  is  0.99  in  this  case.  With  the  use 
of  Eqs.  (2.21)  and  (2.22),  the  cost  (C)  and  the  loss  of  lives  (L)  are  calculated  using  f4(«)  and  f5(*)  at  all 
the  terminal  nodes  for  each  of  the  decision  arcs.  Note  that  each  of  the  risk  functions  f4(«)  and  ^(•)  is 
composed  of  two  components  —  cost  and  loss  of  lives. 

Calculation  of  the  Loss  Vectors  For  the  First  Period 

Chance  node  Cl 

Assuming  that  the  possible  pdf  of  the  water  flow  (W)  is  LNj  with  probability  1/3,  i  =  1,2,3  and 
that  the  flood  stage  is  at  W  =  50,000  cfs,  two  outcomes  are  considered  at  the  end  of  the  second  period:  a 
flood  or  no  flood  event  (see  Fig.  2-5).  The  values  of  the  components  of  f4(«)  and  f5(*)  for  node  Cl  are 
calculated  using  Eqs.  (2.21)  and  (2.22),  respectively.  The  value  of  the  loss  vector  for  Cl  using  f4(*)  is 
shown  in  Fig.  2-10. 


Chance  nodes  C2  and  C3 

Four  possible  outcomes  at  the  beginning  of  the  second  period  are  investigated  at  nodes  C2  and  C3: 
a  flood  event,  the  water  flow  is  higher,  the  water  flow  is  the  same,  and  the  water  flow  is  lower  (see  Fig. 
2-20).  Similarly  to  the  discrete  case,  the  probabilities  of  these  outcomes  are  calculated  using  Eqs.  (2.12), 
(2.13),  and  (2.14). 


Calculation  of  Loss  Vectors  for  Second  Period 


Regardless  of  whether  a  watch  (WAl)  was  issued  or  a  do-nothing  (DN2)  action  was  followed  at 
the  first  period,  the  same  three  possible  actions  are  evaluated  at  the  second  period:  evacuate,  issue  another 
flood  watch,  or  do  nothing.  Depending  on  the  actions  taken  at  the  first  and  second  periods  and  the  water 
flow  level  at  the  second  period,  different  values  of  losses  are  generated  for  each  terminal  chance  node. 
There  are  three  equally  probable  underlying  pdfs  for  the  water  flow  for  the  first  period.  After  measuring 
the  water  flow  Wj  at  the  end  of  the  first  period,  the  posterior  probabilities  are  calculated  using  Eq.  (2.15). 
The  required  value  of  the  loss  vector  [of  f4(«)  and  f5(*)]  is  then  calculated  using  Eqs.  (2.21)  and  (2.23)  for 
f5(*)  and  Eqs.  (2.22)  and  (2.24)  for  f4(*): 

3 


f5(«  1 

II 

[f5(*) 

1  LNi]  Pr(LNi  1  Wj)) 

(2.23) 

U*  1 

II 

[f4(*) 

1  LNJ  Pr(LNi  1  Wj)) 

(2.24) 

For  example. 

f4(*  1 

higher)  = 

f4(*  1 

LN,)  *  prob.(LNi  |  higher)  -1-  f4(«  1  LN2)  * 

prob.(LN2  1  higher)  H-  f4(«  |  LN3)  *  prob.(LN3  1  higher) 
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(0;  88,574) 

EVI 

(O.OL;  0.20 

Figure  2-10.  Averaging  out  Chance  Node  Cl  Using  (Continuous  Case) 

The  values  of  prob.(LNi  |  higher)  (i=l,  2,  3)  are  calculated  using  Eq.  (2.15)  and  the  values  of 
f4(*  I  higher)  are  calculated  using  Eq.  (2.22).  Therefore,  Eq.  (2.24)  yields, 

f4(*  I  higher) 

=  500,697.95*0.603129  +  136.456.11*0.351699  +  37,188.63*0.045172 
=  351,657.32  (2.25) 

The  values  for  f4(«  |  same),  f4(«  |  lower),  f5(»  |  higher),  f5(*  |  same),  and  f5(*  |  lower)  are  calculated 
in  a  similar  way.  The  loss  vector  is  then  computed  by  multiplying  these  results  by  the  ratio  to  the  maximum 
damage  and  Lf  or  Cf,  as  the  case  may  be.  For  example,  the  components  of  the  loss  vectors  for  arc  EV2 
corresponding  to  decision  node  D2  are 

^Eva  I  D2,f4(*)  “  35.1657  *  0.1  =  3.5166 

Cev2  I  D2,f4(*)  =  35,165,732  *  0.4  +  3,000,000  =  17,066,293  (2.26) 

Table  2-4  summarizes  the  values  of  the  loss  vectors  f5(*)  and  f4(*)  for  the  decision  arcs  corresponding  to 
the  second  period.  Once  these  values  are  calculated,  the  noninferior  decisions  for  each  node  are  calculated 
by  folding  back  the  same  way  as  in  the  discrete  case.  Table  2-5  yields  the  noninferior  decisions  for  the 
second-period  decision  arcs.  Averaging-out  at  the  chance  nodes  for  the  first  period  follows  the  same 
procedure  used  in  the  discrete  case.  Consider,  for  example,  action  WAl.  There  are  27  different 
combinations  when  using  the  expected  value  f5(*),  and  four  different  combinations  when  using  f4(*).  Table 
2-6  yields  the  values  of  the  loss 
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(0.2L;  0.8C) 


FLOOD 


10.0798;  3.319,034) 
(0.IS9S;  696.792) 
(0.2393;  717,6  26) 


(0.0798;  319,034) 


(0.1595;  398,792) 


(0.1  L;  0.4C) 


HIGHER 


LOWER' 


(0.2393;  717,826) 


(0.2L:  0.5C) 


(0.3L;  0.9C) 


(0.1  L;  0.4C) 


(0.2L:  0.5C) 


(0.3L:  0.9C) 


(0.1  L;  0.4C) 


/  WA2  /L,.. 
D4  - - - (C11 


(0.2L;  0.5C) 


(0.3L;  0.9C) 


Figure  2-11.  Second  Stage  Corresponding  to  Chance  Node  Cw  Using  fs  (Continuous  Case) 
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Table  2-4.  Loss  Vectors  for  the  Second-period  Decision  Arcs  (Continuous  Case) 


Node 

Arc 

f 

5<*> 

^4 

(•) 

L 

C 

L 

C 

D2 

•  *  EV2 

0.0798 

3,319,034 

3.5166 

17,066,293 

•  WA2 

0.1595 

898,792 

7.0332 

18,082,866 

•  .  DN2 

0.2393 

717,826 

10.5497 

31,649,159 

D3 

•  *  EV2 

0.0312 

3,124,816 

1.9595 

10,838,047 

•  *  WA2 

0.0624 

656,020 

3.9190 

10,297,559 

•  DN2 

0.0936 

280,835 

5.8785 

17,653,606 

D4 

•  *  EV2 

0.0040 

3,016,172 

0.7598 

6,039,120 

•  *  WA2 

0.0.081 

520,215 

1.5196 

4,298,901 

•  DN2 

0.0121 

36,387 

2.2793 

6,838,021 

D5 

•  *  EV2 

0.1595 

3,478,550 

7.0332 

24,099,439 

•  WA2 

0.2393 

1,058,309 

10.5497 

25,116,012 

•  DN2 

0.7976 

797,584 

35.1657 

35,165,732 

D6 

•  *  EV2 

0.0624 

3,187,223 

3.9190 

14,757,071 

•  *  WA2 

0.0936 

718,427 

5.8785 

14,216,583 

•  DN2 

0.3120 

312,039 

19.5951 

19,595,118 

D7 

•  *  EV2 

0.0081 

3,024,258 

1.5196 

7,558,681 

•  *  WA2 

0.0121 

528,301 

2.2793 

5,818,461 

•  DN2 

0.0404 

40,430 

7.5978 

7,597,801 

C2 

F 

0.0886 

354,298 

4.4956 

17,982,472 

C3 

F 

0.4429 

442,872 

22.4781 

22,478,090 

•  noninferior  decisions  using 

*  noninferior  decisions  using  f^(*) 


Table  2-5.  Noninferior  Decisions  for  the  Second-period  Decision  Nodes  (Continuous  Case) 


Node 

!  Noninferior  decisions  | 

f5(.) 

D2 

EV2,  WA2.  DN2 

EV2 

D3 

EV2,  WA2,  DN2 

EV2,  WA2 

D4 

EV2,  WA2,  DN2 

EV2,  WA2 

D5 

EV2,  WA2,  DN2 

EV2 

D6 

EV2.  WA2.  DN2 

1  EV2 ,  WA2 

D7 

EV2,  WA2,  DN2 

EV2 ,  WA2 
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vectors  for  the  first-period  decision  node  using  and  Table  2-7  yields  the  values  of  the  loss  vectors 
using  f4(*).  Note  from  Table  2-8  that  for  action  WAl  there  are  a  total  of  10  noninferior  decisions  by  self¬ 
comparison.  Similarly  for  action  DNl,  there  are  8  noninferior  solutions  by  self-comparistm  of  all  vectors 
for  action  DNl,  and  6  noninferior  solutions  after  comparison  of  all  decisions  for  all  actions  using  fsC*)  (see 
Fig.  2-12).  Figure  2-13  depicts  the  graph  of  all  noninferior  solutions  using  fj(*).  Note  from  Table  2-8  that 
there  is  only  one  noninferior  action.  The  action  EVl  yields  the  most  conservative  action  from  the  point  of 
view  of  extreme  events.  When  the  decisionmaker  considers  the  possible  extreme  event,  the  potential  loss 
of  property  dominates  the  cost  of  the  warning  system.  Thus,  the  two  objective  functions  do  not  conflict  at 
this  case. 


Conclusions 

Multiobjective  decision-tree  analysis  is  an  extension  of  the  single-objective-based  (tecision-tree 
analysis  formally  introduced  more  than  two  decades  ago  by  Howard  Raiffa  [1968].  This  extension  is  made 
possible  by  making  a  synthesis  of  the  traditional  method  and  more  recently  developed  approaches  used  for 
multiobjective  analysis  and  for  the  risk  of  extreme  and  catastrophic  events.  Success^  applications  of  single¬ 
objective  decision-tree  analysis  to  numerous  business,  engineering,  and  governmental  decisionmaking 
problems  over  the  years  have  made  the  methodology  into  an  important  and  valuable  tool  in  systems 
analysis.  Its  extension  -  incorporating  multiple  noncommensurate  objectives,  impact  analysis,  and  the 
conditional  expected  value  for  extreme  and  catastrophic  events  -  might  be  viewed  as  an  indicator  of  growth 
in  the  broader  field  of  systems  analysis  and  in  decisionmaking  under  risk  and  uncertainty.  Undoubtedly, 
there  remain  several  theoretical  challenges  that  must  be  addressed  to  fully  realize  the  strengths  and 
usefulness  of  the  extended  methodology.  In  this  sense,  the  multiobjective  decisiem-tree  analysis  proposed 
here  constitutes  the  first,  albeit  important,  step  in  the  direction  of  developing  improved  and  more 
representative  models  and  decisionmaking  tools. 
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Table  2-6.  Decisions  for  the  First-period  Node  Using  fs  (Continuous  Case) 


First- 


*  EVl 

*  WAl 

*  WAl 


Second-period  decision 


Loss  vector 


C 


0.0408 

0.0423 


*  WAl 

EV2 

EV2 

DN2 

0.0437 

WAl 

EV2 

WA2 

EV2 

0.0492 

*  WAl 

EV2 

WA2 

WA2 

0.0507 

*  WAl 

EV2 

WA2 

DN2 

0.0521 

WAl 

F.V2 

DN2 

EV2 

0.0575 

WAl 

EV2 

DN2 

WA2 

0.0590 

*  WAl 

EV2 

DN2 

DN2 

0.0604 

WAl 

WA2 

EV2 

EV2 

0.0604 

WAl 

WA2 

EV2 

WA2 

0.0619 

WAl 

WA2 

EV2 

DN2 

0.0633 

WAl 

WA2 

WA2 

EV2 

0.0688 

*  WAl 

WA2 

WA2 

WA2 

0.0703 

*  WAl 

WA2 

WA2 

DN2 

0.0717 

WAl 

WA2 

DN2 

EV2 

0.0771 

WAl 

WA2 

DN2 

WA2 

0.0786 

*  WAl 

WA2 

DN2 

DN2 

0.0800 

WAl 

DN2 

EV2 

EV2 

0.0801 

WAl 

DN2 

EV2 

WA2 

0.0816 

WAl 

DN2 

EV2 

DN2 

0.0830 

WAl 

DN2 

WA2 

EV2 

0.0885 

WAl 

DN2 

WA2 

WA2 

0.0900 

WAl 

DN2 

WA2 

DN2 

0.0914 

WAl 

DN2 

DN2 

EV2 

0.0968 

WAl 

DN2 

DN2 

WA2 

0.0983 

*  WAl 

DN2 

DN2 

DN2 

0.0997 

DNl 

EV2 

EV2 

EV2 

0.1153 

DNl 

EV2 

EV2 

WA2 

0.1167 

DNl 

EV2 

EV2 

DN2 

0.1270 

DNl 

EV2 

WA2 

EV2 

0.1236 

*  DNl 

EV2 

WA2 

WA2 

0.1250 

*  DNl 

EV2 

WA2 

DN2 

0.1353 

DNl 

EV2 

DN2 

EV2 

0.1823 

DNl 

EV2 

DN2 

WA2 

0.1837 

DNl 

EV2 

ON2 

DN2 

0.1940 

DNl 

WA2 

EV2 

EV2 

0.1350 

DNl 

WA2 

EV2 

WA2 

0.1364 

DNl 

WA2 

EV2 

DN2 

0.1467 

5,088,574 

3,781,716 

2,888,912 

2,715,847 

3,118,597 

2,225,793 

2,052,728 

3,017,822 

2,125,018 

1,951,953 

3,184,884 

2,292,080 

2,119,015. 

2,521,765 

1,628,961. 

1,455,896 

2,420,990 

1,528,186 

1,355,121 

3,140,258 

2,247,454 

2,074,389 

2,477,139 

1,584,335 

1,411,270 

2,376,364 

1,483,560 

1,310,495 

2,851,964 

1,959,160 

1,784,649 

2,188,846 

1,296,042 

1,121,531 

2,079,690 

1,186,886 

1,012,375 

2,255,133 

1,362,329 

1,187,818 
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First 

Period 

Decision 

DNl 

*  DNl 

*  DNl 
DNl 
DNl 

*  DNl 
DNl 
DNl 
DNl 
DNl 
DNl 
DNl 
DNl 
DNl 

*  DNl 


Table  2-6  (continued) 

Second  period  decision _ I 


Higher 

Sane 

Lower 

WA2 

WA2 

EV2 

WA2 

WA2 

WA2 

WA2 

WA2 

DN2 

WA2 

DN2 

EV2 

WA2 

DN2 

WA2 

WA2 

DN2 

DN2 

DN2 

EV2 

EV2 

DN2 

EV2 

WA2 

DN2 

EV2 

DN2 

DN2 

WA2 

EV2 

DN2 

WA2 

WA2 

DN2 

WA2 

DN2 

DN2 

DN2 

EV2 

DN2 

DN2 

WA2 

DN2 

DN2 

DN2 

Loss  vector 
I  C 


0.1433 

0.1447 

0.1550 

0.2020 

0.2034 

0.2137 

0.2727 

0.2741 

0.2844 

0.2810 

0.2824 

0.2927 

0.3397 

0.3411 

0.3514 


1,592,015 

699,211 

524,700 

1.482,859 

590,055 

415,544 

2,190,838 

1,298,034 

1,123,523 

1,527,720 

634,916 

460,405 

1,418,564 

525,760 

351,249 


*  noninferior  decisions 

Table  2-7.  Decisions  for  the  First-period  Decision  Node  Using  f4  (Continuous  Case) 


First- 

period 

decision 

*  EVl 
WAl 
WAl 
WAl 
WAl 
DNl 
DNl 
DNl 
DNl 


Second-period  decision 
Higher  I  Same  Lower 


Loss  vector 

I  c 


0.0000  9,495,618 
EV2  2.2367  12,565,412 
WA2  2.5085  11,942,936 
EV2  2.7630  12,420,237 
WA2  3.0348  11,797.761 
EV2  6.1876  15,467,376 
WA2  6.4593  14,844,900 
EV2  6.7140  15,322,201 
WA2  6.9854  14,699,725 

-  Noninferior  decisions 
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Figure  2-12.  Decision  Tree  for  the  Second  Stage  Using  fj  (Continuous  Case) 
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Figure  2-13.  Pareto  Optimal  Frontier  Using  fs  (Continuous  Case) 
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Parts 

Performance  Characteristics  of  a 
Flood  Warning  System: 

Technical 

Introduction 

From  the  utilitarian  point  of  view,  rooted  in  the  Bayesian  principles  of  rationality,  the  ultimate 
measure  of  performance  of  a  flood  warning  system  is  the  ex  ante  economic  value.  From  the  engineering 
point  of  view,  there  remains  the  need  for  auxiliary  measures  that  characterize,  perhaps  only  partially,  the 
performance  of  various  components  of  a  flood  warning  system.  The  purpose  of  such  measures  is  to  aid  the 
engineer  in  the  process  of  planning  and  design. 

One  aspect  of  the  performance  of  a  flood  warning  system  is  its  reliability .  The  following  presents 
an  overview  of  a  model  that  outputs  two  measures  of  system  reliability: 

•  the  relative  operating  characteristic  (ROC),  which  shows  a  relationship  among  (i)  the 
probability  of  detection,  (ii)  the  probability  of  a  false  warning,  and  (iii)  the  expected  lead 
time  of  a  warning,  and 

•  the  performance  tradeoff  characteristic  (PTC),  which  shows  a  relationship  among  (i)  the 
expected  number  of  detections  per  year,  (ii)  the  expected  number  of  false  warnings  per 
year,  and  (iii)  the  expected  lead  time  of  a  warning. 

Each  characteristic,  the  ROC  and  PTC,  can  be  displayed  graphically  in  the  form  of  a  family  of  curves.  The 
displays  offer  an  aid  to  engineering  planning  and  design  of  flood  warning  systems.  The  concept  and 
interpretation  of  these  displays  are  illustrated  with  a  case  study  of  the  flood  warning  system  for  Milton, 
Pennsylvania. 

System  Model 

Structure 


The  model  is  tailored  to  a  class  of  local  warning  systems  which  can  be  conceptualized  as  a  cascade 
coupling  of  three  components,  shown  in  Figure  3-1 :  monitor,  forecaster,  and  decider .  The  operation  of  such 
a  system  is  idealized  as  follows. 

Floods  occur  intermittently.  For  economic  reasons,  a  flood  data  collection  network,  forecasting 
procedure,  and  emergency  management  do  not  operate  continuously.  Rather,  their  operation  is  triggered 
only  when  potential  flood  conditions  are  detected.  To  enable  such  detections,  a  system  monitoring 
hydrometeorologic  conditions  operates  continuously.  When  a  set  of  predefined  conditions  is  observed,  the 
monitor  triggers  operation  of  Ae  forecast  system.  The  flood  data  collection  network  is  activated,  and  a 
forecast  of  the  flood  hydrograph  is  prepared.  This  forecast  is  supplied  to  the  decision  system  --  a  flood 
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TRIGGER  ^ 

FORECAST. 

DECIDER 

WARNING  . 

rwiitlwAw  1  cn 

Figure  3-1.  Functional  Structure  of  a  Flood  Warning  System 

preparedness  organization,  or  a  floodplain  manager  --  who  must  then  decide  whether  or  not  to  issue  a 
warning  to  the  public. 

Assumptions 

Principal  definitions  and  assumptions  underlying  our  model  of  a  flood  warning  system  are  as 

follows: 


1.  A  flood  is  the  portion  of  a  hydrograph  above  a  flood  stage,  officially  specified  for  a  given  river 
gauging  station. 

2.  If  a  flood  forecast  is  prepared,  it  is  issued  at  a  well-defined  instant,  consistently  for  every  flood. 
The  performance  of  a  warning  system  is  evaluated  based  on  this  one  forecast. 

3.  The  decision  whether  or  not  to  issue  a  warning  to  the  public  is  based  on  the  forecasted  flood 

crest. 


4.  The  flood  plain  is  divided  into  elevation  zones.  A  flood  warning  is  issued  for  a  zone.  Thus, 
depending  on  the  forecast,  it  may  be  optimal  to  issue  a  warning  for  a  lower  zone,  but  not  for  an  upper  zone. 
Consequently,  the  performance  characteristics  are  defined  for  a  zone. 

Mathematical  models  of  the  three  system  components  are  described  below.  Section  3  defines  the 
performance  measures.  Section  4  illustrates  them  with  numerical  examples. 

Monitor 

An  all-important  design  decision  is  the  choice  of  a  forecast  trigger  ~  an  observable  state  that  is 
likely  to  precede  every  flood  and  that,  once  observed,  will  trigger  preparation  of  flood  forecasts.  Here  are 
three  examples  of  triggers: 

(river  stage)  >  (threshold) 

(rainfall  intensity  and  duration)  >  (threshold) 

(meteorologic  situation)  e  {potential  flood  situations} 
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A  good  monitor  is  particularly  critical  to  local  warning  systems  for  flash  floods  in  headwater  areas 
with  small  watersheds  and  short  concentration  times.  For  instance,  a  flood  developing  rapidly  during 
nighttime  may  occur  undetected  because  of  an  equipment  failure;  a  trigger  may  be  false  because  the 
oncoming  storm  suddenly  changes  its  track  and  bypasses  the  watershed.  Consequently,  the  performance  of 
the  monitor  may  limit  the  performance  of  the  total  warning  system,  no  matter  how  sophisticated  its  flood 
data  collection  network,  forecasting  procedure,  and  emergency  management. 

In  order  to  characterize  the  performance  of  a  monitor,  introduce  the  following  variables: 

T  ~  trigger  indicator:  trigger  is  not  observed  (T  =  0),  trigger  is  observed  (T  =  1); 

0  -  flood  indicator:  flood  does  not  occur  (0  =  0),  flood  occurs  (0  =  1). 

Next,  define  two  conditional  probabilities: 

Y  =  p(0  =  1  I  T  =  1)  ,  (3-1) 

p  =  P(T  =  1  I  0  =  1)  .  (3.2) 

Probability  y  characterizes  the  diagnosticity  of  the  monitor.  For  example,  y  =  1  means  that  every 
trigger  is  followed  by  a  flood;  in  other  words,  the  monitor  provides  a  perfect  diagnosis  of  a  flood  situation. 
Probability  p  characterizes  the  reliability  of  the  monitor.  For  example,  p  =  1  means  that  every  food  is 

preceded  by  a  trigger;  in  other  words,  the  monitor  never  fails  to  signal  the  oncoming  flood.  Probabilities 

Y  and  p  are  independent. 

Forecaster 

The  objective  of  modeling  is  to  obtain  a  stochastic  characterization  of  floods  and  forecasts  in  the 
form  requisite  for  decisionmaking  and  perfoimance  evaluation.  Toward  this  end,  a  Bayesian  processor  of 
forecasts  is  formulated  following  the  principles  laid  down  in  earlier  works  of  Krzysztofowicz  [1983a, 
1983b,  1985,  1987].  The  inputs  into  the  processor  are  a  prior  distribution  describing  natural  flood  events 
and  a  likelihood  function  describing  the  stochastic  dependence  between  forecasted  and  actual  flood  events. 
The  principal  output  from  the  processor  is  the  posterior  probability  of  flooding  a  given  zone  elevation, 
conditional  on  the  forecast.  This  probability  provides  a  basis  for  deciding  the  warning.  In  addition,  the 
processor  outputs  several  other  probability  distributions  needed  for  system  performance  evaluation.  The 
remainder  of  this  section  outlines  our  approach  to  modeling  the  prior  distribution  and  the  likelihood 
functions. 

Model  of  floods:  prior  distribution.  A  flood  is  described  in  terms  of  the  actual  flood  crest,  h, 
measured  from  the  flood  stage.  Before  a  forecast  is  prepared,  the  uncertainty  about  the  magnitude  of  &e 
flood  crest,  conditional  on  the  hypothesis  that  a  flood  will  occur,  0=  1,  is  described  in  terms  of  the  prior 
density: 


g(h  I  0  =  1). 


(3.3) 
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This  density  should  be  estimated  from  a  partial  duration  flood  series. 

Model  of  forecasts:  likelihood  function.  When,  and  only  when,  a  trigger  is  observed,  T  =  1,  the 
forecaster  is  activated,  hydrometeorologic  observations  are  collected,  and  a  flood  forecast  is  prepared.  A 
categorical  forecast  is  assumed  to  specify  a  point  estimate  s  of  h. 

If  a  flood  does  not  occur  after  the  forecast,  0=0,  then  the  forecasted  flood  crest  s  cannot  be 
verified.  Let 

Ko(s  I  0=  0,  T  =  1)  (3.4) 

denote  the  density  of  s  on  those  occasions.  If  a  flood  occurs  after  the  forecast,  0=1,  then  the  forecasted 
flood  crest  s  can  be  verified  against  the  actual  flood  crest  h.  Let 

f(s  I  h,  0  =  1,  T  =  1)  (3.5) 

denote  the  density  of  s,  conditional  on  h,  on  those  occasions. 

For  a  fixed  s,  functions  Ko(s  |  •)  and  f(s  |  •)  are  termed  the  likelihood  functions  of  the  respective 
events.  These  likelihood  functions  constitute  a  model  of  the  forecaster.  They  can  be  estimated  from  a  joint 
record  of  forecasted  and  actual  floods. 

Decider 


When  the  trigger  is  observed,  T  =  1 ,  and  the  forecast  s  of  the  flood  crest  is  prepared;  the  manager 
must  then  decide  whether  or  not  to  issue  a  flood  warning  for  a  zone  of  the  floodplain;  {w  =  0,  do  not  issue 
warning;  w  =  1,  issue  warning}.  Thereafter  the  event  takes  place:  {0  =  0,  zone  is  not  flooded;  6  =  1, 
zone  is  flooded}.  Each  decision-event  vector  (w,6)  leads  to  outcomes  whose  undesirability  (as  they  are 
mostly  losses  rather  than  gains)  is  evaluated  in  terms  of  a  disutility  function.  The  arguments  of  the  disutility 
function  are  the  actual  flood  crest  and  the  lead  time  of  the  warning  (to  be  defined  precisely  later). 

Let  W  denote  a  warning  rule  which  for  every  forecast  specifies  decision  w  =  W(s)  for  a  given 
zone.  The  objective  of  decision  analysis  is  to  find  the  optimal  ■warning  rule  W*.  According  to  the  Bayesian 
postulates  of  rationality,  the  rule  W*  should  minimize  the  posterior  expected  disutility  of  outcomes. 

For  a  statistical,  as  contrasted  with  the  economic,  evaluation  of  system  performance,  it  is  not 
necessary  to  find  the  exact  form  of  W*.  It  suffices  to  know  its  general  structure.  Under  certain  monotonicity 
conditions,  W*  is  of  the  threshold  type: 


W*(s)  = 


10  if  qis) 

[1  if  q(.s) 


< 

> 


q* 


(3.6) 


where  q(s)  =  P(0  =  1  |  s,  T  =  1)  is  the  posterior  probability  of  a  flood  in  a  given  zone,  and  q*  is  a 
threshold  dependent  upon  the  disutility  function  and  the  density  functions  (3.3)  -  (3.5).  The  optimal  warning 
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rule  states,  then,  that  a  warning  should  be  issued  for  a  given  zone  whenever  the  posterior  probability  q(s) 
of  flooding  that  zone,  conditional  on  forecast  s,  exceeds  threshold  q*. 

Performance  Measures 

Performance  Probabilities 


The  vector  (T,w,0,e)  of  binary  indicators  of  the  status  of  the  trigger  T,  warning  w,  flood  0,  and 
zone  flood  0  can  take  on  nine  values  which  define  four  performance  states  of  the  warning  system,  as  shown 
in  Figure  3-2.  The  states  are  as  follows: 

missed  flood:  M  =  (w  =  O|0=l,0=l) 

false  warning:  F  =  (w  =  1  |  0  =  0,  T  =  1) 

detection:  D  =  (w  =  1  |  0  =  1,  0  =  1) 

quiet:  Q  =  (w  =  0  |  0  =  0,  T  =  1)  (3.7) 

These  states  are  observable  in  the  sense  that  one  could  count  their  occurrences  over  a  period  of  time.  In  the 

limit,  this  count  would  give  rise  to  conditional  probabilities  of  incorrect  system  performance,  P(M)  and 
P(F),  and  correct  system  performance,  P(D)  and  P(Q). 

Since  P(M)  =  1  -  P(D)  and  P(Q)  =  1  -  P(F),  it  suffices  to  find  the  probability  of  detection,  P(D), 
and  the  probability  of  false  warning,  P(F).  The  objective  of  modeling,  then,  is  to  express  these  probabilities 
in  terms  of  parameters  and  functions  which  characterize  the  warning  system.  The  main  result  en  route  to 
deriving  such  expressions  is  the  following  factorizations: 


P(D)  =  pP(w  =  1  I  0  =  1,  0  =  1,  T  =  1)  (3.8) 

P(F)  =  (1  -  Y)P(w  =  1  I  0  =  0,  T  =  1)  yP(w  =  1  1  0  =  0,  0  =  1,  T  =  1)  (3.9) 

The  diagnosticity  y  and  reliability  p  characterize  the  monitor,  while  the  remaining  conditional  probabilities 
depend  upon  the  prior  distribution  (3.3),  the  likelihood  function  (3.4)-(3.5),  and  the  optimal  warning  rule 
(3.6). 

Relative  Operating  Characteristic 

Different  disutility  functions  may  result  in  different  threshold  values  q*  in  (3.6).  By  varying  the 
threshold  q*  throughout  its  range  (0,1),  one  can  generate  all  possible  warning  rules  that  could  result  from 
various  disutility  functions.  With  each  threshold  value,  there  is  associated  a  probability  of  detection  P(D) 
and  a  probability  of  false  warning  P(F).  A  plot  of  P(D)  versus  P(F),  obtained  by  varying  the  threshold  q*, 
is  called  the  relative  operating  characteristic  (ROC). 
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Figure  3-2.  Tree  of  Events  Leading  to  One  of  the  Four  Performance  States  of  a  Flood 

Warning  System 

(Missed  Flood  (M),  False  Warning  (F),  Detection  (D),  and  Quiet  (Q)) 
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The  ROC  curve  conveys  the  essential  information  about  the  tradeoffs  that  a  given  system  offers 
between  the  probability  of  detection  and  the  probability  of  false  warning.  However,  the  intuitive 
interpretation  of  these  performance  probabilities  is  not  straightforward  for  they  are  conditional  probabilities. 
In  fact,  P(D)  and  P(F)  are  the  probabilities  of  the  same  event  (w  =  1),  but  each  is  conditional  on  a  different 
vector  of  two  events,  as  shown  in  (3.7).  Human  intuition  does  not  grasp  easily  such  conditional  events. 
Moreover,  human  cognition  is  generally  not  well  trained  in  understanding  and  processing  probabilities. 
Evidence  of  numerous  and  large  biases  in  judgments  involving  probabilities  is  plentiful. 

Performance  Tradeoff  Characteristic 

In  order  to  overcome  the  interpretive  difficulties  associated  with  the  ROC,  we  propose  to  transform 
the  probabilities  of  various  states  into  the  expected  number  of  states  per  year.  Given  the  expected  number 
of  floods  per  year,  N,  the  following  quantities  can  readily  be  obtained: 

expected  number  of  zone  floods  per  year: 


n  =  N  •  P(0  =  1  I  0  =  1) 


(3.10) 


expected  number  of  detections  per  year  for  a  zone: 

ND  =  N  •  P(D)  •  P(e  =  1  1  0  =  1)  (3-11) 


expected  number  of  false  warnings  per  year  for  a  zone: 

NF  =  •  N  •  P(F)  •  [1  -  YP(e  =  1  I  0  =  1)]  (3-12) 

Through  equations  (3.11)  and  (3.12),  the  ROC  curve  can  be  rescaled  into  a  function  between  the 
expected  number  of  detections  and  the  expected  number  of  false  warnings  per  year.  This  function  will  be 
called  the  performance  tradeoff  characteristic  (PTC). 

Expected  Lead  Time 

The  forecast  time  is  the  instant  up  to  which  hydrometeorologic  observations  for  preparing  the 
forecast  are  collected.  The  lead  time,  ,  of  a  warning  for  a  given  zone,  conditional  on  the  hypothesis  that 
the  zone  will  be  flooded,  6  =  1,  is  the  time  interval  elapsed  from  the  forecast  time  to  the  instant  at  which 
the  flood  waters  reach  the  zone  elevation.  Let 

g(A  I  Y  =  1)  <3.13) 

denote  the  density  of  X,  conditional  on  event  0  =  1.  The  expected  lead  time  is  thus  given  by 


LT 


OO 

/ 


Xg{X  I  0  =  l)dX 


(3.14) 
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The  designer  of  a  warning  system  can  affect  the  lead  time  indirectly,  through  the  definitions  of  (i) 
the  forecast  trigger  and  (ii)  the  forecast  time.  Each  of  these  specifications  may  affect  all  three  probability 
densities  given  by  (3.4),  (3.5),  and  (3.13).  Consequently,  any  change  in  the  design  specifications  may 
simultaneously  affect  P(D),  P(F),  ND,  NF,  and  LT. 

Case  Study 

General  Description 

Properties  of  the  ROC  and  PTC  curves,  and  their  potential  role  as  aids  to  design  analysis,  will  be 
illustrated  through  a  case  study  of  the  flood  warning  system  for  Milton,  Pennsylvania.  The  town  has  a 
population  of  about  8000  and  is  located  on  the  West  Branch  of  the  Susquehanna  River  in  northeastern 
Pennsylvania.  The  data  used  in  the  smdy  were  collected  by  Krzysztofowicz  and  Davis  [1983].  The  source 
of  the  flood  and  forecast  data  is  the  U.S.  National  Weather  Service,  River  Forecast  Center  at  Harrisburg, 
Pennsylvania.  The  forecast  data  are  from  the  period  1959-1975.  Thus  the  case  studies  reported  herein  are 
representative  of  the  system  performance  during  that  period. 

In  all  specifications  of  the  parameters,  the  units  of  time  are  hours  and  the  units  of  elevation  are  feet 
above  the  zero  of  the  river  gauge.  The  flood  stage  is  at  19  ft,  but  almost  all  structures  are  located  above  the 
elevation  of  22  ft.  The  probability  densities  in  the  model  of  the  forecaster  are  assumed  to  follow  the 
Gaussian  law,  denoted  N(M,S^),  where  M  is  the  mean  and  is  the  variance. 

Input  Models  and  Parameters 

Record  offloads.  The  record  of  floods  from  the  period  1885-1975  contains  20  flood  events.  From 
this  record,  we  estimated: 

expected  number  of  floods  per  year:  N  =  0.53 

prior  density:  g(h  |  0  =  1)  =  N(|ih,  o^),  where  ^i,  =  24.9  and  =  4.82 

Models  of  Monitor  and  Forecaster.  The  parameters  that  must  be  estimated  from  the  joint  record  of 
forecasted  and  actual  floods  are  as  follows: 

diagnosticity  of  the  monitor:  y 

reliability  of  the  monitor:  p 

likelihood  functions:  Ko(s  |  0  =  0,  T  =  1)  =  N(/is,  a^) 

f(s  I  h,  0  =  1,  T  =  1)  =  N(ah  +  b,  o^) 

expected  lead  time  (for  each  zone  elevation  y): 

LT 
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The  likelihood  function  f  arises  from  the  linear  relation: 


s  =  ah  +  b  +  e 

where  e  is  a  random  variable,  stochastically  independent  of  h,  and  having  density  N(0,s ),  so  that 


E(s  I  h)  =  ah  +  b 
and 

Var(s  I  h)  =  ol 

Thus,  the  parameters  a,  b,  may  be  estimated  via  a  regression  analysis  performed  on  a  record  of  the 
forecasted  and  actual  flood  crests. 


Record  of  forecasts.  The  forecast  verification  reports  for  the  period  1959-1975  contain  a  record  of 
9  floods  and  37  forecasts.  The  record  does  not  contain  information  sufficient  for  the  estimation  of  a 
parameters  via  statistical  methods.  Consequently,  parameters  y  and  p  of  the  monitor  and  parmeters  of  Ae 
likelihood  function  Kq  had  to  be  estimated  subjectively  based  on  a  plausible  interpretation  and  interpolation 
of  the  available  information.  On  the  other  hand,  the  parameters  of  the  likelihood  function  f  and  the  expected 
lead  times  LT  were  estimated  statistically. 


System  designs.  The  monitor  is  assumed  to  trigger  the  forecaster  when  the  river  stage  ex^eds  a 
specified  threshold.  Three  alternative  system  designs  are  analyzed,  in  which  the  forecast  trigger  is  define 
as  follows: 


System  SI:  river  stage  >  11  ft 

System  S2:  river  stage  >  15  ft 

System  S3:  river  stage  >  19  ft 

The  likelihood  function  Kq  is  assumed  to  be  the  same  for  each  system;  the  estimates  of  its  parameters  are: 


=  17.8 


Os  =  117 

The  remaining  parameters  vary  with  the  design.  Table  3-1  lists  estimates  of  y,  P,  a,  b,  and  o.  Table  3-1 
lists  the  estimated  expected  lead  times,  LT,  and  the  calculated  expected  number  of  zone  floods,  n,  for  four 
zones  of  the  floodplain  extending  upwards  from  the  following  elevations: 


y  =  19,  22,  25,  28 

Since  the  expected  numbers  per  year  are  relatively  small,  we  rescale  them  in  the  discussion  to  a  100-year 
period. 
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Table  3-1.  Parameters  of  Three  Alternative  Designs  of  a 

for  Milton,  Pennsylvania 

Flood  Warning  System 

System 

Design 

Monitor 

Likelihood  Function 

Forecast 

Sufficiency 

Characteristic 

FSC 

Diagnosticity  Reliability 

Y  P 

Slope 

a 

Intercept 

b 

St.  Dev. 

0 

SI 

0.80 

1.00 

0.44 

10.65 

3.06 

6.95 

S2 

0.90 

0.89 

0.45 

12.10 

1.90 

4.22 

S3 

1.00 

0.83 

0.64 

8.48 

2.21 

3.45 

Table  3-2.  Expected  Number  of  Zone  Floods  and  Expected  Lead  Times  of  Flood  Warnings 

for  Milton,  Pennsylvania 

Zone 

Elevation 

yfft] 

Expecl^ed 
Number  of 
Zone  Floods 

n 

Expected  Lead  Time 

LT  [hrsj 

System 

SI 

System 

S2 

System 

S3 

19 

47.1 

9 

5 

-3 

22 

38.4 

15 

11 

4 

25 

26.0 

21 

17 

11 

28 

13.7 

27 

24 

18 

The  expected  numbers  are  for  the  period  of  100  years. 
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Interpretation.  Figure  3-3  shows  the  expected  lead  time  LT  plotted  as  a  function  of  die  elevation 
y  for  each  of  the  three  systems.  Clearly,  when  the  threshold  stage  for  triggering  the  forecaster  is  raised,  &e 
expected  lead  time  LT  is  reduced  uniformly  for  all  elevations.  Table  3-1  reveals  further  imphcations.  When 
the  lead  time  decreases,  the  diagnosticity  of  the  monitor  g  increases,  since  a  higher  threshold  stage  is  always 
more  diagnostic  of  the  incoming  flood.  On  the  other  hand,  when  the  lead  time  LT  decreases,  the  reliability 
r  also  decreases.  This  is  so  because  the  observations  of  the  river  are  made  in  6-hour  intervals,  and  it  is 
possible  for  a  rapidly  rising  river  to  exceed  both  the  threshold  stage  and  the  flood  stage  withm  the  6-hour 
interval.  In  such  an  instance,  flooding  occurs  prior  to  the  preparation  of  a  forecast.  The  likelihood  of  such 
an  event  increases  as  the  threshold  stage  is  raised  closer  to  the  flood  stage;  hence  the  rehabihty  r  decreases. 


When  the  expected  lead  time  LT  decreases,  one  also  anticipates  an  increase  in  the  quality  of  the 
flood  crest  forecasts.  Table  3-1  reveals  that  the  parameters  a,  b,  and  s  of  the  likelihood  function  change  their 
values  with  LT.  But  do  these  changes  imply  anything  about  the  forecast  quality?  The  answer  to  this  question 
may  be  obtained  via  the  forecast  sufficiency  characteristic: 


FSC  = 

This  measure  is  sufficient  for  comparing  any  two  forecasters  who  produce  forecasts,  say  Sm  s„,  of  the 
same  variate  h.  A  theorem  of  Krzysztofowicz  [1987]  states  that  FSC„  <  FSC„  if,  and  only  if ,  forecast  s^ 
has  the  economic  value  at  least  as  high  as  forecast  s„,  for  every  rational  decisionmaker.  In  other  words  the 
FSC  enables  us  to  order  forecasts  in  terms  of  their  economic  values.  The  FSCs  calculated  m  the  last  co  umn 
of  Table  3-1  confirm  our  hypothesis:  when  LT  decreases,  the  quality  of  the  flood  crest  forecasts  mcreases. 


Properties  of  the  ROC  and  PTC 

The  ROC  and  PTC  curves  for  design  SI  are  displayed  in  Figures  3-4  and  3-5.  We  shall  highlight 
some  general  properties  of  these  curves. 


1.  For  a  fixed  zone  elevation,  the  associated  ROC  is  a  concave  function  specifying  a  unique 
relationship  between  the  probability  of  false  warning,  P(F),  and  the  probability  of  detection,  P(D). 
Probability  P(F)  may  vary  from  0  to  1,  but  probability  P(D)  is  bounded  from  above  by  the  reliability  of  the 
monitor  r  which  for  design  SI  happened  to  be  1.0.  For  a  fixed  zone  elevation  and  a  prior  distribution  of 
the  flood  crest,  the  shape  of  the  ROC  curve  depends  solely  upon  the  design  specifications  for  the  nionitor 
(via  diagnosticity  g  and  reliability  r),  and  the  design  specifications  for  the  forecaster  (via  the  likelihood 

functions  Kq  and  f). 


2.  By  mapping  each  point  from  the  ROC  in  Figure  3-4  through  relations  (3.11)-(3.12),  we  obtain 
the  PTC  shown  in  Figure  3-5.  The  PTC  is  also  a  concave  function,  increasing  from  the  origin,  which 
corresponds  to  P(F)  =  P(D)  =  0,  to  a  point  which  corresponds  to  P(F)  =  1  and  P(p)  =  r. 
number  of  detections  per  year,  ND,  never  exceeds  the  expected  number  of  zone  floods  n.  On  Ae  oto 
hand,  the  expected  number  of  false  warnings  per  year,  NF,  may  exceed  n.  The  PTC  for  zone  elevations 
y  =  25  and  y  =  28  do  just  that. 


EXPECTED  LEAD  TIME  LT  [hrs] 
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ZONE  ELEVATION  y  [ft] 


Figure  3-3.  Expected  Lead  Time  of  a  Flood  Warning  Versus  the  Elevation  of 
the  Floodplain  for  Three  Warning  Systems  in  Milton,  Pennsylvania 


PROBABILITY  OF  DETECTION  P(D) 


PROBABILITY  OF  FALSE  WARNING  P(F) 


Figure  3-4.  Relative  Operating  Characteristics  (ROC)  of  Warning  System  SI 
for  Four  Zone  Elevations,  Y,  in  Milton,  Pennsylvania. 

(Symbol  LT  Denotes  the  Expected  Lead  Time  of  a  Flood  Warning) 


NUMBER  OF  DETECTIONS  ND 
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n 

47.1 


38.4 


26.0 


13.7 


Figure  3-5.  Performance  Tradeoff  Characteristics  (PTC)  of  Warning  System  SI  for 
Four  Zone  Elevations,  Y,  in  Milton,  Pennsylvania. 

(Symbol  n  Denotes  the  Expected  Number  of  Zone  Floods.  All  Numbers,  NF,  ND,  and  n  are  for 

the  Period  of  100  Years.) 


84 


Flood  Wamitii 


Risk-Based  Evaluation  of 
;  and  Preparedness  Systems 


Performance  Differences  Between  Zones 

1  The  ROC  curves  for  different  elevation  zones  have  generally  similar  shapes  and  cross  each  other. 

In  other  words,  when  the  performance  of  a  warning  system  is  measured  in  terms  of  the  probability  of  false 
warning  P(F)  and  the  probability  of  detection  P(D),  all  zones  seem  to  be  served  equally  well.  However,  the 
third  performance  measure  -  the  expected  lead  time  of  the  warning,  illuminates  the  differences  bemeen 
the  low-lying  zones  and  the  high-lying  zones:  LT  =  9  hours  for  y  =  19  and  LT  =  27  hours  for  y  -  28, 
a  threefold  difference. 

2.  The  PTC  curves  are  quite  dissimilar,  underscoring  the  fact  that  they  convey  different  information 
than  the  ROC  curves  do.  There  are  two  main  distinctions  between  the  zones.  First,  diere  is  the  obvious 
distinction  resulting  from  the  elevations  difference:  the  expected  number  of  floods  n  in  100  years  is  47.1 
for  y  =  19  and  only  13.7  for  y  =  28.  Second,  there  is  a  remarkable  difference  in  terms  of  the  expected 
number  of  false  warnings  NF  associated  with  the  maximum  expected  number  of  detections  ND  =  n.  This 
NF  is  equal  to  19.2  for  y  =  19,  and  it  increases  to  52.5  for  y  =  28.  In  other  words,  to  reach  the  upper 
limit  of  expected  detections  for  the  higher  zone,  one  must  accept  a  rate  of  false  warnings  NF  =  52.5,  which 
is  3.8  times  higher  than  the  rate  of  floodings  n  =  13.7. 

3 .  To  place  these  results  in  proper  perspective,  one  has  only  to  realize  that  floods  reaching  zone 
y  =  28  are  more  extreme  and  rare  than  floods  reaching  only  zone  y  =  19.  The  PTC  curves  in  Figure  3-5 
inform  us  that  a  high  detection  rate  for  rare  events  comes  at  the  price  of  a  high  rate  of  false  warnings.  This 
appears  to  be  an  inescapable  tradeoff. 

Operating  Points 

1.  A  point  on  the  ROC,  or  PTC,  is  called  an  operating  point.  In  Figure  3-4,  we  fixed  an  operating 
point  for  each  zone  such  that  for  all  zones  the  probability  of  false  warning  P(F)  =  Q.25.  The  probability 
of  detection  P(D)  is  different  for  each  zone,  but  the  differences  are  small.  Table  3-3  lists  the  exact 
coordinates  of  these  operating  points  on  the  PTC.  Figure  3-5  depicts  these  points,  each  of  which  has  distmct 
NF  and  ND  coordinates. 

2.  In  general,  the  mapping  between  the  operating  points  of  the  ROC  and  PTC  is  one-to-one,  with 
the  following  properties:  (i)  The  operating  points  which  have  the  same  P(D)  coordinate  on  the  ROC,  have 
also  the  same  ND  coordinate  on  the  PTC.  (ii)  The  operating  points  which  have  the  same  P(F)  coordinate 
on  the  ROC,  may  have  different  NF  coordinates  on  the  PTC.  We  shall  say  that  the  mapping  between  the 
ROC  and  PTC  is  nonorthogonal. 

3  The  nonorthogonality  of  the  mapping  between  the  ROC  and  PTC  should  be  taken  as  a  caution: 
judgmental  analysis  of  tradeoffs  on  the  ROC,  or  PTC,  is  not  a  simple  cognitive  task!  We  recommend  using 
the  PTC  as  the  primary  aid  to  planning  and  design  because  the  expected  number  ND  and  NF  are  easier  to 
interpret  and  understand  than  the  probabilities  P(D)  and  P(F),  which,  one  should  recall  here  expression 
(3.7),  are  conditional  probabilities. 

4.  With  each  operating  point  on  the  PTC,  or  ROC,  there  is  associated  a  unique  threshold  q  in  the 
warning  rule  (3.6).  Thus,  a  specification  of  the  operating  point  is  equivalent  to  a  specification  of  the  rule 
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Table  3-3.  Coordinates  of  Operating  Points  on  the  ROC  and  PTC  Curves  That  Give  the  Same 
Probability  of  False  Warning  P(F)  for  Each  Zone  Elevation;  System  Design  SI  for  Milton,  Pennsylvania 


Operating 

Point 

Zone 

Elevation 

yfft] 

Probability  of 

Expected  Number*  of 

Expected 
Lead  Time 
LT  [hrs] 

Detection 

P(D) 

False  W, 
P(F) 

Detections 

ND 

False  W. 
NF 

A 

19 

0.75 

0.25 

35.1 

4.8 

9 

B 

22 

0.71 

0.25 

27.3 

6.9 

15 

C 

25 

0.69 

0.25 

18.1 

10.0 

21 

D 

28 

0.72 

0.25 

9.9 

13.3 

27 

A 

The  expected  nuunbers 

are  for  the  period  of  100 

years . 

for  deciding  warnings.  To  specify  an  operating  point  on  the  PTC  for  a  given  zone,  one  should  consider  a 
tradeoff  between  the  expected  number  of  detections  ND  and  the  expected  number  of  false  warnings  NF. 
This  tradeoff  should  encapsulate  one's  preferences  for  outcomes  of  all  possible  decision-event  vectors  for 
this  particular  zone.  It  follows  that  it  would  be  irrational  to  fix  the  operating  point  based  solely  on  a 
displayed  PTC  or  ROC,  without  an  in-depth  analysis  of  all  socioeconomic  outcomes  of  every  decision-event 
vector.  That  is  why  the  PTC  and  ROC  curves  should  be  viewed  only  as  aids  to  the  planning  and  design 
process,  rather  than  as  a  means  of  specifying  the  warning  rule.  The  optimal  warning  rule  should  be  found 
by  minimizing  the  expected  disutility  of  outcomes  resulting  from  all  possible  decision-event  vectors. 

Performance  Tradeoffs 

From  a  purely  statistical  point  of  view,  which  ignores  the  economic  and  social  decision  criteria,  the 
engineer  could  consider  the  design  process  as  an  optimization  problem  with  three  criteria:  maximize  the 
expected  number  of  detections  ND,  minimize  the  expected  number  of  false  warnings  NF,  and  maximize 
the  expected  lead  time  LT.  The  ideal  solution  is  an  operating  point  having  the  coordinates  (ND,  NF,  LT) 
=  (n,  0,  oo).  In  the  absence  of  the  ideal  solution,  tradeoffs  must  be  made. 

The  kinds  of  tradeoffs  that  one  may  encounter  are  illustrated  for  the  three  alternative  system 
designs,  SI,  S2,  and  S3.  The  ROC  and  PTC  curves  of  these  systems  are  compared  in  Figures  3-6  and  3-7 
for  zone  elevation  y  =  22  and  in  Figures  3-8  and  3-9  for  zone  elevation  y  =  28.  An  immediate  observation 
is  that  designs  SI  and  S2  offer  distinct  performance  characteristics.  On  the  other  hand,  designs  S2  and  S3 
have  similar  ROC  and  PTC  curves  over  a  range  of  operating  points,  while  over  the  remaining  range  S2 
outperforms  S3.  Together  with  the  fact  that  S3  offers  much  shorter  expected  lead  times  LT  than  S2  does, 
it  is  unlikely  that  decisionmakers  would  prefer  S3  over  S2.  This  example  illustrates  then  a  screening  analysis 
that  may  be  performed  on  a  large  set  of  alternative  designs  before  a  few  are  selected  for  a  detailed  analysis 
of  tradeoffs. 
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The  ensuing  discussion  highlights  the  nature  of  performance  tradeoffs  that  the  PTC  allows  one  to 
analyze  between  designs  SI  and  S2.  The  discussion  concentrates  on  three  operating  points,  labeled  A,  B, 
C,  in  Figure  3-7.  Their  coordinates  (ND,  NF,  LT)  are  listed  in  Table  3-4. 

1 .  A  good  way  to  start  the  analysis  is  to  fix  the  expected  number  of  false  warnings  NF  at  a  level 
that  appears  acceptable,  at  least  initially,  say  NF  =  5.0,  which  means  that  one  would  be  willing  to  accept 
5.0  false  warnings  in  100  years,  on  the  average.  At  this  level  of  NF,  design  SI  ensures  the  expected 
detection  of  ND  =  23.8  floods  out  of  the  expected  n  =  38.4  floods  in  100  years.  The  expected  lead  tiine 
of  a  warning  for  each  detected  flood  is  LT  =  15  hours.  The  difference,  n  -  ND  =  38.4  -  23.8  —  14.6,  is 
the  expected  number  of  floods  in  100  years  that  will  arrive  undetected,  and  thus  will  not  be  preceded  by 
a  warning  to  the  public. 

2.  At  the  same  level  of  NF  =  5.0,  design  S2  ensures  the  expected  detection  of  28.9  floods  in  100 
years,  with  the  expected  lead  time  of  a  warning  equal  to  11  hours;  the  expected  number  of  missed  floods 
in  lo'o  years  is  38.4  -  28.9  =  9.5.  Thus,  when  comparing  the  operating  points  A  and  B,  the  following 
tradeoff  should  be  considered:  is  it  preferable  or  not  to  reduce  LT  from  15  to  11  hours  in  order  to  increase 
ND  from  23.8  to  28.9  (or,  equivalently  to  reduce  the  expected  number  of  missed  floods  from  14.6  to  9.5)? 

3.  A  similar  analysis  of  tradeoffs  may  be  carried  out  for  a  fixed  expected  number  of  detections  ND, 
say  28.9  in  100  years.  At  this  level  of  ND,  the  number  of  false  warnings  expected  in  100  years  is  5.0  for 
design  S2  and  8.1  for  design  SI;  the  accompanying  expected  lead  times  of  a  warning  are,  respectively,  11 
and  15  hours.  Thus,  when  comparing  the  operating  points  B  and  C,  the  following  tradeoff  should  be 
considered:  is  it  preferable  or  not  to  reduce  LT  from  15  to  11  hours  in  order  to  decrease  NF  from  8.1  to 
5.0? 


4.  The  right  endpoints  of  the  PTC  curves  indicate  that,  given  the  present  specifications  for  the 
monitor,  design  SI  can  detect  all  n  =  38.4  floods  expected  in  100  years.  However,  design  S2  has  an  upper 
limit  of  34.2  expected  detections  in  100  years;  thus,  the  minimum  expected  number  of  missed  floods  in  100 
years  is  38.4  -  34.2  =  4.2.  The  upper  limit  of  ND  is  achieved  by  each  design  at  a  different  level  of  the 
expected  number  of  false  warnings  NF,  which  is  27.8  for  design  SI,  and  18.2  for  design  S2. 

Closure 

The  relative  operating  characteristic  (ROC)  and  the  performance  tradeoff  characteristic  (PTC )  are 
a  part  of  a  general  theory  of  flood  warning  systems  that  is  being  developed.  A  number  of  questions  are  still 
awaiting  answers.  Among  them  is  the  connection  between  these  statistical  measures  of  performance  and  the 
ex  ante  economic  value  of  a  warning  system.  Such  a  connection  is  well  known  within  the  classical  detection 
paradigm,  but  it  remains  to  be  investigated  whether  or  not  it  extends  to  a  much  more  complex  paradigm 
of  a  flood  warning  system. 

Applications  of  ROC  and  PTC  concepts  to  other  flood  warning  systems  are  also  awaiting  us.  It 
would  be  desirable  to  make  a  number  of  applications  to  systems  with  distinct  hydrologic  regimes,  such  as 
flash-flood  streams  and  main-stem  rivers,  and  distinct  technologies,  such  as  found  in  local  warning  systems 
and  the  forecast  offices  of  the  National  Weather  Service.  Collectively,  results  of  such  case  studies  would 
offer  useful  guidance  to  engineers  who  plan  and  design  flood  warning  systems. 
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Figure  3-6.  Relative  Operating  Characteristics  (ROC)  of  Three  Warning  Systems,  SI,  S2,  and  S3 
for  Zone  Elevation  y  =  22  ft.  in  Milton,  Pennsylvania. 

(Symbol  p  Denotes  the  Reliability  of  the  Monitor.) 
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Figure  3-7.  Performance  Tradeoff  Characteristics  (PTC)  of  Three  Warning  Systems,  SI,  S2,  and 
S3,  for  Zone  Elevation  y  =  22  ft.  in  Milton,  Pennsylvania. 
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Figure  3-8.  Relative  Operating  Characteristics  (ROC)  of  Three  Warning  Systems, 
SI,  S2,  and  S3,  for  Zone  Elevation  y  =  28  ft  in  Milton,  Peimsylvania. 
(Symbol  p  Denotes  the  Reliability  of  the  Monitor.) 
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Figure  3-9.  Performance  Tradeoff  Characteristics  (PTC)  of  Three  Warning  Systems,  SI,  S2, 
and  S3,  for  Zone  Elevation  y  =  22  ft  in  Milton,  Pennsylvania 


Table  3-4.  Coordinates  of  Three  Alternative  Points  on  the  ROC  and  PTC  Curves 

for  Zone  Elevation  y  =  22  ft 


System 

Design 

Operating 

Point 

Probability  of 

Expected  Number*  of 

Expected 
Lead  Time 
LT  [hrs] 

Detection 

P(D) 

False  W. 
P(F) 

Detections 

ND 

False  W. 
NF 

SI 

A 

0,62 

0.18 

23.8 

5.0 

15 

SI 

B 

0.75 

0.29 

28.9 

8.1 

15 

S2 

C 

0.75 

0.27 

28.9 

5.0 

11 

*The  expected  numbers  are  for  the  period  of  100  years. 
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Part  4 

Selection  of  Optimal 
Flood  Warning  Threshold: 

Technical 

Introduction 

In  this  chapter  flood  warning  systems  are  studied  in  a  two-level  hierarchical  system  framework.  The 
interactions  between  the  forecast  subsystem  and  the  response  subsystem  are  investigated.  Emphasis  is  placed 
on  exploring  the  impact  of  the  current  selected  flood  warning  threshold  on  the  future  response  fraction  of 
a  flood  warning.  The  probabilistic  evaluation  of  a  forecast  system  coupled  with  a  stochastic  dynamic  model 
of  the  evolvement  of  the  response  fraction  in  a  community  reveals  that  the  desire  for  high  present  flood-loss 
reduction  must  be  balanced  with  the  possibility  of  high  future  flood  loss.  Multiobjective  dynamic 
programming  is  used  to  select  the  optimal  flood  warning  threshold.  The  proposed  methodology  is  applied 
to  the  case  study  in  Milton,  Pennsylvania. 


Features  of  the  Model 

Description  of  the  Methodology 

In  general,  the  overall  flood  warning  system  can  be  viewed  as  and  modeled  in  a  two-level 
hierarchical  system  framework  (see  Figure  4-1).  There  are  two  subsystems  at  the  lower  level.  One  is  the 
forecasting  subsystem,  which  issues  a  flood  forecast  based  on  hydrological  and  climatic  information.  The 
other  is  the  response  subsystem,  which  includes  decisionmaking  and  action  implementation  of  a  community 
in  response  to  flood  warning.  At  the  upper  level,  it  is  assumed  that  a  regional  agency  exists  whose  functions 
are  to  set  a  warning  threshold,  disseminate  a  flood  warning  to  the  community,  provide  transportation  during 
the  evacuation  process,  and  collect  statistical  data  of  the  warning  system. 

Performance  Measures  of  a  Warning  System 

Define  H  to  be  a  random  variable  which  represents  the  actual  flood  crest  and  S  to  be  a  random 
variable  which  represents  the  forecasted  flood  crest.  If  the  prior  probability  density  function  of  the  flood 
crest  is  denoted  by  g(h)  and  the  conditional  probability  density  function  of  s,  given  h,  is  denoted  by  f(s  [ 
h),  then  the  posterior  probability  density  function  of  h,  given  forecast  s,  is 

f(h  I  s)  =  f(s  I  h)  g(h)/k(s),  (4.1) 

where  k(s)  is  the  marginal  probability  density  function  of  forecast  s, 

k(s)  =  /  f(s  I  h)g(h)dh  (4.2) 

0 
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Figure  4-1.  Multilevel  Structure  of  Flood  Warning  System 


In  this  report  both  prior  distribution  and  the  likelihood  function  are  assumed  to  be  of  normal 


distributions  with  the  forms  [Krzysztofowicz  1987] 

g(h)  ~  N(Mh,ah^)  (4.3) 

and 

f(s  I  h)  ~  N(ah  +  b,a^).  (4.4) 

The  marginal  distribution  of  s  is  then  of  a  normal  distribution  [Krzysztofowicz  1987], 

k(s)  ~  N(a/ih  +  b,a^  +  (4.5) 

The  posterior  distribution  of  h  given  a  forecast,  s,  is  also  normal  [Krzysztofowicz  1987], 

f(h  I  s)  ~  N(As  +  B,C^),  (4.6) 

where 

A  =  ao\/(a^  +  a^oh^)  (4.7) 

B  =  (^*1,0^  -  abah^)/(o^  +  a^a,,^)  (4.8) 

=  0^0^ lip^  +  a^Oh^)-  (4-9) 
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A  missed  forecast  is  a  flood  event  which  is  not  preceded  by  a  warning.  The  probability  of  a  missed  flood 
warning  is 

Poi  (s*,y)  =  /P(h  ^  y  I  s)k(s)ds 
0 

=  P(h  >  y)  -  P,i  (s*,y) 

=  0.5{1  -  erf[y  -  -  P„  (s’,y) 

(4.14) 

A  correct  quiet  is  the  case  of  no  warning  and  no  flood  (h  <  y).  The  probability  of  acting  correctly 
in  not  issuing  a  flood  warning  is 

Poo  (s*,y)  =  (h  <  y  I  s)k(s)ds 

=  /fc(s)ds-Poi  (s*,y) 

0 

=  0.5(1  +  erfKs*  -  a/n^  -  b)/^2(o^  +  a^a^) ]}  -  Pq,  (s*,y)  (4.15) 

These  four  probabilistic  measures  are  related  to  each  other.  Knowing  one  of  them  and  the  prior  flood 
probability,  P(h  >  y),  and  the  probability  density  function  of  k(s),  the  other  three  can  be  calculated. 

There  are  two  types  of  prediction  errors  of  a  forecast  system  --  Type  I  and  Type  II  errors.  Type 
I  errors  are  those  of  missed  predictions.  Type  II  errors  are  those  of  false  alerts.  It  is  clear  from  Eqs.  (4.12)  - 
(4.15)  that  the  value  of  the  selected  threshold,  s*,  plays  a  key  role  in  determining  the  probabilities  of  T)q)e 
I  and  Type  II  errors.  If  the  threshold,  s*,  is  set  lower,  the  forecast  will  have  a  lower  probability  value  of 
a  Type  I  error  and  a  higher  probability  value  of  a  Type  II  error.  If  the  threshold,  s*,  is  set  higher,  the 
forecast  will  have  a  higher  probability  value  of  a  Type  I  error  and  a  lower  probability  value  of  a  Type  II 
error. 


Type  I  and  Type  II  errors  have  different  impacts  on  flood-loss  reduction.  A  Type  I  error  will  result 
in  an  immediate  flood  loss.  Thus,  it  has  a  short-term  impact.  On  the  other  hand,  a  Type  II  error  will  reduce 
the  credibility  of  the  forecast  system.  This  cry-wolf  consequence  has  a  long-term  impact.  Such  errors  do 
not  cause  a  flood  loss  at  the  present  stage,  but  will  discourage  the  response  to  flood  warnings  for  future 
flood  events.  The  present  fraction  of  people  who  respond  to  a  flood  warning  is  certainly  an  indispensable 
factor  in  constructing  the  flood-loss  function  for  a  community.  It  thus  affects  the  selection  of  the  flood 
warning  threshold.  On  the  other  hand,  the  response  fraction  fluctuates  as  time  passes,  based  on  the  past 
performance  of  the  warning  system.  The  coupling  between  successive  flood  events  is  carried  by  dynamic 
evolvement  of  the  fraction  of  people  who  respond  to  a  flood  warning.  The  past  performance  of  a  flood 
warning  system  affects  the  present  fraction  of  people  who  respond  to  the  warning  system.  Therefore,  it  is 
necessary  to  investigate  the  operation  of  flood  warning  systems  in  a  dynamic  framework.  Figure  4-2 
presents  the  interconnection  between  the  flood  warning  threshold  and  the  response  fraction  in  the  decision 
logic. 
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Figure  4-2.  Interaction  Between  Forecast  and  Response  Subsystems 


A  Model  of  the  Response  System 

In  general,  the  response  of  a  community  to  a  flood  warning  system  is  affected  by  people's 
experience  of  flooding  and  their  subjective  evaluation  of  the  past  performance  of  the  forecasting  system. 

The  general  interaction  between  a  forecasting  system  and  a  response  system  can  be  described  from 
the  following  considerations.  The  effectiveness  of  a  forecasting  system  can  be  judged  from  its  performance 
measures  of  Type  I  and  Type  n  errors.  The  response  of  a  community  to  a  flood  warning  can  be  described 
by  a  state  variable,  that  is,  the  fraction  of  people  in  the  community  who  respond  to  a  call  for  evacuation 
when  warned.  This  fraction  is  denoted  by  Oj  for  the  Tth  flood  event.  If  a  past  flood  event  has  been 
predicted,  then  confidence  in  the  flood  forecasting  system  will  increase,  and  thus,  future  rates  of  response 
will  also  increase.  On  the  other  hand,  a  cry-wolf  (Type  n  error)  event  will  decrease  confidence  in  the  flood 
forecasting  system,  thereby  decreasing  future  rates  of  response.  People  tend  to  have  decreased  confidence 
in  a  flood  warning  system  when  they  have  experienced  a  missed  warning.  However,  the  experience  of 
flooding  will  increase  people's  alertness  to  the  possibilities  of  future  floods.  For  simplicity,  it  is  reasonable 
to  assume  that  the  response  fraction  will  remain  unchanged  after  a  missed  warning  has  been  experienced. 
It  is  also  assumed  that  a  correct  quiet  does  not  change  the  response  fraction  in  the  future.  In  view  of  the 
above  discussion,  the  fraction  otj  can  be  assumed  to  evolve  dynamically,  being  governed  by  the  following 
equation: 


aT+1  = 


ccT  +  Cj(l  -  aJ) 
'  aT 
c^aT 


with  prob.  Pjj(5j,y) 

with  prob.  Poo(sj,y)  +  Poi(s^,y,) 

with  prob.  P^ffs^y) 


(4.16) 


where  c,  and  Cj  are  constants  or  functions  of  otj.  Their  range  is  in  (0, 1),  and  they  can  be  determined  using 
identification  methods  based  on  historic  data. 
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The  response  fraction  in  the  community  is  described  here  as  a  controlled  stochastic  process. 
Knowing  the  present  value  of  the  response  fraction,  three  possible  transitions  exist  with  given  probabilities. 
The  actual  transition  depends  on  the  real  outcome  associated  with  the  present  warning  decision.  This 
stochastic  system  can  be  controlled  in  the  sense  of  the  expected  value. 

Note  here  that  the  feedback  loop  that  encompasses  the  forecast  and  the  response  subsystems  is 
closed  only  when  the  next  flood  event  occurs.  The  present  performance  of  a  forecast  system  does  not  affect 
the  fraction  at  the  present  flood  event,  but  it  does  affect  the  fraction  ccj+i  at  the  next  flood  event. 

Multiobjective  Multistage  Optimization  Model 

A  key  aspect  of  flood  warning  systems  is  that  the  selection  of  the  flood  warning  threshold  cannot 
be  viewed  in  isolation  at  each  single  flood  event  since  the  decisionmaker  must  balance  the  desire  for  high 
present  flood-loss  reduction  with  the  possibility  of  high  future  flood  loss.  A  multiobjective  multistage 
optimization  model  has  been  proposed  by  Environmental  Systems  Modeling,  Inc.  [1990]  for  finding  the  best 
value  for  the  flood  warning  threshold  of  a  flood  warning  system.  Evaluating  the  trade-off  between  short- 
and  long-term  effects  yields  to  an  acceptable  balance  between  the  expected  loss  reduction  at  the  current 
stage  and  the  fraction  of  people  who  respond  at  the  next  flood  event. 

Assume  that  there  are  N  successive  flood  events  in  the  time  horizon  under  consideration.  At  each 
flood  event,  the  maximization  of  two  noncommensurate  and/or  conflicting  objective  functions  is  considered. 
At  the  Tth  flood  event,  one  objective  is  to  maximize  the  expected  property-loss  reduction;  the  second 
objective  is  to  increase  the  system's  credibility  by  reducing  the  cry-wolf  effect  (i.e.,  increasing  the  fraction 
of  people  who  respond  in  the  future). 

At  the  Tth  flood  event,  assume  that  the  expected  fraction  of  people  that  respond  to  the  warning 
system  is  The  decision  logic  for  property  loss  is  described  in  Figure  4-3  for  a  given  forecast  s,  where 
F  is  the  event  of  flooding,  w  is  the  action  of  flood  warning,  and  the  notation  denotes  the  negation  of 
an  event.  Denote  Dqo  (which  is  assumed  to  be  zero)  to  be  the  expected  loss  when  no  warning  is  given  and 
no  flood  occurs,  Doi(s,y)  to  be  the  expected  community  property  loss  without  a  warning,  D,o(aT>y)  to  be 
the  cost  of  the  evacuation  in  the  community,  and  Di,(s,aT.y)  to  be  the  expected  community  property  loss 
with  a  warning.  It  is  evident  that  the  response  fraction  affects  loss  functions  Dk,  and  Dj,.  For  a  specific 
value  of  forecast  s,  the  probability  that  a  zone  with  elevation  y  is  flooded  is  q(s,y)  =  P(h  >  y  |  s).  The 
expected  property  loss  with  no  warning  issued  is  q(s,y)Do,(s,y),  and  the  expected  property  loss  with  a 
warning  issued  is  q(s,y)D,i(s,0£T,y)  +  [1  -  q(s,y)]Dio(«T>y)-  Assume  that  the  flood  warning  is  issued  when 
the  forecast  flood  crest,  s,  exceeds  a  preassigned  threshold  level  s  j  .  The  expected  property  loss  with  a 

warning  system  then  becomes 

* 

Srj^ 

R{sj,aT,y)  =  /  {q(s,y)Doi(s,y)  +  [1  -  q(s,y)]Doo}k(s)ds 

00 

+ 1  {q(s,y)Dii(s,  ax.y)  +  [1  -  q(s,y)Dio(aT,y)}k(s)ds.  (4.17) 

Sij< 

Recall  that  the  expected  property  loss  without  a  warning  system  is 
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Ro(«T>y)  =  /q(s,y)Do,(s,y)k(s)ds.  (4.18) 

0 

Thus,  the  expected  property-loss  reduction,  / 1',  is 

00 

=  I  Q(s,y)(Doi(s,y)  -  q(s,y)D,i(s,ax,y)  -  [1  -  q(s,y)]D,o(aT.y)}k(s)ds  (4.19) 

* 

Sj’ 

To  construct  a  reasonable  loss  function,  the  concept  of  category-unit  loss  functions  proposed  by 
Krzysztofowicz  and  Davis  [1983]  is  adopted.  The  main  modification  is  that  the  notation  a,  which  was 
originally  used  in  Krzysztofowicz  and  Davis  [1983]  as  the  response  degree  of  an  individual,  is  used  in  this 
study  to  represent  the  fraction  of  people  who  respond  to  the  warning.  The  cost  function  of  evacuation  is 
assumed  to  be  a  linear  function  of  the  response  fraction 


d,o  =  MCa, 


(4.20) 


where  MC  is  the  maximum  evacuation  cost  for  the  community  when  a  full  response  is  present. 

The  flood-loss  function  without  a  warning  is  given  by 

doi  =  MD6(h  -  y),  (4.21) 

where  MD  is  the  maximum  possible  damage  due  to  flooding  of  the  highest  magnimde  when  no  response 
is  made,  and  6(h  -  y)  is  the  unit-damage  function  specifying  the  fraction  of  MD  which  occurs  when  the 
depth  of  flooding  is  (h  -  y). 


99 


Risk-Based  Evaluation  of 

Flood  Warning  and  Preparedness  Systems 


The  loss  fiinction  with  a  warning  is 

d„  =  MCa  +  MD[1  -  aMR(h  -  y)]6(h  -  y),  (4.22) 

where  MR(h  -  y)  is  the  unit  reduction  fiinction  specifying  the  reduction  of  MD  when  the  depth  of  flooding 
is  (h  -  y)  and  full  response  of  the  community  is  made,  a  =  1. 

The  expected  cost  of  the  community  evacuation  conditioned  on  a  given  forecast,  s,  is 


Djo  =  d,o  =  MCa. 


(4.23) 


The  expected  property  loss  without  a  warning  conditioned  on  a  given  forecast,  s,  is 

Doi=^  /  do,(h)f(h  I  s)dh. 

q(s,y)  y 

The  expected  property  loss  with  a  warning  conditioned  on  a  given  forecast,  s,  is 
1 


D„  = 


/  d„(a,h)f(h  I  s)dh. 
q(s,y)  y 


=  MCa  + 


Do,  -  [  MR(h  -  y)6(h  -  y)f(h  |  s)dh. 

q(s,y) 


In  this  case,  the  expected  property-loss  reduction  at  stage  T  can  be  simplified  into 

00  00 

=  af  [MD  I  MR(h  -  y)6(h  -  y)(h  |  s)dh  -  MC]  k(s)ds. 


(4.24) 


(4.25) 


(4.26) 


Note  here  that  the  expected  loss  reduction  is  a  linear  fiinction  of  the  response  fraction. 
When  MR(h  -  y)  and  6(h  -  y)  are  polynomial  functions  of  h,  the  integration  of 


/ 


MR(h  -  y)6(  h  -  y)f(h  |  s)dh. 


can  be  solved  using  the  following  formulas: 

f(h  I  s)dh  =  0.5[erf(^  ~  ~  -  erf{^  - 


v/2C 


V^c 


(4.27) 


/  hf(h  I  s)dh  =  [exp[-^^  -  As  -Bf 

o  2C^ 


]  -  exp[- 


(b  -  As  -  By 
2C^ 


]] 


-I-  (As  +  B)  f  f(h  I  s)dh 


(4.28) 


100 


Risk-Based  Evaluation  of 
Flood  Warning  and  Preparedness  Systems 


/  hJQi  I  s)dh  =  |[a"-^exp[-^^ — 

^  yJlL'K  y  2C7  2(7 

b  b 

+  (n  -  1)C^I  h''-^f{h  I  s)dh  +  (As  +  B)j  \  s)dh  (4.29) 

a  a 

m2 


The  N-stage  multiobjective  optimization  problem  of  flood  warning  systems  is  to  maximize  the  sum  of 
the  expected  property-loss  reductions  of  all  flood  events  over  the  time  horizon  under  consideration,  and  to 
maximize  the  forecast  system's  credibility,  which  is  implicitly  expressed  by  E(aN+i),  the  expected  fraction 
of  people  who  respond  to  the  warning  beyond  the  time  horizon  under  consideration.  Mathematically,  this 
overall  multiobjective  optimization  problem  can  be  posed  as: 


max 


(4.30a) 


S.t.  ^X+1 


ttj.  +  Cj(l  -  af)  with  prob.  P^fsp,y) 


a. 


with  prob.  Poo(sj,y)  +  Pofsp^y) 
with  prob.  P^Q(sp,y) 


(4.30b) 


where  is  given  in  Eq.  (4.26).  It  may  be  useful  to  consider  also  the  maximization  of  a  third  objective 
function  in  Eq.  (4.30a),  namely,  the  expected  life-loss  reductions  of  all  flood  events  over  the  time  horizon 
under  consideration.  The  life-loss  reductions  will  not  be  discussed  further  in  this  report. 


Solving  the  multiobjective  multistage  optimization  problem  in  Eq.  (4.30)  yields  the  set  of  noninferior 
solutions.  A  decision  sequence  s*  =  [suS2>---^Nli  is  said  to  be  noninferior  if  there  does  not  exist  another 
decision  sequence  s*  =  [si,S2,---Sji\f]  such  that  f(s*)  ^f(s*),  i  =  1,2,  with  at  least  one  strict 
inequality. 

Equation  (4.30)  can  be  solved  by  the  weighting  method  and  dynamic  programming.  At  the  final 
stage,  the  following  optimization  problem  is  solved  for  a  given  weighting  coefficient,  0,  and  a  given  value 
of  (Xn: 

V0,a^)  =  ”1"  0  [fY  I  “at]  +  (1  -  0)  \  a4  (4.31) 

subject  to  Eq.  (4.30b). 
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At  stage  T  (1  ^  T  ^  N  -  1),  the  following  optimization  problem  is  solved  for  a  given  value  of  the 
weighting  coefficient,  0,  and  a  given  value  of  cct^ 


4)/e,  a,)  =  0  l/f  1  OL,]  +  £{(j)^,i(0,  I  a^} 


(4.32) 


where  %  satisfies  Eq.  (4.30b). 

Case  Studies 


Application  to  Milton,  Pennsylvania 

System  design  S2,  described  in  this  report,  Part  3,  Performance  Characteristics  of  a  Flood  Warning 
System,  developed  in  a  study  of  Milton,  Pennsylvania,  is  selected  as  the  illustrative  application  of  the 
methodology  described  in  the  above  sections.  In  this  case,  the  flood  crest  H  is  of  a  normal  distribution 
N(24. 88,4.8216387^)  and  the  conditional  probability  density  function  of  the  forecasted  flood  crest  S,  given 
h,  is  of  a  normal  distribution  N(0.4503h  +  12.1044,1.8973^).  It  can  be  shown  that  (1)  the  marginal 
probability  density  function  of  the  forecast  k(s)  is  N(23.307864,2. 8833638^)  and  (2)  the  posterior 
distribution  density  function  of  h  given  s,  f(h  |  s),  is  N(l. 2591933s  -  4.4691054,3.1727162^).  The 
probability  of  flooding,  given  forecast  s,  is  obtained  by 
00 

q(s,y)  +/  f(s  I  h)dh 

y 

=  1  -  «I>[(y  -  1.2591933s  +  4.4691054)/3. 1727162],  (4.33) 

where  $(•)  is  the  cumulative  distribution  function  of  the  standard  normal.  In  our  computer  program,  the 
cumulative  distribution  function  0(«)  is  calculated  through  its  relationship  with  the  error  function, 

4>(x)  =  0.5[1  +  erf(x/\/2)],  (4.34) 


where  the  error  function  is  defined  by 


erf(x)  =  -^  /«■“' 
\/n  0 


du 


(4.35) 


Substituting  k;(s)  and  q(s,y)  into  Eqs.  (4.12)  -  (4.15),  the  four  probabilistic  measures  of  a  forecasting 
system  can  be  calculated  for  given  warning  threshold  s*  and  zone  elevation  y.  Table  4-1  shows  those 
measures  for  two  values  of  y  and  various  values  of  s*.  A  tradeoff  between  Type  I  and  II  errors  can  be 
clearly  seen  from  Fig.  4-4.  Different  values  of  s*  associate  different  values  of  the  probabilistic  measures, 
P,i(s*),  Pio(s*),  Poi(s*),  and  Poo(s*).  They  thus  yield  different  impacts  on  the  response  fraction  at  the 
subsequent  stage. 


For  simplicity,  we  only  consider  property  losses  in  this  illustrative  example.  To  construct  a  reasonable 
loss  function,  the  concept  of  category-unit  loss  functions  proposed  by  Krzysztofowicz  and  Davis  [1983d] 
is  adopted  in  the  following.  The  main  modification  is  that  the  notation  a,  which  was  originally  used  in 
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Krzysztofowicz  and  Davis  [1983d]  as  the  response  degree  of  an  individual,  is  used  here  to  represent  the 
fraction  of  people  who  respond  to  the  warning.  The  cost  function  is  assumed  to  be  of  the  form 

d,o  =  MCa,  (4.36) 


Table  4-1.  Probabilistic  Measures  of  the  Warning  System 
y  -  19 


'k 

'k 

k 

k 

S 

Pio(s  .y) 

Poi(s  ,y) 

Pqo^"  .y) 

15.0 

0.8885870 

0.1094329 

0.0000862 

0.0018939 

15.5 

0.8884591 

0.1081553 

0.0002140 

0.0031715 

16.0 

0.8881721 

0.1061974 

0.0005011 

0.0051294 

16.5 

0.8875671 

0.1033221 

0.0011061 

0,0080047 

17.0 

0.8863676 

0.0992858 

0.0023056 

0.0120411 

17.5 

0.8841263 

0.0938829 

0.0045469 

0.0174439 

IS.O 

0.8801783 

0.0870001 

0.0084949 

0.0243263 

18.5 

0.8736137 

0.0786737 

0.0150595 

0,0326531 

19.0 

0.8632919 

0.0691255 

0.0253813 

0.0422013 

19.5 

0.8479231 

0.0587649 

0.0407501 

0.0525619 

y  -  22 


k 

S 

Pll(s*,y) 

Pio(s*.y) 

Poi(s*,y) 

19.0 

0.7171233 

0.2152941 

0.0077255 

0.0598570 

19.5 

0.7108782 

0.1958098 

0.0139707 

0.0793414 

20.0 

0.7009019 

0.1734535 

0.0239469 

0.1016977 

20.5 

0.6858719 

0.1490539 

0.0389770 

0.1260972 

21.0 

0.6644858 

0.1237766 

0.0603631 

0.1513745 

21.5 

0.6357059 

0.0989638 

0.0891429 

0.1761873 

22.0 

0.5990129 

0.0759250 

0.1258360 

0.1992262 

22.5 

0.5546101 

0.0557205 

0.1702388 

0.2194306 

23.0 

0.5035065 

0.0390087 

0.2213424 

0.2361424 

23.5 

0.4474492 

0.0259864 

0.2773997 

0.2491647 
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Figure  4-4.  Trade-off  Between  Type  I  and  Type  II  Errors 
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where  MC  =  4  ($  ten  thousand)  is  the  evacuation  cost  for  the  community. 


The  flood-loss  function  without  a  warning  is  given  by 

do,  =  MD6(h  -  y),  <4.37) 

where  MD  =  100  ($  ten  thousand)  is  the  maximum  possible  damage  due  to  flood  of  the  highest  magnitude 
when  no  response  is  made,  and  6(h  -  y)  is  the  unit  damage  function  specifying  the  fraction  of  MD  which 
occurs  when  the  depth  of  flooding  is  (h  -  y).  We  assume  6(h  -  y)  to  be  equal  to  0.4(h  -  y)  when  h  -  y  ^  2.5 
and  to  be  equal  to  1  when  h  -  y  >  2.5  (h  and  y  are  in  the  unit  of  10  feet). 


The  loss  function  with  a  warning  is 

d,,  =  MCa  -I-  MD[1  -  aMR(h  -  y)]6(h  -  y),  (4-38) 

where  MR(h  -  y)  is  the  unit  reduction  function  specifying  the  reduction  of  MD  when  the  depth  of  flooding 
is  (h  -  y)  and  full  response  of  community  is  made,  a  =  1.  We  assume  MR(h  -  y)  to  be  equal  to  0.25  -I- 
0.4(h  -  y)  -  0.333(h  -  yf  when  h  -  y  <  1.2  and  to  be  equal  to  0.25  when  h  -  y  >  1.2. 

The  expected  cost  of  the  community  evacuation  conditioned  on  a  given  forecast  s  is 

Dio  =  dio  =  MCa.  <4.39) 

The  expected  property  loss  without  a  warning  conditioned  on  a  given  forecast  s  is 

Do,  =  -7^/  doi(h)f(h  I  s)dh 

qis,y)  y 


nAT^  y+2.5 

=  — —  /  QAQi  -  fjfQi  1  s)dh  +  !  f(h  \  s)dh  . 

q{s,y)  [  y  y-2.5  J 

The  expected  property-loss  with  a  warning  conditioned  on  a  given  forecast  s  is 


(4.40) 


D,,  =  d„(a,h)f(h  I  s)dh 

qis,y)  y 


aMD 


q(s,y) 


y^2.5 


>+1.2 


/  [0.25  +  OAQi  -  y)  -  0.333(h  -  yf]0.4(h  -  y)fih  1  s)dh 


+  /  0.25-  OAQi  -  yW  \  s)dh  +  /  0.25  f(h  \  s)dh 

y*1.2  >’■►2.5 


(4.41) 
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Since  f(h  |  s)  is  of  normal  distribution  N(l. 2591933s  -  4.4691054,  3.1727162^), 
integrations  can  be  performed  using  the  following  equations: 

b  -  1.2591933J  +  4.4691054 
3.1727162  ^ 

^.a  -  1.2591933^  +  4.4691054. 

^  3.1727162  ^ 


fih  I  s)dh  =  a>( 


L,  I  3.1727162 

jf(h  I  s)dh  =  - — — 

a  fin 


(a  -  1.25919335  +  4.4691054)^ 
2  •  3.17271622 


-  Exp\ 


(b  -  1.2591933.?  +  4:4691054)2 
2  •  3.17271622 


b 

+  (1.2591933s  +  4.4691054)  /  f(h  |  s)dh 

a 


Ih^f^h  I  s)dh  - 
a  ^Jln 


(a  -  1.25919335  +  4.4691054)2 
2  •  3.17271622 


.  tFxpl  ~  ^-^^919335  +  4.4691054)2 
2  •  3.17271622  J 

b  b 

+  3.17271622/  f(h  I  s)dh  +  (1.2591933s  +  4.4691054)/  hf(h  |s)dh 

a  a 


b 

jh^fih  I  s)dh 

a 


3.1727162 


n^Exp\ 


(a 


-  1.25919335  +  +4.4691054)2 
2  •  3.17271622 


b  2£Tjr[  ~  1 -2591933.?  +  4.4691054)2 
2  •  3.17271622 


b  b 

+  2  •  3. 17271622 /hf(h  |  s)dh  +  (1.2591933s  +  4.4691054)  /h2f(h  |  s)dh 


a 


a 


the  above 


(4.42) 


(4.43) 


(4.44) 


(4.45) 
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The  objective  function  /|^can  now  be  posed  as 

00 

/f  =  I  {q(s,y)Doi(s.y)  -  q(s,y)Dii(s,aT,y)  -  (1  -  q(s,y))Dio(aT.y)}k(s)ds  (4.46) 

ST 

Table  4-2  gives  the  calculated  values  of  the  expected  flood-loss  reduction  for  various  response 
fractions  and  preselected  warning  thresholds  where  y  is  equal  to  19  or  22.  The  relationship  between  the 
expected  flood-loss  reduction  and  the  warning  threshold  is  also  depicted  in  Fig.  4-5  for  various  response 
fractions  and  elevation  levels.  In  this  illustrative  example,  we  assume  that  the  parameters  in  the  response 
fraction  dynamic  model,  Eq.  (4.16),  take  values  of  Ci  =  0.1  and  g  =  0.9.  We  consider  a  five-stage 
problem  with  the  initial  responding  fraction  a,  equal  to  0.7.  Two  values  of  the  elevation  y  are  considered, 
y  =  19  (feet)  and  y  =  22  (feet). 


The  overall  problem  can  now  be  posed  as  follows: 


[f2  =  E  {%} 


(4.47a) 


s.t. 


Kj,  +  0.1(1 

®T+i  ~ 

0.9a7. 


-  aT)  with  prob.  Pjj(j'j’,y) 

with  prob.  Pooisj,y)  +  Poiispj) 
with  prob.  P^Q(s^,y) 


tti  =  0.7 


T  =  1,2,3,4,5 


(4.47b) 


Equation  (4.47)  can  be  solved  by  the  weighting  method  and  dynamic  programming.  At  the  fifth  stage,  the 
following  optimization  problem  is  solved  for  given  weighting  coefficient  0  and  given  value  of  Us: 

®5(0.«5)  =  I  %]  +  (1  -  0)  E  {a,  I  ttj},  (4.48) 

■^5 


subject  to  Eq.  (4.47b)  with  T  =  5. 

At  stage  T,  T  =  4,3,2, 1,  the  following  optimization  problem  is  solved  for  given  value  of  the 
weighting  coefficient  0  and  given  value  of 

e\fT  1  aJ  -t-  E{$T-hl  (0,  a^+l)  |  ccr}, 
where  aj  satisfies  Eq.  (4.47b),  and  =  0.7. 
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5(b) 

Figure  4-5.  Relationship  Between  the  Expected  Flood-loss  Reduction  and 
the  Warning  Threshold  for  Various  Response  Fraction  and  Elevation  Levels 
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5(d) 

Figure  4-5.  (continued) 
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Table  4-2.  Expected  Flood-loss  Reduction 
y  -  19 


a 

s 

0.50 

0.60 

0.70 

0.80 

16.5 

6.013321 

7.215985 

8.418648 

9.621314 

17.0 

6.014395 

9.623033 

17.5 

6.011573 

7.213888 

8.416201 

9.618517 

18.0 

5.998908 

7.198690 

8.398471 

9.598253 

18.4 

5.975090 

7.170107 

8.365125 

9.560144 

y  -  22 


a 

s 

0.50 

0.60 

0.70 

0.80 

19.0 

3.592137 

4.310565 

5.028992 

5.747422 

19.5 

3.595634 

4.314761 

5.033888 

5.753016 

20.0 

3.574387 

4.289264 

5.004142 

5.719021 

20.5 

3.511868 

4.214242 

4.916617 

5.618992 

20.9 

3.420732 

4.104879 

4.789026 

5.473174 

To  make  the  computational  procedure  feasible,  the  state  space  of  oct  and  the  control  space  of  s^  are 
quantized  by  the  grid  sizes  0.01  and  0.1,  respectively.  Table  4-3  (a-h)  provides  noninferior  solutions  at 
various  stages  and  various  ax's  for  four  different  weighting  coefficients  and  two  different  values  of  zone 
elevation.  It  can  be  seen  from  Table  4-3  that: 

(1)  the  lower  the  weighting  coefficient  0  associated  with  the  first  objective  (loss  reduction),  the 
higher  the  value  of  the  flood  warning  threshold  will  be  set  to  avoid  possible  high  Type  II  errors; 

(2)  in  order  to  select  a  decision  that  maximizes  the  sum  of  flood-loss  reductions,  the  flood  warning 
threshold  is  set  higher  at  the  earlier  stage  than  at  the  later  stage  (with  respect  to  the  same  value  of  the 
response  fraction)  in  order  to  reduce  the  probability  of  high  loss  at  the  later  stages;  and 

(3)  the  higher  the  present  response  fraction,  the  more  cautious  the  selection  of  the  threshold  is.  That 
means  that  a  higher  value  of  threshold  is  set  for  a  higher  value  of  the  present  response  fraction  in  order  to 
avoid  losing  a  larger  number  of  the  response  population.  We  should  note  here  that  the  third  conclusion  may 
be  model-specific. 

Figures  4-6,  4-7,  and  4-8  show  the  optimal  flood  warning  threshold  as  functions  of  the  weighting 
coefficient,  the  stage,  and  the  response  fraction,  respectively. 
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Q 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

17.6 

1.084007 

17.6 

0.9148158 

17.5 

0.7508677 

17.5 

0.5927224 

0.50 

17.7 

1.158004 

17.6 

0.9834811 

17.6 

0.8136330 

17.6 

0.6489545 

0.55 

17.8 

1.232019 

17.7 

1.052161 

17.7 

0.8764063 

17.6 

0.7051916 

0.60 

17.8 

1.306050 

17.8 

1.120851 

17.7 

0.9391865 

17.7 

0.7614322 

0.65 

17.9 

1.380095 

17.8 

1.189550 

17.8 

1.001974 

17.7 

0.8176754 

0.70 

! 

18.0 

1.6522150 

17.9 

1.454151 

17.9 

1.258261 

17.8 

1.064767 

17.8 

0.8739213 

0.75 

18.0 

1.528219 

17.9 

1.326976 

17.9 

1.127564 

17.8 

0.9301701 

0.80 

18.0 

1.602293 

18.0 

1.395698 

17.9 

1.190364 

17.8 

0.9864190 

)The  upper  number  denotes  the  selected  flood  warning  threshold 
and  the  lower  number  denotes  the  expected  loss  reduction  rom 
the  present  stage  to  the  final  stage. 
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Table  4-3b.  Optimal  Flood  Warning  Thresholds 


Y  -  19  9  -  0.06 


a 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

17.5 

2.045543 

17.4 

1.603284 

17.3 

1.183760 

17.2 

0.7894264 

0.50 

17.5 

2.201129 

17.5 

1.735747 

17.4 

1.290610 

17.2 

0.8678907 

0.55 

17.6 

2.356729 

17.5 

1.868221 

17.4 

1.397466 

17.3 

0.9463577 

0.60 

17.6 

2.512346 

17.6 

2.000699 

17.4 

1.504322 

17.3 

1.024826 

0.65 

17.7 

2.667974 

17.6 

2.133191 

17.5 

1.611183 

17.3 

1.103294 

0.70 

17.8 

3.3908580 

17.7 

2.823610 

17.6 

2.265684 

17.5 

1.718046 

17.3 

1.181762 

0.75 

17.8 

2.979251 

17.6 

2.398177 

17.5 

1.824910 

17.3 

1.260230 

0.80 

17.8 

3.134909 

17.7 

2.530680 

17.5 

1.931774 

17.4 

1.338699 

*)The  upper  number  denotes  the  selected  flood  warning  threshold; 
and  the  lower  number  denotes  the  expected  loss  reduction  from 
the  present  stage  to  the  final  stage. 
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Table  4-3c.  Optimal  Flood  Warning  Thresholds 


Y  =  19  ^  -  0.1 


Q 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.A5 

17.4 

3.007229 

17.4 

2.291853 

17.3 

1.616717 

17.1 

0.9861612 

0.50 

17.5 

3.244441 

17.4 

2.488155 

17.3 

1.767683 

17.1 

1.086875 

0.55 

17.5 

3.481679 

17.5 

2.684460 

17.3 

1.918650 

17.2 

1.187590 

0.60 

17.6 

3.718935 

17.5 

2.880785 

17.4 

2.069619 

17.2 

1.288307 

0.65 

17.6 

3.956197 

17.5 

3.077110 

17.4 

2.220595 

17.2 

1.389023 

0.70 

17.8 

5.1299780 

1 

17.7 

4.193476 

17.6 

3.273436 

17.4 

2.371570 

17.2 

1.489740 

0.75 

17.7 

4.430770 

17.6 

3.469775 

17.4 

2.522546 

17.2 

1.590457 

0.80 

17.7 

4.668063 

17.6 

3.666113 

17.4 

2.673521 

17.2 

1.691174 

)The  upper  number  denoree  the  selected  flood  warning  threshold 
and  the  lower  number  denotes  the  expected  loss  reduction  from 
the  present  stage  to  the  final  stage. 
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Table  4-3d.  Optimal  Flood  Warning  Thresholds 


Y  -  19  -  0.5 


a 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

lIBi 

liH 

17.0 

2.953618 

0.50 

17.4 

13.67814 

17.3 

10.01264 

17.2 

6.538699 

17.0 

3.276874 

0.55 

17.5 

14.73185 

17.4 

10.84742 

17.2 

7.130880 

17.0 

3.600129 

0.60 

17.5 

15.78559 

17.4 

11.68223 

17.2 

7.723061 

17.0 

3.923386 

0.65 

17.6 

16.83937 

17.4 

12.51704 

17.3 

8.315259 

! 

17.0 

4.246640 

0.70 

17.7 

22.522440 

17.6 

17.89324 

17.5 

13.35188 

17.3 

8.907464 

17.0 

4.569896 

0.75 

17.6 

18.94711 

17.5 

14.18674 

17.3 

9.499667 

17.0 

4.893150 

0.80 

17.6 

20.00097 

17.5 

15.02160  ; 

17.3 

LO. 09187 

17.0 

5.216409 

*)The  upper  number  denotes  the  selected  flood  warning  threshold; 
and  the  lower  number  denotes  the  expected  loss  reduction  from 
the  present  stage  to  the  final  stage. 
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Table  4-3e.  Optimal  Flood  Warning  Thresholds 


Y  -  22  8  -  0.02 


a 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

19.8 

0.8289434 

19.8 

0.7296401 

19.8 

0.6317911 

19.8 

0.5355290 

0.50 

19.9 

0.8879178 

19.9 

0.7864377 

19.9 

0.6862146 

19.9 

0.5873594 

0.55 

20.0 

0.9469303 

20.0 

0.8432720 

20.0 

0.7406604 

19.9 

0.6392034 

0.60 

20.1 

1.005980 

20.0 

0.9001282 

20.0 

0.7951348 

20.0 

0.6910616 

0.65 

20.1 

1.065078 

20.1 

0.9570246 

20.1 

0.8496162 

20.0 

0.7429260 

0.70 

20.2 

1.2348560 

20.1 

1.124184 

20.1 

1.013937 

20.1 

0.9041280 

20.1 

0.7947990 

0.75 

20.2 

1.183335 

20.2 

1.070861 

20.1 

0.9586413 

20.1 

0.8466820 

0.80 

20.2 

1.242499 

20.2 

1.127814 

20.2 

1.013166 

20.1 

0.8985649 

*)The  upper  number  denotes  the  selected  flood  warning  threshold, 
and  the  lower  number  denotes  the  expected  loss  reduction  from 
the  present  stage  to  the  final  stage. 


Risk-Based  Evaluation  of 

Flood  Warning  and  Preparedness  Systems _ 

Table  4-3f.  Optimal  Flood  Warning  Thresholds 
Y  -  22  e  ~  0.06 


a 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

19.7 

1.373070 

19.6 

1.121546 

19.6 

0.8787987 

19.5 

0.6456339 

0.50 

19.7 

1.480653 

19.7 

1.216034 

19.6 

0.9588957 

19.6 

0.7099622 

0.55 

19.8 

1.588269 

19.7 

1.310529 

19.7 

1.039008 

19.6 

0.7742918 

0.60 

19.8 

1.695913 

19.8 

1.405033 

19.7 

1.119128 

19.6 

0.8386213 

0.65 

19.8 

1.803563 

19.8 

1.499559 

19.7 

1.199249 

19.6 

0.9029509 

0.70 

19.9 

2.2306700 

19.9 

1.911242 

19.8 

1.594090 

19.7 

1.279371 

19.6 

0.9672803 

0.75 

19.9 

2.018938 

19.8 

1.688627 

19.7 

1.359500 

19.7 

1.031617 

0.80 

19.9 

2.126638 

19.8 

1.783173 

19.8 

1.439644 

19.7 

1.095958 

*)The  upper  number  denotes  the  selected  flood  warning  threshold; 
and  the  lower  number  denotes  the  expected  loss  reduction  from 
the  present  stage  to  the  final  stage. 
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Table  4-3g.  Optimal  Flood  Warning  Thresholds 


Y  -  22  5-0.1 


Q 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

19.6 

1.917686 

19.6 

1.513842 

19.5 

1.126039 

19.5 

0.7558652 

0.50 

19.7 

2.074062 

19.6 

1.646128 

19.6 

1.231924 

19.5 

0.8327414 

0.55 

19.7 

2.230452 

19.7 

1.778437 

19.6 

1.337814 

19.5 

0.9096177 

0.60 

19.8 

2.386854 

19.7 

1.910764 

19.6 

1.443704 

19.5 

0.9864939 

0.65 

19.8 

2.543298 

19.7 

2.043091 

19.6  ^ 
1.549595 

19.5 

1.063370 

0.70 

19.9 

3.2281410 

19.8 

2.699744 

19.7 

2.175418 

19.6 

1.655485 

19.5 

1.140247 

0.75 

19.8 

2.856199 

19.7 

2.307755 

19.7 

1.761377 

19.5 

1.217123 

0.80 

19.8 

3.012672 

19.8 

2.440112 

19.7 

1.867286 

19.5 

1.293999 

*)The  upper  number  denotes  the  selected  flood  warning  threshold; 
and  the  lower  number  denotes  the  expected  loss  reduction  from 
the  present  stage  to  the  final  stage. 
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Table  4-3h.  Optimal  Flood  Warning  Thresholds 


Y  =  22  6  -0.5 


a 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

19.6 

7.365396 

19.5 

5.437704 

19.5 

3.599206 

19.4 

1.858485 

0.50 

19.6 

8.009911 

19.6 

5.948506 

19.5 

3.963165 

19.4 

2.061026 

0.55 

19.7 

8.654481 

19.6 

6.459373 

19.5 

4.327123 

19.4 

2.263566 

0.60 

19.7 

9.299195 

19.6 

6.970240 

19.5 

4.691081 

19.4 

2.466107 

0.65 

19.7 

9.943911 

19.6 

7.481107 

19.5 

5.055040 

19.4 

2.668648 

0.70 

19.8 

13.207410 

19.7 

10.58863 

19.6 

7.991975 

19.5 

5.418999 

19.4 

2.871188 

0.75 

19.7 

11.23334 

19.6 

8.502842 

19.5 

5.782957 

19.4 

3.073730 

0.80 

19.7 

11.87806 

19.6 

9.013708 

19.5 

6.146915 

19.4 

3.276271 

*)The  upper  number  denotes  the  selected  flood  warning  threshold; 
and  the  lower  number  denotes  the  expected  loss  reduction  from 
the  present  stage  to  the  final  stage. 
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stage  4 
a  -0.55 


Figure  4-6.  Relationship  Between  the  Optimal  Warning  Threshold 
and  the  Weighting  Coefficient  for  a  =  0.55  at  Stage  4 
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Figure  4-8.  Relationship  Between  the  Optimal  Warning  Threshold  and  the 
Response  Fraction  for  0  =  0.02  at  Stage  3 
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Input  Information  for  Eldred  and  Connellsville 

Flood  Forecasts  for  Eldred 


Location 

Eldred  is  a  small  community  situated  on  the  upper  Allegheny  River  in  northern  Pennsylvania,  about 
three  miles  south  of  the  state  border  with  New  York.  The  river  gauge  has  a  datum  at  1417  feet  and  closes 
a  drainage  area  of  50  square  miles  of  mountainous  terrain  with  highest  ranges  towering  at  2500  feet.  The 
river  flow  at  Eldred  is  essentially  unimpaired  natural  runoff. 

Data  Records  and  Parameter  Estimation 

Historical  flood  and  forecast  data  were  retrieved  from  the  archives  of  the  National  Weather  Service 
Forecast  Office  in  Pittsburgh.  The  prior  distribution  of  flood  crests  was  estimated  from  a  record  of  floods 
spanning  1942-1989.  During  these  48  years,  36  flood  crests  exceeded  the  gauge  height  of  1 1  feet  (3  floods 
in  every  4  years,  on  the  average),  and  14  had  crests  above  the  official  stage  of  17  feet  (about  one  flood  in 
every  3.5  years,  on  the  average).  The  highest  flood  on  record  occurred  in  June  1972  and  reached  29  feet. 

The  likelihood  functions  were  estimated  from  a  historical  joint  record  of  forecasted  and  actual  flood 
crests.  This  record  contained  12  floods  that  occurred  in  the  period  1984-1988. 

Flood  Forecasts  for  Connellsville ' 


Location 

Connellsville,  a  town  in  southwestern  Pennsylvania,  embraces  the  banks  of  the  Youghiogheny 
River-a  tributary  of  the  Monongahela  River.  The  river  gauge  has  a  datum  at  860  feet  and  closes  a  drainage 
area  of  1326  square  miles.The  terrain  varies  from  hilly  to  moimtainous,  especially  in  the  eastern  part  of  the 
basin  where  Mt.  Davis  reaches  3213  feet-the  highest  point  in  Pennsylvania. 

Reservoirs 

The  river  flow  in  Connellsville  is  partly  regulated  by  storage  reservoirs.  The  Deep  Creek  Reservoir, 
completed  in  January  1925,  is  used  for  hydro-electric  power  generation.  It  is  owned  an  operated  by  the 
Pennsylvania  Electric  Company.  The  reservoir  has  a  capacity  of  93,000  acre-feet  and  closes  a  drainage  area 
of  65  square  miles,  or  about  5%  of  the  total  basin.  Thus  its  influence  on  flood  flows  at  Connellsville  is 
insignificant. 

The  Youghiogheny  Reservoir,  downstream  of  the  Deep  Creek  Dam,  was  completed  in  October 
1943.  It  serves  multiple  purposes  and  is  operated  by  the  U.S.  Army  Corps  of  Engineers.  The  reservoir  has 
a  capacity  of  254,000  acre-feet,  which  equals  42%  of  the  average  annual  runoff  at  the  dam,  and  controls 
a  drainage  area  of  434  square  miles,  which  constitute  33  %  of  the  total  basin.  The  length  of  the  river 
between  the  dam  and  Connellsville  is  29.4  miles.  All  these  facts  together  suggest  that  the  reservoir  can  only 
partially  control  floods  at  Connellsville. 
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Case  Studies  in  Connellsville 

Two  cases  of  the  flood  forecast  system  were  analyzed  and,  accordingly,  two  sets  of  parameter 
estimates  had  to  be  constructed.  The  first  case  describes  the  present  system  which  is  composed  of  the 
Youghiogheny  Dam  and  the  National  Weather  Service  (NWS)  river  forecasting  technology.  A  flood  forecast 
for  Connellsville  is  prepared  by  routing  the  project  regulated  outflow  from  the  dam  and  superimposing  on 
it  the  predicted  runoff  from  the  drainage  area  between  the  dam  and  the  forecast  point.  The  second  case 
describes  a  hypothetical  system  composed  of  the  NWS  river  forecasting  technology  but  without  any 
influence  of  the  Youghiogheny  Dam  on  flood  flows.  Thus  runoff  from  the  entire  basin  must  be  predicted 
as  flow  at  Connellsville  is  unregulated. 

Data  Records  and  Parameter  Estimation 

Historical  flood  and  forecast  data  were  retrieved  from  the  archives  of  the  National  Weather  Service 
Forecast  Office  in  Pittsburgh.  For  the  present  system,  with  the  Youghiogheny  Dam,  the  prior  distribution 
of  the  flood  crest  was  estimated  from  a  record  spanning  1943-1986.  In  these  44  years,  22  flood  crests 
exceeded  the  official  flood  stage  of  12  feet  (thus,  a  flood  occurred  every  two  years,  on  the  average).  The 
highest  flood  on  record  occurred  in  October  1954  and  reached  22  feet.  The  likelihood  functions  were 
estimated  from  a  historical  joint  record  of  6  forecasted  and  acmal  crests  in  the  period  1984-1986. 

For  the  hypothetical  system,  without  the  Youghiogheny  Dam,  the  prior  distribution  of  flood  crests 
was  estimated  from  a  record  spanning  1910-1942.  During  these  43  years,  22  flood  crests  were  observed 
above  the  flood  stage  of  12  feet.  The  highest  flood  during  that  period  occurred  in  March  1936  and  exceeded 
20  feet.  Estimation  of  the  likelihood  functions  presented  a  challenge  since  there  is  no  historical  joint  record 
of  forecasted  and  acmal  flood  crests-a  record  that  would  correspond  to  the  modem  forecasting  technology 
yet  without  any  influence  of  the  Youghiogheny  Dam.  The  theory  of  sufficient  comparisons  of  forecasts 
systems  [Krzysztofowicz  1992]  came  to  the  rescue  here.  It  seemed  reasonable  to  assume  that  systems 
utilizing  the  same  forecasting  technology  and  operated  for  rivers  with  similar  geomorphologic,  hydrologic, 
and  climatic  characteristics  should  exhibit  similar  statistical  characteristics  of  performance.  In  particular, 
their  standardized  sufficiency  characteristics  (SSC)  should  be  similar.  [For  a  definition  and  properties  of 
the  SSC  see  Krzysztofowicz  1992.]  A  flood  forecast  system  for  Milton,  Pennsylvania,  where  river  flows 
are  unregulated,  was  taken  as  an  analog.  Its  SSC  was  estimated  from  a  historical  joint  record  of  forecasted 
and  acmal  flood  crests;  this  record  contained  8  forecasts  of  floods  that  occurred  in  the  period  1959-1975. 
Next,  the  variance  estimate  in  the  likelihood  functions  for  Connellsville,  the  case  with  the  dam,  was  adjusted 
so  as  to  give  the  SSC  for  Connellsville  the  same  rnagnimde  as  the  SSC  for  Milton.  In  a  sense,  we  have  done 
a  "statistical  transfer"  of  a  forecast  system  from  Milton  to  Connellsville. 

Limitations  of  Models 


Historical  data  records  obtained  from  the  National  Weather  Service  were  used  to  estimate  moments 
of  the  acmal  and  forecasted  crests;  these  estimates  were  next  employed  as  parameter  values  in  the  normal- 
linear  model  of  a  forecast  system.  It  must  be  stressed  that  this  mode,  while  convenient  analytically,  offers 
at  best  an  approximate  representation  of  uncertainties  in  flood  crests  and  their  forecasts.  Moreover,  the 
models  employed  in  this  report  for  computing  probabilities  of  correct  and  false  warnings  are  simplified 
versions  which  ignore  several  significant  sources  of  risk.  (For  example,  they  do  not  account  for  the 
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possibility  of  the  forecasted  crest  s  exceeding  the  threshold  s*,  when  the  actual  crest  h  never  reached  the 
flood  stage-the  case  of  a  false  warning.)  As  a  consequence  of  this  and  other  simplifications,  analyses  based 
on  these  models  cannot  be  taken  as  representative  of  performance  of  flood  forecasts  issued  by  the  National 
Weather  Service  for  Eldred  and  Connellsville.  Rather,  the  results  presented  in  this  report  should  be  treated 
as  hypothetical  examples  having  only  some  (but  not  all)  realistic  features. 

Flood  Damages  in  Connellsville 

Distribution  of  Damages 

The  stage-damage  function  for  Connellsville  was  estimated  according  to  the  methodology  of 
Krzysztofowicz  and  Davis  [1983a,b].  A  crude  inventory  of  establishments  located  at  various  elevations  of 
the  floodplain  was  extracted  from  the  River  Stage  Data  form  prepared  by  the  National  Weather  Service  and 
the  U.S.  Geological  Survey  in  May  1990.  About  212  establishments  were  coimted  in  the  floodplain  and  the 
maximum  possible  damage  for  the  community  was  estimated  to  be  MD  =  $10,400,000  at  the  1991  price 
level. 


In  order  to  construct  the  stage-damage  function,  the  floodplain  was  discretized  into  five  steps, 
whose  elevations  are  listed  in  Table  4-4,  and  the  establishments  were  grouped  intro  three  structural 
categories,  which  are  defined  in  Table  4-5.  The  distribution  {T|(mr);  m  =  1,...,5;  r  =  1,2,3},  partitioning 
the  maximum  possible  damage  MD  among  location  steps  m  and  structural  categories  r,  is  shown  in  Table 
4-6. 


Table  4-4.  Discretization  of  the  Floodplain  in  Connellsville,  Pennsylvania 


Location  step  ra  12345 

Elevation  in  ft  y(m)  12  14  16  18  20 


Flood  stage  is  at  the  elevation  of  12  ft. 
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Table  4-5.  Categories  of  Structures 


r 

Structural  Category 

1 

Two-story  house 

2 

Commercial,  garage  type 

3 

Commercial,  store  type 

Table  4-6.  Distribution  {Ti(mr)}  Partitioning  the  Maximum  Possible  Damage  MD  among 
Location  Steps  m  and  Structural  Categories  r  in  Connellsville,  Pennsylvania 


Location 
Step  m 

Structural  Category  r 

n(m) 

1 

2  3 

1 

.13 

.13 

2 

.18 

.18 

3 

.13 

.13 

4 

.04 

.05  .05 

.14 

5 

.40 

CM 

O 

• 

.42 

n(r) 

• 

00 

00 

.07  .05 

MD  =  $10,400,000 
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Stage-Damage  Function 

For  each  structural  category  r,  there  is  a  unit  damage  function  6r(z)  specifying  the  fraction  of 
maximum  possible  damage  to  an  establishment  which  occurs  when  the  depth  of  flooding,  measured  from 
the  first-floor  level,  is  z.  The  general  form  of  the  unit  damage  function  is  polynomial: 

n 

Sr(z) 

for  0  <  z  <  Zr  and  6r(z)  =  1  for  z  >  ,  where  the  depth  of  flooding  z  is  measured  in  feet.  Table  4-7 

specifies  the  polynomial  coefficients  c^  (i  =  0,...,n)  and  the  domain  limit  7^  for  each  structural  category 
r  =  1,2,3. 

Given  all  this  information,  the  stage-damage  function  D(h)  for  a  community  is  constructed  as 
follows.  For  h  k  y(l), 

M(h)  3 

D(h)  =  MD  £  £  r|(mr)  6r(h  -  y(m)), 

m=l  r=l 

where 

M(h)  =  max  {m  :  h  >  y(m)}. 

For  h  <  y(l),  we  set  D(h)  =  0. 

Application  to  Eldred,  Pennsylvania 

From  the  historical  data,  the  flood  crest,  H,  is  fitted  by  a  normal  distribution,  N(16. 29, 17.38),  and 
the  conditional  probability  density  function  of  the  forecasted  flood  crest,  S,  given  h,  is  fitted  by  a  normal 
distribution,  N(0.63h  -I-  5.57,  1.45).  It  can  be  shown  using  Eqs.  (4.5)-(4.9)  that 

(a)  the  marginal  probability  density  function  of  the  forecast,  k(s),  is  N(15. 8327,8. 348122),  and 

(b)  the  posterior  distribution  density  function  of  h  given  s,  f(h  |  s),  is  N(1.3116s  - 
4.476175,3.018763). 

Substituting  k(s)  and  f(h  |  s)  into  Eqs.  (4.12)  -  (4.15),  the  four  probabilistic  measures  of  the 
forecasting  system  can  be  calculated  for  a  given  warning  threshold,  s*,  and  zone  elevation,  y.  Table  4-8 
shows  those  measures  for  two  values  of  y  and  various  values  of  s*.  Different  values  of  s*  associate  different 
values  of  the  probabilistic  measures,  Pll(s*),  P10(s*),  P01(s*),  and  P00(s*).  They  thus  yield  different 
impacts  on  the  response  fraction  at  the  subsequent  flood  stages. 

The  flood-loss  information  is  not  available  in  the  case  study.  The  following  flood-loss  relationship 
is  assumed.  The  evacuation  cost  for  the  community,  MC,  is  assumed  to  be  equal  to  4  ($10  thousand).  The 
flood-loss  function  without  a  warning  is  assumed  to  be 
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2  Commercial,  garage  type  0  .439931  x  10  ^  11 

1  .707324  X  lO'l 

2  .157361  X  lO'l 

3  -.302723  X  lO"^ 

4  -.576608  X  lO"^ 

5  .978475  x  lO"'^ 

6  .887390  X  10'-'’ 

7  -.143767  X  10"^ 

8  -.476253  x  10'^ 

9  .741636  X  10‘® 

3  Commercial,  store  type  0  .402845  11 

1  .138426 

2  .899010  X  10"^ 

3  -.220052  X  10'2 

4  .506582  X  10"^ 

5  .143909  X  10"^ 

6  -.648618  X  10'^ 
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doi  =  MD6(h  -  y), 


where  MD  =  100  ($10  thousand)  and  6(h  -  y)  equals  0.4(h  -  y)  when  h  -  y  ^  2.5,  and  equals  1  when  h  - 
y  >  2.5  (h  and  y  are  in  the  unit  of  10  feet). 

The  loss  function  with  a  warning  is 

d„  =  MCa  +  MD[1  -  aMR(h  -  y)]6(h  -  y), 

where  MR(h  -  y)  is  assumed  to  be  equal  to  0.25  +  0.4(h  -  y)  -  0.333(h  -  yf  when  h  -  y  ^  1.2,  and  to  be 
equal  to  0.25  when  h  -  y  >  1.2. 

Table  4-9  gives  the  calculated  values  of  the  expected  flood-loss  reduction  for  various  response 
fractions  and  preselected  warning  thresholds  where  y  is  equal  to  16  or  19. 

In  this  case  smdy,  it  is  assumed  that  the  parameters  in  the  response  fraction  dynamic  model,  Eq. 
(4.16),  take  values  of  Ci  =  0.1  and  C2  =  0.9.  A  five-stage  problem  is  considered  with  the  initial  response 
fraction  al  equal  to  0.7.  Two  values  of  the  elevation,  y,  are  considered,  y  =  16  (feet)  and  y  =  19  (feet). 

To  make  the  computational  procedure  feasible,  the  state  space  of  ax  and  the  control  space  of  Sj  are 
quantized  by  the  grid  sizes  0.01  and  0.1,  respectively.  Table  4-10  (a-h)  provides  noninferior  solutions  at 
various  stages  and  various  aj's  for  four  different  weighting  coefficients  and  two  different  values  of  zone 
elevation.  It  can  be  seen  from  Table  4-10  that 

(a)  the  lower  the  weighting  coefficient  0  associated  with  the  first  objective  (loss  reduction),  the 
higher  the  value  of  flood  warning  threshold  will  be  set  to  avoid  possible  high  Type  II  errors; 

(b)  in  order  to  select  a  decision  that  maximizes  the  sum  of  flood-loss  reductions,  the  flood  warning 
threshold  is  set  higher  at  the  earlier  stage  than  at  the  later  stage  (with  respect  to  the  same  value  of  the 
response  fraction)  to  reduce  the  probability  of  high  loss  at  the  later  stages;  and 

(c)  the  higher  the  present  response  fraction,  the  more  cautious  is  the  selection  of  the  threshold.  That 
means  that  a  higher  threshold  value  is  set  for  a  higher  value  of  the  present  response  fraction  to  avoid  losing 
a  larger  number  of  the  response  population.  Note  here  that  the  third  conclusion  may  be  model-specific. 

Application  to  Connellsville,  Pennsylvania 

In  the  case  study  of  Coimellsville,  Pennsylvania,  four  different  situations  with  structural  and 
nonstructural  flood  prevention  measures  are  investigated.  The  expected  flood  losses  in  the  following  four 
cases  are  calculated: 

(1)  expected  flood  loss  in  the  case  with  neither  a  dam  nor  a  flood  warning  system, 

(2)  expected  flood  loss  in  the  case  with  a  dam  and  without  a  flood  warning  system. 
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Table  4-8.  Probabilistic  Measures  of  the  Warning  System  (Eldred,  Pennsylvania) 


y  “  16 


S 

Pll(s*.y) 

Pio(^*.y) 

Poi(^*.y) 

Poo(=*.y) 

14.7 

0.5059415 

0.1465411 

0.0217876 

0.3257299 

14.8 

0.5026233 

0.1369883 

0.0251057 

0.3352826 

14.9 

0.4989391 

0.1276417 

0.0287900 

0.3446293 

15.0 

0.4948698 

0.1185335 

0.0328592 

0.3537375 

15.1 

0.1096941 

0.0373292 

0.3625769 

15.2 

0.4855153 

0.1011519 

0.0422138 

0.3711191 

15.3 

0.4802053 

0.0929329 

0.0475238 

0.3793381 

15.4 

0.4744629 

0.0850599 

0.0532662 

0.3872111 

15.5 

0.4682839 

0.0775526 

0.0594452 

0.3947184 

15.6 

0.4616681 

0.0704275 

0.0660610 

0.4018435 

15.7 

0.4546191 

0.0636972 

0.0731100 

0.4085738 

15.8 

0.4471441 

0.0573709 

0.0805849 

0.4149001 

15.9 

0.4392552 

0.0514533 

0.0884738 

0.4208177 

16.0 

0.4309671 

0.0459459 

0.0967619 

0.4263251 

16.1 

0.4222985 

0.0408466 

0.1054305 

0.4314244 

16.2 

0.4132717 

0.0361495 

0.1144573 

0.4361214 

16.3 

0.4039121 

0.0318458 

0.1238169 

0.4404252 

16.4 

0.3942471 

0.0279234 

0.1334819 

0.4443476 

16.5 

0.3843071 

0.0243680 

0.1434219 

0.4479030 

16.6 

0.3741245 

0.0211626 

0.1536045 

0.4511083 
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Table  4-8.  (continued) 
y  -  19 


ic 

S 

Pll(s*.y) 

Plo(s  ,y) 

Poi(^*.y) 

Poo(^*’y> 

16.5 

0.2443706 

0.1643045 

0.0134605 

0.5778644 

16.6 

0.2423047 

0.1529825 

0.0155265 

0.5891864 

16.7 

0.2400105 

0.1420112 

0.0178206 

0.6001576 

16.8 

0.2374777 

0.1314159 

0.0203534 

0.6107529 

16.9 

0.2346971 

0.1212193 

0.0231341 

0.6209495 

17.0 

0.2316615 

0.1114427 

0.0261697 

0.6307262 

17.1 

0.2283659 

0.1021038 

0.0294652 

0.6400651 

17.2 

0.2248077 

0.0932180 

0.0330235 

0.6489509 

17.3 

0.2209867 

0.0847972 

0.0368444 

0.6573716 

17.4 

0.2169053 

0.0768500 

0.0409259 

0.6653188 

17.5 

0.2125688 

0.0693815 

0.0452624 

0.6727874 

17;6 

0.2079852 

0.0623935 

0.0498459 

0.6797754 

17.7 

0.2031652 

0.0558843 

0.0546659 

0.6862845 

17.8 

0.1981222 

0.0498489 

0.0597090 

0.6923199 

17.9 

0.1928717 

0.0442788 

0.0649595 

0.6978901 

18.0 

0.1874320 

0.0391625 

0.0703992 

0.7030064 

18.1 

0.1818229 

0.0344861 

0.0760083 

0.7076827 

18. -2 

0.1760663 

0.0302327 

0.0817649 

0.7119361 

18.3 

0.1701853 

0.0263839 

0.0876458 

0.7157850 

18.4 

0.1642038 

0.0229186 

0.0936274 

0.7192503 
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Table  4-9.  Expected  Flood-loss  Reduction 

y  =  16 


^  a 

s 

0.50 

0.60 

14.8 

2.274187 

2.729024 

15.0 

2.281817 

2.738181 

l5.2 

2.278504 

2.734205 

15.4 

2.262371 

15.6 

2.231861 

2.678233 

15.8 

2.185879 

2.623055 

16.0 

2.123912 

2.548693 

16.2 

2.046122 

2.455346 

16.4 

1.953364 

2.344036 

16.6 

1.847164 

2.216596 

y  -  19 

★  a 

s 

0.50 

0.60 

16.6 

0.4732468 

0.5678961 

16.8 

0.4981555 

0.5977865 

17.0 

0.5155946 

0.6187135 

17.2 

0.5250767 

0.6300920 

17.4 

0.5263813 

0.6316575 

17.6 

0.5196016 

0.6235219 

17.8 

0.5051445 

0.6061733 

18.0 

0.4837199 

0.5804638 

18.2 

0.4563089 

0.5475706 

18.4 

0.4240967 

0.5089160 

3.124605 


3.060230 


2.973476 


2.864571 


2.734709 


2.586029 


3.638699 

3.650908 

3.645607 


3.619794 


3.5709.77 


3.497406 


3.398258 


3.273795 


3.125381 


2.955462 


0.70 


0.6625457 


0.7218326 


0.7351076 


0.7274424 


0.7072024 


0.6772079 


0.6388326 


0.80 


0.7571951 


0.8249515 


0.8401229 


0.8313628 


0.8082313 


0.7739519 


0.7300944 
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Table  4-lOa.  Optimal  Flood  Warning  Threshold  (Eldred,  Pennsylvania) 


Y  -  16  ^  -  0.02 


a 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

15.2 

0.6962703 

15.2 

0.6313940 

15.2 

0.5671985 

15.2 

0.5037215 

0.50 

15.4 

0.7515748 

15.3 

0.6852435 

15.3 

0.6195143 

15.3 

0.5544193 

0.55 

15.4 

0.806968A 

15.3 

0.7391637 

15.3 

0.6718839 

15.3 

0.6051552 

0.60 

15.5 

0.8624514 

15.4 

0.7931497 

15.4 

0.7242879 

15.4 

0.6558943 

0.65 

15.5 

0.9179950 

15.4 

0.8471968 

15.4 

0.7767507 

15.4 

0.7066774 

0.70 

15.5 

1.0461380 

15.5 

0.9735616 

15.4 

0.9012520 

15.4 

0.8292162 

15.4 

0.7574605 

0.75 

15.6 

1.0292320 

15.5 

0.9553727 

15.5 

0.8817053 

15.4 

0.8082436 

0.80 

15.6 

1.084935 

15.5 

1.009541 

15.5 

0.9342493 

15.5 

0.8590660 

The  upper  number  denotes  the  selected  flood  warning  threshold, 
the  lower  number  denotes  the  expected  loss  reduction  from 
the  present  stage  to  the  final  stage. 
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Table  4-lOb.  Optimal  Flood  Warning  Threshold  (Eldred,  Peimsylvania) 


Y  -  16  S  -  0.06 


a 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

15.2 

1.026391 

15.1 

0.8688733 

15.1 

0.7156359 

15.1 

0.5669228 

0.50 

15.2 

1.113314 

15.2 

0.9466786 

15.2 

0.7837586 

15.1 

0.6247922 

0.55 

15.2 

1.200264 

15.2 

1.024525 

15.2 

0.8519347 

15.2 

0.6826890 

0.60 

15.2 

1.287223 

15.2 

1.102380 

15.2 

0.9201196 

15.2 

0.7406027 

0.65 

15.3 

1.374238 

15.2 

1.180238 

15.2 

0.9883067 

15.2 

0.7985163 

0.70 

15.3 

1.6659270 

15.3 

1.461304 

15.3 

1.258129 

15.2 

1.056494 

15.2 

0.8564301 

0.75 

15.3 

1.548389 

15.3 

1.336045 

15.2 

1.124681 

15.2 

0.9143438 

0.80 

15.3 

1.635497 

15.3 

1.413986 

15.3 

1.192896 

15.2 

0.9722574 

The  upper  number  denotes  the  selected  flood  warning  threshold; 
the  lower  number  denotes  the  expected  loss  reduction  from 
the  present  stage  to  the  final  stage. 
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Table  4-lOc.  Optimal  Flood  Warning  Threshold  (Eldred,  Pennsylvania) 


Y  -  16  ^  -  0.1 


a 

Stage  1 

Stage  3 

Stage  4 

Stage  5 

0.45 

15.1 

1.356822 

15.1 

1.106607 

15.1 

0.8642271 

15.1 

0.6301808 

0.50 

15.2 

1.475427 

15.1 

1.208406 

15.1 

0.9482521 

15.1 

0.6952969 

0.55 

15.2 

1.594128 

15.2 

1.310294 

15.1 

1.032277 

15.1 

0.7604131 

0.60 

15.2 

1.712853 

15.2 

1.412213 

15.2 

1.116337 

15.1 

0.8255290 

0.65 

15.2 

1.831588 

15.2 

1.514143 

15.2 

1.200418 

15.1 

0.8906451 

0.70 

15.3 

2.2871940 

15.2 

1.950328 

15.2 

1.616077 

15.2 

1.284499 

15.1 

0.9557613 

0.75 

15.3 

2.069109 

15.2 

1.718019 

15.2 

1.368585 

15.1 

1.020877 

0.80 

15.3 

2.187940 

15.2 

1.819972 

15.2 

1.452688 

15.2 

1.086017 

The  upper  number  denotes 
the  lower  number  denotes 
the  present  stage  to  the 


the  selected  flood  warning  threshold; 
the  expected  loss  reduction  from 
final  stage. 
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Table  4-lOd.  Optimal  Flood  Warning  Threshold  (Eldred,  Pennsylvania) 


Y  -  16  6-0.5 


a 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

15.1 

4.662013 

15.1 

3.484156 

15.1 

2.350139 

15.1 

1.262761 

0.50 

15.1 

5.098404 

15.1 

3.827040 

15.1 

2.593547 

15.1 

1.400343 

0.55 

15.1 

5.534795 

15.1 

4.169925 

15.1 

2.836957 

15.1 

1.537926 

0.60 

15.2 

5.971472 

15.1 

4.512809 

15.1 

3.080365 

15.1 

1.675508 

0.65 

15.2 

6.408165 

15.1 

4.855693 

15.1 

3.323773 

15.1 

1.813090 

0.70 

15.2 

8.5052660 

15.2 

6.844906 

15.1 

5.198578 

15.1 

3.567182 

15.1 

1.950673 

0.75 

15.2 

7.281742 

15.2 

5.541626 

15.1 

3.810591 

15.1 

2.088255 

0.80 

15.2  ^ 

7.718588 

15.2 

5.884681 

15.1 

4.054000 

15.1 

2.225837 

The  upper  number  denotes  the  selected  flood  warning  threshold, 
the  lower  number  denotes  the  expected  loss  reduction  from 
the  present  stage  to  the  final  stage. 
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Table  4-lOe.  Optimal  Flood  Warning  Threshold  (Eldred,  Pennsylvania) 


Y  -  19  5-0.02 


a 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

17.5 

0.5120240 

17.5 

0.4943402 

17.5 

0.4766091 

17.5 

0.4588294 

0.50 

17.6 

0.5598703 

17.6 

0.5424591 

17.6 

0.5250116 

17.6 

0.5075260 

0.55 

17.6 

0.6078050 

17.6 

0.5906430 

17.6 

0.5734548 

17.6 

0.5562404 

0.60 

17.7 

0.6558577 

17.7 

0.6389248 

17.7 

0.6219716 

17.7 

0.6049971 

0.65 

17.7 

0.7039466 

17.7 

0.6872281 

17.7 

0.6704972 

17.7 

0.6537544 

0.70 

17.8 

0.7686545 

17.8 

0.7521415 

17.8 

0.7356198 

17.8 

0.7190891 

17.8 

0.7025492 

0.75 

17.8 

0.8003519 

17.8 

0.7840200 

17.8 

0.7676843 

17.8 

0.7513444 

0.80 

17.9 

0.8486295 

17.9 

0.8324716 

17.9 

0.8163140 

17.9 

0.8001567 

The  upper  number  denotes  the  selected  flood  warning  threshold; 
the  lower  number  denotes  the  expected  loss  reduction  from 
the  present  stage  to  the  final  stage. 
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Table  4-lOf.  Optimal  Flood  Warning  Threshold  (Eldred,  Pennsylvania) 


Y  -  19  B  ~  0.06 


a 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

17.4 

0.5703409 

17.4 

0.5330652 

17.4 

0.4960785 

17.4 

0.4593878 

0.50 

17.5 

0.6242197 

17.5 

0.5852870 

17.5 

0.5466008 

17.5 

0.5081689 

0.55 

17.5 

0.6781828 

17.5 

0.6375848 

17.5 

0.5971845 

17.5 

0.5569876 

0.60 

17.5 

0.7321548 

17.5 

0.6898881 

17.5 

0.6477705 

17.5 

0.6058063 

0.65 

17.5 

0.7861278 

17.5 

0.7421919 

17.5 

0.6983567 

17.5 

0.6546251 

0.70 

17.6 

0.8857769 

17.5 

0.8401028 

17.5 

0.7944958 

17.5 

0.7489429 

17.5 

0.7034438 

0.75 

17.6 

0.8941438 

17.6 

0.8468298 

17.6 

0.7995355 

17.5 

0.7522625 

0.80 

17.6 

0.9482240 

17.6 

0.8992121 

17.6 

0.8501706 

17.6 

0.8108692 

The  upper  number  denotes  the  selected  flood  warning  threshold; 
the  lower  number  denotes  the  expected  loss  reduction  from 
the  present  stage  to  the  final  stage. 
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Table  4-lOg.  Optimal  Flood  Warning  Threshold  (Eldred,  Pennsylvania) 


Y  -  19  ^  -  0.5 


a 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

17.4 

1.215052 

17.4 

0.9613236 

17.3 

0.7116798 

17.3 

0.4662075 

0.50 

17.4 

1.337938 

17.4 

1.060809 

17.4 

0.7870544 

17.3 

0.5167807 

0.55 

17.4 

1.460841 

17.4 

1.160317 

17.4 

0.8624612 

17.3 

0.5673540 

0.60 

17.4 

1.583747 

17.4 

1.259827 

17.4 

0.9378675 

17.3 

0.6179269 

0.65 

17.4 

1.706653 

17.4 

1.359338 

17.4 

1.013274 

17.3 

0.6685001 

0.70 

17.4 

2.2007990 

1 

17.4 

1.829560 

17.4 

1.458848 

17.4 

1.088681 

17.3 

0.7190733 

0.75 

17.4 

1.952466 

17.4 

1.558359 

17.4 

1.164087 

17.3 

0.7696464 

0.80 

17.4 

2.075372 

17.4 

1.657869 

17.4 

1.239494 

17.4 

0.8202195 

The  upper  number  denotes 
the  lower  number  denotes 
the  present  stage  to  the 


the  selected  flood  warning  threshold; 
the  expected  loss  reduction  from 
final  stage. 
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Table  4-101i.  Optimal  Flood  Warning  Threshold  (Eldred,  Pennsylvania) 


Y  -  19  e  -  0.9 


a 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

17.4 

1.801429 

17.3 

1.350909 

17.3 

0.9078862 

17.3 

0.4725028 

0.50 

17.4 

1.987141 

17.4 

1.493429 

17.3 

1.005951 

17.3 

0.5247576 

0.55 

17.4 

2.172883 

17.4 

1.635987 

17.3 

1.104016 

17.3 

0.5770127 

0.60 

17.4 

2.358627 

17.4 

1.778545 

17.3 

1.202080 

17.3 

0.6292672 

0.65 

17.4 

2.544372 

17.4 

1.921103 

17.3 

1.300145 

17.3 

0.6815220 

0.70 

17.4 

3.3975510 

17.4 

2.730119 

17.4 

2.063663 

17.3 

1.398209 

17.3 

0.7337769 

0.75 

17.4 

2.915882 

17.4 

2.206238 

17.4 

1.496290 

17.3 

0.7860316 

0.80 

17.4 

3.101650 

17.4 

2.348823 

17.4 

1.594384 

17.3 

0.8382865 

The  upper  number  denotes 
the  lower  number  denotes 
the  present  stage  to  the 


the  selected  flood  warning  threshold; 
the  expected  loss  reduction  from 
final  stage. 
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(3)  expected  flood  loss  in  the  case  without  a  dam  and  with  a  flood  warning  system,  and 

(4)  expected  flood  loss  in  the  case  with  both  a  dam  and  a  flood  warning  system. 


When  there  is  a  dam,  from  the  historical  data  the  flood  crest,  H,  is  fitted  by  a  prior  normal 
distribution,  N(13.72,5.06),  and  the  conditional  probability  density  function  of  the  forecasted  flood  crest, 
S,  given  h,  is  fitted  by  a  normal  distribution,  N(1.16h  -  1.77,  0.45).  It  can  be  shown  using  Eqs.  (4.5)-(4.9) 
that  (a)  the  marginal  probability  density  function  of  the  forecast,  k(s),  is  N(14.1452,  7.258735)  and  (b)  the 
posterior  distribution  density  fonction  of  h  given  s,  f(h  |  s),  is  N(0.808626s  -  2.281829,  0.313691). 

When  there  is  no  dam,  from  the  historical  data  the  flood  crest,  H,  is  fitted  by  a  prior  normal 
distribution  N(14.53,5.14)  and  the  conditional  probability  density  function  of  the  forecasted  flood  crest,  S, 
given  h,  is  fitted  by  a  normal  distribution,  N(1.16h  -  1.77,  3.51).  It  can  be  shown  using  Eqs.  (4.5)-(4.9) 
that  (a)  the  marginal  probability  density  function  of  the  forecast,  k(s),  is  N(15.0848,  10.42638)  and  (b)  the 
posterior  distribution  density  function  of  h  given  s,  f(h  |  s),  is  N(0. 571857s  -  5.903653,  1.73036). 

Substituting  k(s)  and  f(h  |  s)  into  Eqs.  (4.12)-(4.15),  the  four  probabilistic  measures  of  the 
forecasting  system  can  be  calculated  for  a  given  warning  threshold,  s*,  and  zone  elevation,  y.  Table  4-11 
shows  those  measures  for  two  values  of  y  and  various  values  of  s*  for  both  cases,  with  a  dam  and  without 
a  dam.  Different  values  of  s*  are  associated  with  different  values  of  the  probabilistic  measures,  Pii(s*), 
Pio(s*),  Poi(s*),  and  Poo(s*).  They  thus  yield  different  impacts  on  the  response  fraction  at  the  subsequent 
flood  events. 

From  the  historical  data,  the  unit  damage  function  is  fitted  in  this  case  study  [unit  damage  function 
for  a  two-story  house]  by 

6(h  -  y)  =  0.110007  +  0.271 166(h  -  y)  -1-  0.137889  (h  -  y)^ 

-  0.0399962(h  -  y)^  -  0.00032665(h  -  y)^when  h  -  y  ^  2.4, 

and  equals  1  when  h  -  y  >  2.4  (h  and  y  are  in  the  unit  of  10  feet).  The  other  parameters  in  this  model  are 
assumed  as  follows. 

The  evacuation  cost  for  the  community,  MC,  is  assumed  to  be  equal  to  4  ($10  thousand). 

The  flood-loss  function  without  a  warning  is  assumed  to  be 

do,  =  MD6(h  -  y), 
where  MD  =  100  ($10  thousand). 

The  loss  function  with  a  warning  is 
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dll  =  MCa  +  MD[1  -  aMR(h  -  y)]6(h  -  y), 

where  MR(h  -  y)  is  assumed  to  be  equal  to  0.25  +  0.4(h  -  y)  -  0.333(h  -  y)^  when  h  -  y  ^  1.2,  and  equal 
to  be  equal  to  0.25  when  h  -  y  >  1.2. 

Table  4-12  gives  the  calculated  values  of  the  expected  flood  loss  without  a  warning  system  for  y 
equal  to  12  and  14  in  both  cases,  with  a  dam  and  without  a  dam.  Table  4-11  also  gives  the  calculated  values 
of  the  expected  flood-loss  reduction  with  full  response  for  various  preselected  warning  thresholds  for  y 
equal  to  12  and  14  in  both  cases,  with  a  dam  and  without  a  dam. 

In  this  case  study,  it  is  assumed  that  the  parameters  in  the  response  fraction  dynamic  model,  Eq. 
(4.16),  take  values  of  Cj  =  0.1  and  C2  =  0.9.  A  five-stage  problem  is  considered  with  the  initial  response 
fraction  aj  equal  to  0.7.  Two  values  of  the  elevation,  y,  are  considered,  y  =  12  (feet)  and  y  =  14  (feet). 

To  make  the  computational  procedure  feasible,  the  state  space  of  Oj  and  the  control  space  of  Sjare 
quantized  by  the  grid  sizes  0.01  and  0.1,  respectively.  Table  4-13  (a-h)  provides  noninferior  solutions  at 
various  stages  and  various  aj's  for  two  different  weighting  coefficients  and  two  different  values  of  zone 
elevation  for  both  cases,  with  a  dam  and  without  a  dam.  It  can  be  seen  from  Table  4-13  that 

(a)  the  lower  the  weighting  coefficient,  0,  associated  with  the  first  objective  (loss  reduction),  the 
higher  the  value  of  the  flood  warning  threshold  will  be  set  to  avoid  possible  high  Type  II  errors; 

(b)  in  order  to  select  a  decision  that  maximizes  the  sum  of  flood-loss  reductions,  the  flood  warning 
threshold  is  set  higher  at  the  earlier  stage  than  at  the  later  stage  (with  respect  to  the  same  value  of  the 
response  fraction)  to  reduce  the  probability  of  high  loss  at  the  later  stages;  and 

(c)  the  higher  the  present  response  fraction,  the  more  cautious  is  the  selection  of  the  threshold.  That 
means  that  a  higher  threshold  value  is  set  for  a  higher  value  of  the  present  response  fraction  to  avoid  losing 
a  larger  number  of  the  response  population.  Note  here  that  the  third  conclusion  may  be  model-specific. 
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Table  4-11.  Probabilistic  Measures  of  the  Warning  System  (Connellsville,  Pennsylvania) 


y  -  12  (with  a  dam) 


k 

k 

k 

k 

S 

.y) 

?^q(s  .y) 

Poi(s  .y) 

Pqo^®  -y) 

11.0 

0.7763117 

0.1021628 

0.0014434 

0.1200821 

0.0879903 

0.0029068 

0.1342546 

11.4 

0.7722843 

0.0735971 

0.0054707 

0.1486478 

11.6 

0.7681055 

0.0594870 

0.0096495 

0.1627579 

11.8 

0.7617464 

0.0462294 

0.0160087 

0.1760156 

12.0 

0.7526836 

0.0343665 

0.0250714 

0.1878784 

12.2 

0.7405316 

0.1979255 

12.4 

0.7251226 

0.0163076 

0.0526324 

0.2059374 

12.6 

0.7065377 

0.0103184 

0.0712173 

0.2119266 

12.8 

0.6850762 

0.0061376 

0.0926788 

0.2161073 

14  (with  a  dam) 
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Table  4-11.  (continued) 
y  -  12  (without  a  darn) 


s 

Pio(=*.y) 

9.6 

0.8598625 

0.0954413 

9.8 

10.0 

0.8552552 

0.0870857 

10.2 

0.8522224 

0.0826117 

10.4 

0.8486115 

0.0779786 

10.6 

0.8443505 

0.0732205 

10.8 

0.0634843 

11.2 

0.8269404 

0.0585908 

11.4 

0.8193607 

0.0537378 

Poi(s*.y) 

^00^^ 

0.0079126 

0.0367837 

0.0099954 

0.0408541 

0.0125198 

0.0451392 

0.0155527 

0.0496133 

0.0191636 

0.0542464 

0.0234246 

0.0590045 

0.0284086 

0.0638494 

0.0341881 

0.0687406 

0.0408347 

0.0736342 

0.0484143 

0.0784871 

y  -  14  (without  a  dam) 


* 

s 

Pii(s*.y) 

Pio(s*,y) 

Poi(»*.y) 

Poo(5  .y) 

13.0 

0.5588912 

0.1818558 

0.0335281 

0.2257249 

13.2 

0.5522861 

0.1680066 

0.0401332 

0.2395740 

13.4 

0.5447636 

0.1543220 

0.0476557 

0.2532586 

13.6 

0.5362751 

0.1409070 

0.0561443 

0.2666736 

13.8 

0.5267833 

0.1278629 

0.0656361 

0.2797178 

14.0 

0.5162637 

0.1152843 

0.0761557 

0.2922964 

14.2 

0.5047076 

0.1032567 

0.0877118 

0.3043239 

14.4 

0.4921228 

0.0918543 

0.1002966 

0.3157264 

14.6 

0.4785339 

0.0811387 

0.1138854 

0.3264419 

14.8 

0.4639837 

0.0711576 

0.1284356 

0.3364230 
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Table  4-12.  Expected  Flood-loss  Reduction  with  Full  Response  (Connellsville,  Pennsylvania) 

WITH  A  DAM 
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Table  4-13a.  Optimal  Flood  Warning  Thresholds  (Connellsville,  Pennsylvania,  with  a  Dam) 


Y  -  12  -  0.02 


a 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

11.9 

0.6457910 

11.9 

0.5969478 

11.9 

0.5466338 

11.9 

0.4947060 

0.50 

12.0 

0.6870015 

12.0 

0.6398434 

12.0 

0.5913768 

12.0 

0.5414749 

0.55 

12.1 

0.7283248 

12.1 

0.6828295 

12.1 

0.6361857 

12.1 

0.5882831 

0.60 

12.2 

0.7697667 

12.2 

0.7259145 

12.2 

0.6810647 

12.2 

0.6351289 

0.65 

12.2 

0.8113291 

12.3 

0.7690976 

12.3 

0.7260182 

12.3 

0.6820114 

0.70 

12.3 

0.8930091 

12.3 

0.8530287 

12.3 

0.8123971 

12.3 

0.7710563 

12.3 

0.7289464 

0.75 

12.4 

0.8948482 

12.4 

0.8557923 

12.4 

0.8161666 

12.4 

0.7759222 

0.80 

12.5 

0.9367898 

12.5 

0.8992963 

12.5 

0.8613571 

12.5 

0.8229374 

The  upper  number  denotes  the  selected  flood  warning  threshold; 
the  lower  number  denotes  the  expected  loss  reduction  from 
the  present  stage  to  the  final  stage. 
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Table  4-13b.  Optimal  Flood  Warning  Thresholds  (Connellsville,  Pennsylvania,  with  a  Dam) 


Y  =  12  ^  =  0.50 


a 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

12.0 

1.945561 

12.0 

1.473788 

12.0 

1.028254 

12.1 

0.6111944 

0.50 

12.1 

2.108311 

12.1 

1.606186 

12.1 

1.127736 

12.1 

0.6749548 

0.55 

12.1 

2.271216 

12.1 

1.738675 

12.1 

1.227249 

12.1 

0.7387153 

0.60 

12.1 

2.434124 

12.1 

1.871165 

1 

12.1 

1.326761 

12.1 

0.8024757 

0.65 

12.2 

2.597144 

12.1 

2.003655 

12.1 

1.426274 

12.1 

0.8662361 

0.70 

12.2 

3.3968910 

12.2 

2.760314 

12.2 

2.136265 

12.1 

1.525787 

12.1 

0.9299966 

0.75 

12.2 

2.923496 

12.2 

2.268907 

12.2 

1.625355 

12.1 

0.9937571 

0.80 

12.2 

3.086682 

12.2 

2.401549 

12.2 

1.724926 

12.1 

1.057517 

The  upper  number  denotes  the  selected  flood  warning  threshold; 
the  lower  number  denotes  the  expected  loss  reduction  from 
the  present  stage  to  the  final  stage. 
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Table  4-13c.  Optimal  Flood  Warning  Thresholds  (Connellsville,  Pennsylvania,  with  a  Dam) 


Y  -  14  ^  -  0.02 


a 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

14.4 

0.5520497 

14.4 

0.5250422 

14.4 

0.4975516 

14.4 

0.4695483 

0.50 

14.5 

0.5962349 

14.5 

0.5703243 

14.5 

0.5439987 

14.5 

0.5172327 

0.55 

14.6 

0.6405975 

14.6 

0.6157462 

14.6 

0.5905436 

14.6 

0.5649695 

0.60 

14.6 

0.6851410 

14.6 

0.6613101 

14.7 

0.6371866 

14.7 

0.6127548 

0.65 

14.7 

0.7298880 

14.7 

0.7070405 

14.7 

0.6839515 

14.7 

0.6606089 

0.70 

14.8 

0.7965054 

14.8 

0.7747915 

14.8 

0.7528939 

14.8 

0.7308034 

14.8 

0.7085098 

0.75 

14.9 

0.8198694 

14.9 

0.7988837 

14.9 

0.7777481 

14.9 

0.7564557 

0.80 

14.9 

0.8651059 

14.9 

0.8449958 

15.0 

0.8247769 

15.0 

0.8044469 

The  upper  number  denotes  the  selected  flood  warning  threshold 
the  lower  number  denotes  the  expected  loss  reduction  from 
the  present  stage  to  the  final  stage. 
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Table  4-13d.  Optimal  Flood  Warning  Thresholds  (Connellsville,  Pennsylvania,  with  a  Dam) 


Y  -  14  6  ^  0.5 


a 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

14.5 

1.133536 

14.5 

0. 8937591 

14.5 

0.6621569 

14.5 

0.4391043 

0.50 

14.5 

1.239075 

14.5 

0.9805041 

14.5 

0.7292412 

14.6 

0.4856232 

0.55 

14.6 

1.344762 

14.6 

1.067357 

14.6 

0.7963985 

14.6 

0.5321780 

0.60 

14.6 

1.450539 

14.6 

1.154261 

14.6 

0.8635748 

14.6 

0.5787328 

0.65 

14.6 

1.556319 

14 . 6 

1.241166 

14.6 

0.9307513 

14 . 6 

0.6252878 

0.70 

14 . 6 

1.9998400 

14.6 

1.662099 

14.6 

1.328072 

14.6 

0.9979277 

j 

14.6 

0.6718426 

0.75 

14.7 

1.767897 

14.6 

1.414977 

14.6 

1.065104 

14.6 

0.7183975 

0.80 

,14.7 

1.873788 

14.6 

1.501882 

14.6 

1.132280 

14.6 

0.7649522 

The  upper  ntunber  denotes  Che  selected  flood  warning  threshold; 
the  lower  number  denotes  the  expected  loss  reduction  from 
the  present  stage  to  the  final  stage. 
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Table  4-13e.  Optimal  Flood  Warning  Thresholds  (Connellsville,  Pennsylvania,  without  a  Dam) 


Y  -  12  -  0.06 


a 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

10.2 

0.8719684 

10.2 

0.7561035 

10.2 

0.6425062 

10.2 

0.5313926 

0.50 

10.3 

0.9300157 

10.3 

0.8117548 

10.3 

0.6955323 

10.3 

0.5815435 

0.55 

10.5 

0.9881452 

10.5 

0.8674700 

10.5 

0.7486029 

10.4 

0.6317162 

0.60 

10.6 

1.046350 

10.6 

0.9232452 

10.6 

0.8017159 

10.6 

0.6819128 

0.65 

10.8 

1.104629 

10.7 

0.9790777 

10.7 

0.8548683 

10.7 

0.7321298 

0.70 

10.9 

1.2920150 

10.9 

1.162982 

10.8 

1.034964 

10.8 

0.9080573 

10.8 

0.7823657 

0.75 

11.0 

1.221405 

11.0 

1.090904 

10.9 

0.9612833 

10.9 

0.8326194 

0.80 

11.1 

1.279884 

11.1 

1.146892 

11.0 

1.014544 

11.0 

0.8828902 

The  upper  number  denotes  the  selected  flood  warning  threshold 
the  lower  number  denotes  the  expected  loss  reduction  from 
the  present  stage  to  the  final  stage. 
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Table  4-13f.  Optimal  Flood  Warning  Thresholds  (Connellsville,  Pennsylvania,  without  a  Dam) 


Y  -  12  5=0.5 


a 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

10.1 

2.866462 

10.1 

2.137300 

10.0 

1.450978 

10.0 

0.8119178 

0.50 

10.2 

3.101412 

10.1 

2.327197 

10.1 

1.591223 

10.0 

0.8973795 

0.55 

10.3 

3.336430 

10.2 

2.517138 

10.1 

1.731475 

10.0 

0.9828411 

0.60 

10.3 

3.571528 

i 

10.3 

2.707108 

10.2 

1.871732 

10.0 

1.068303 

0.65 

10.4 

3.806681 

10:3 

2.897112 

10.2 

2.012008 

10.0 

1.153764 

0.70 

10.5 

5.0148220 

10.5 

4.041880 

10.4 

3.087120 

10.2 

2.152283 

10.0 

1.239226 

0.75 

10.5 

4.277140 

10.4 

3.277164 

10.2 

2.292558 

10.0 

1.324688 

0.80 

10.6 

4.512427 

10.4 

3.467223 

10.3 

2.432843 

10.0 

1.410149 

The  upper  number  denotes  the  selected  flood  warning  threshold; 
the  lower  number  denotes  the  expected  loss  reduction  from 
the  present  stage  to  the  final  stage. 
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Table  4-13g.  Optimal  Flood  Warning  Thresholds  (Connellsville,  Pennsylvania,  without  a  Dam) 


Y  -  14  ^  0.06 


a 

Stage  1 

Stage  2 

Stage  3 

Stage  4 

Stage  5 

0.45 

13.7 

0.6979931 

13.7 

0.6279304 

13.7 

0.5587157 

13.7 

0.4903911 

0.50 

13.9 

0.7524474 

13.9 

0.6806543 

13.9 

0.6096215 

13.9 

0.5393894 

0.55 

14.1 

0.8071754 

14.1 

0.7335954 

14.1 

0.6606817 

14.1 

0.5884703 

0.60 

14.3 

0.8621708 

14.2 

0.7867416 

14.2 

0.7118861 

14.2 

0.6376284 

0.65 

14.4 

0.9174379 

14.4 

0.840098 

14.4 

0.7632285 

14.4 

0.6868537 

0.70 

14.6 

1.0525830 

14.5 

0.9729285 

' 

14.5 

0.8936280 

14.5 

0.8146937 

14.5 

0.7361451 

0.75 

14.7 

1.028641 

14.7 

0.9473197 

14.6 

0.8662665 

14.6 

0.7854905 

0.80 

14.8 

1.084541 

14.8 

1.001158 

14.7 

0.9179361 

14.7 

0.8348857 

The  upper  number  denotes  the  selected  flood  warning  threshold 
the  lower  number  denotes  the  expected  loss  reduction  from 
the  present  stage  to  the  final  stage. 
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W  denotes  plans  incorporating  flood  warning  systems 

y  elevation  of  the  floodplain  zone  under  consideration 

a  nonexceedance  probability  that  partitions  the  range  of  extreme  events;  used  in  the  definition 

of  the  conditional  expected  value  f^ 

6(h  -  y)  unit  damage  function  specifying  the  fraction  of  MD  that  occurs  when  the  depth  of  flooding 

is  (h  -  y) 

6  fraction  of  the  community  that  responds  to  a  flood  warning;  response  fraction 


Part  2.  Multiobjective  Decision-Tree  Analysis 


an 

C 

C1,C2,C3 

dj™ 

DN1,DN2 

E[*] 

E'[] 

EV1,EV2 

EVE 

f. 

f2,f3 

f4 

f4* 


action,  or  alternative,  or  option,  at  a  decision  node  n 

maximum  possible  loss  of  property  (discrete  case);  possible  loss  of  lives  given  no  flood 
warning  -  linear  function  of  discharge  W  (continuous  case) 

chance  nodes  in  the  decision  tree 

number  of  elements  in  the  the  set  rj™ 

do-nothing  option  in  the  first  and  second  decision  periods,  respectively 
expected  value 

the  5th  averaging-out  strategy;  for  example  E'*  denotes  the  conditional  expected  value  of 
extreme  events  f4 

evacuation  order  in  the  first  and  second  decision  periods,  respectively 

expected  value  of  experimentation;  difference  between  expected  loss  without 
experimentation  and  expected  loss  with  experimentation 

standard  normal  distribution  function 

cost  objective  function;  balanced  with  the  risk  functions  fjthru  fj  in  the  PMRM 
conditional  expected  values 

conditional  expected  value  of  the  (damage)  risk  of  extreme  events 
optimal  value  of  f,,  see  Equation  (2.17) 
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fs 

k 

L 

LN 

Px 

Px 


W 

WA1,WA2 

X 

a 

Ki 

Pu 

All 


o 


overall  expected  value  of  damage 
dimension  of  the  objective  function  vector 

maximum  possible  loss  of  lives  (discrete  case);  possible  loss  of  lives  given  no  flood 
warning  —  linear  function  of  discharge  W  (continuous  case) 

lognormal  distribution 

cumulative  distribution  function  of  X 

probability  density  function  of  X 

the  vector  of  objective  functions  in  the  decision  tree  [r,,...,rj 

set  of  Pareto  optimum  alternatives  associated  with  each  branch  emerging  from  chance  node 
m 

actual  flood  level  (cfs) 

issuing  a  flood  watch  in  the  first  and  second  decision  periods,  respectively 
random  variable  of  damage  or  loss 

partitioning  nonexceedance  for  the  conditional  expected  value 

values  of  nonexceedance  probability  that  partition  the  ranges  of  risk  in  the  PMRM 

values  of  damage  that  partition  the  severity  of  risk  in  the  PMRM  for  the  7th  policy 

tradeoffs  between  the  cost  objective  function  and  the  ith  risk  function 

mean  of  the  discharge  W 

state  of  nature  at  node  n  of  the  decision  tree  (also  used  in  unrelated  context  as  parameters 
in  the  PMRM  defined  by  Equation  2.4) 

standard  deviation  of  the  discharge  W 


Part  3.  Performance  Characteristics  of  a  Flood  Warning  System 


a,b  parameters  of  the  normal-linear  likelihood  model  f 
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D 

F 

f(s  I  h,0=l,T=l) 
FSC 

g(h  I  0  =  1) 
g(A  I  0  =  1) 
h 

LT 

M 

N 

N 

n 

ND 

NF 

PTC 

q 

q(s) 

q* 

ROC 

s 

T 

W 


detection  (Equation  3.7) 
false  warning  (Equation  3.7) 

probability  density  of  s  conditional  on  the  actual  crest  h,  0  =  1,  T  =  1 

forecast  sufficiency  characteristic,  a  measure  sufficient  for  comparing  any  two 
forecasters  who  produce  forecasts  of  the  same  variate 

prior  probability  density  function  of  flood  crest  given  flood  occurs 

probability  density  of  X  conditional  on  0  =  1 

height  of  actual  flood  crest 

expected  lead  time 

missed  flood  (Equation  3.7) 

expected  number  of  floods  per  year 

normal  probability  distribution 

expected  number  of  zone  floods  per  year 

expected  number  of  detections  per  year  for  a  zone 

expected  number  of  false  warnings  per  year  for  a  zone 

performance  tradeoff  characteristic,  a  plot  of  ND  versus  NF 

quiet  (Equation  3.7) 

P(0  =  1  I  s,  T  =  1),  posterior  probability  of  a  flood  in  a  given  zone 

optimal  threshold  associated  with  warning  rule  W* 
relative  operating  characteristic,  a  plot  of  P(D)  versus  P(F) 

forecasted  flood  crest 

trigger  indicator:  trigger  is  not  observed  (T  =  0),  trigger  is  observed  (T  =  1) 

warning  rule,  w  =  W(s),  where  w  =  0  and  w  =  1  denote  "do  not  issue  warning" 
and  "do  issue  warning,"  respectively 
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W*(s)  optimal  warning  rule  (of  the  threshold  type)  minimizes  expected  disutility  of 

outcomes 

y  zone  elevation 

Y  P(0  =  1  1  T  =  1),  diagnosticity  conditional  probability 

Ko(s  I  ©=0,T=1)  probability  density  of  s  conditional  on  the  forecast  0  =  0  and  T  =  1 

A  lead  time  of  a  warning  for  a  given  zone,  conditional  on  hypothesis  that  zone  will  be 

flooded 

/i^,0h  mean  and  standard  deviation  of  the  prior  density  g(h  |  0  =  1) 

jUj.Oj  mean  and  standard  deviation  of  the  likelihood  function  ko 

0  zone  flood  indicator:  zone  flood  does  not  occur  (0  =  0),  zone  flood  occurs  (0  =  1) 

p  P(T  =  1  I  0  =  1),  reliability  conditional  probability 

0  flood  indicator:  flood  does  not  exist  (0  =  0),  flood  occurs  (0  =  1) 

Part  4.  Selection  of  Optimal  Food  Warning  Threshold 

A,B,C  parameters  used  in  the  normal-linear  likelihood  model  (Equations  4.7-4. 9) 

Ci,C2  constants  governing  the  evolution  of  the  response  fraction  a 

D(h)  stage-damage  function  for  a  community 

Dqo  expected  loss  when  no  warning  is  given  and  no  flood  occurs  (zero) 

Doi  expected  community  property  loss  without  a  warning 

Dpi  expected  property  loss  without  a  warning  conditioned  on  forecast  s 

Dio  cost  of  evacuation  in  the  community 

Djo  expected  cost  of  community  evacuation  conditioned  on  forecast  s 

djo  cost  function  of  evacuation,  linear  function  of  response  fraction  a 

Dll  expected  community  property  loss  with  a  warning 
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D„ 

cl„ 

erf 

f(h  I  s) 
f(s  I  h) 
f, 

f2 

g(h) 

h 

k(s) 

MC 

MD 

MR(h  -  y) 
N 

N(/i,o) 

PflO 

Poi 

P.0 

p.. 

q(s,y) 
r  =  1,2,3 


expected  property  loss  with  a  warning  conditioned  on  forecast  s 

loss  function  with  a  warning 

standard  error  function 

posterior  distribution  of  h  given  a  forecast  s 

conditional  density  of  s  given  h 

sum  of  the  expected  property  loss  reductions  over  the  planning  horizon 
expected  property  loss  reduction  (difference  made  by  warning  system) 
expected  property  loss  reduction  at  stage  T 

objective  function  representing  credibility  of  forecast  system:  E{aN+i} 
prior  probability  density  of  flood  crest  h 
flood  crest 

marginal  probability  density  of  forcast  s 

maximum  evacuation  cost  with  a  full  response 

maximum  possible  damage  due  to  highest  flooding  with  no  response 

unit  reduction  function— reduction  of  MD  when  the  depth  of  flooding  is  (h  -  y)  and  a 

number  of  successive  flood  events  on  planning  horizon 

normal  distribution  with  mean  jj.  and  variance  s^ 
probability  of  a  correct  quiet 

probability  of  a  missed  forecast  (Type  I  error) 

probability  of  a  false  warning  (Type  II  error) 

probability  of  a  correct  warning 

probability  that  zone  of  elevation  y  will  be  flooded  conditioned  on  forecast  s 
structural  categories  in  the  floodplain 


1 
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s  forecasted  flood  crest 

s*  flood  warning  threshold;  warning  issued  when  s  ^  s* 

<{)  noninferior  decision  sequence  consisting  of  the  set  of  warning  thresholds  for  all  decision 

periods  in  the  planning  horizon 

y  elevation  of  a  zone  in  the  floodplain 

Ct  response  fraction  of  the  community  in  period  T 

6(h  -  y)  unit  damage  function  specifies  the  fraction  of  MD  when  flood  depth  is  (h  -  y) 

6Xz)  fraction  of  maximum  possible  damage  to  an  establishment  that  occurs  when  the  depth  of 

flooding  measured  from  the  first  floor  level  is  z 

standard  normal  distribution  function 

mean  and  standard  deviation  of  the  distribution  g(h) 
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